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Preface 



In every department of physical science there is only so much science, properly 
so-called, as there is mathematics. 

Immanuel Kant 

Most air defense systems in use or under development today, employ homing 
guidance to effect intercept of the target. By virtue of the use of onboard data 
gathering, the homing guidance system provides continually improving quality of 
target information right up to the intercept point. More than any single device, the 
guided missile has shaped the aerospace forces of the world today. Combat aircraft, 
for example, are fitted with airborne weapons that can be launched against enemy 
aircraft, ground forces, or strategic targets deep inside enemy territory. Also, the 
guided missile can be employed as a diversionary weapon to confuse ground and 
air forces. Ground-based missile systems have various range capabilities from a few 
miles to several thousand miles. These ground-based missiles are ballistic or nonbal- 
listic types, depending on their mission requirements. The design of a guided weapon 
(i.e., a missile) is a large undertaking, requiring the team effort of many engineers 
having expertise in the areas of aerodynamics, flight controls, structures, and propul- 
sion, among others. The different design groups must work together to produce the 
most efficient weapon in terms of high accuracy and low cost. 

The intent of this book is to present the fundamental concepts of guided 
missiles, both tactical, and strategic and the guidance, control, and instrumenta- 
tion needed to acquire a target. In essence, this book is about the mathematics of 
guided flight. This book differs from similar books on the subject in that it presents a 
detailed account of missile aerodynamic forces and moments, the missile mathemati- 
cal model, weapon delivery, GPS (global positioning system) and TERCOM (terrain 
contour matching) guidance, cruise missile mechanization equations, and a detailed 
analysis of ballistic guidance laws. Moreover, an attempt has been made to give 
each subject proper emphasis, while at the same time special effort has been put 
forth to obtain simplicity, both from the logical and pedagogical standpoint. Typi- 
cal examples are provided, where necessary, to illustrate the principles involved. 
Numerous figures give the maximum value of visual aids by showing important 
relations at a glance and motivating the various topics. Finally, this book will be 
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of benefit to engineers engaged in the design and development of guided missiles and 
to aeronautical engineering students, as well as serving as a convenient reference for 
researchers in weapon system design. 

The aerospace engineering field and its disciplines are undergoing a revolutionary 
change, albeit one that is difficult to secure great perspective on at the time of this 
writing. The author has done his best to present the state of the art in weapons systems. 
To this end, all criticism and suggestions for future improvement of the book are 
welcomed. 

The book consists of seven chapters and several appendices. Chapter 1 presents 
a historical background of past and present guided missile systems and the evolu- 
tion of modern weapons. Chapter 2 discusses the generalized missile equations of 
motion. Among the topics discussed are generalized coordinate systems, rigid body 
equations of motion, D’Alembert’s principle, and Lagrange’s equations for rotat- 
ing coordinate systems. Chapter 3 covers aerodynamic forces and coefficients. Of 
interest here is the extensive treatment of aerodynamic forces and moments, the vari- 
ous types of missile seekers and their function in the guidance loop, autopilots, and 
control surface actuators. Chapter 4 treats the important subject of the various types 
of tactical guidance laws and/or techniques. The types of guidance laws discussed 
in some detail are homing guidance, command guidance, proportional navigation, 
augmented proportional navigation, and guidance laws using modern control and 
estimation theory. Chapter 5 deals with weapon delivery systems and techniques. 
Here the reader will find many topics not found in similar books. Among the numer- 
ous topics treated are weapon delivery requirements, the navigation/weapon delivery 
system, the fire control computer, accuracies in weapon delivery, and modern topics 
such as situational awareness/situation assessment. Chapter 6 is devoted to strate- 
gic missiles, including the classical two-body problem and Lambert’s theorem, the 
spherical Earth hit equation, explicit and implicit guidance techniques, atmospheric 
reentry, and ballistic missile intercept. Chapter 7 focuses on cruise missile theory and 
design. Much of the material in this chapter centers on the concepts of cruise missile 
navigation, the terrain contour matching concept, and the global positioning system. 
Each chapter contains references for further research and study. Several appendices 
provide added useful information for the reader. Appendix A lists several fundamental 
constants, Appendix B presents a glossary of terms found in technical publications 
and books, Appendix C gives a list of acronyms, Appendix D discusses the standard 
atmosphere, Appendix E presents the missile classification. Appendix F lists past 
and present missile systems. Appendix G summarizes the properties of conics that 
are useful in understanding the material of Chapter 6, Appendix H is a list of radar 
frequencies, and Appendix I presents a list of the most commonly needed conversion 
factors. 

Such is the process of learning that it is never possible for anyone to say exactly 
how he acquired any given body of knowledge. My own knowledge was acquired 
from many people from academia, industry, and the government. Specifically, my 
knowledge in guided weapons and control systems was acquired and nurtured during 
my many years of association with the Department of the Air Force’s Aeronautical 
Systems Center, Wright-Patterson AFB, Ohio, while participating in the theory, 
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design, operation, and testing (i.e., from concept to fly-out) the air-launched cruise 
missile (ALCM), SRAM II, Minuteman III, the AIM-9 Sidewinder, and other programs 
too numerous to list. 

Obviously, as anyone who has attempted it knows, writing a book is hardly a soli- 
tary activity. In writing this book, I owe thanks and acknowledgment to various people. 
For obvious reasons, I cannot acknowledge my indebtedness to all these people, and so 
1 must necessarily limit my thanks to those who helped me directly in the preparation 
and checking of the material in this book. Therefore, I would like to acknowledge 
the advice and encouragement that I received from my good friend Dr. Guanrong 
Chen, formerly Professor of Electrical and Computer Engineering, University of 
Houston, Houston, Texas, and currently Chair Professor, Department of Electronic 
Engineering, City University of Hong Kong. In particular, I am thankful to Professor 
Chen for suggesting this book to Springer- Verlag New York and working hard to see 
that it received equitable consideration. Also, I would like to thank my good friend 
Dr. Victor A. Skormin, Professor, Department of Electrical Engineering, Thomas J. 
Watson School of Engineering and Applied Science, Binghamton University (SUNY), 
Binghamton, New York, for his encouragement in this effort. To Dr. Pravas R. 
Mahapatra, Professor, Department of Aerospace Engineering, Indian Institute of 
Science, Bangalore, India, I express my sincere thanks for his commitment and 
painstaking effort in reviewing Chapters 2-4. His criticism and suggestions have 
been of great service to me. Much care has been devoted to the writing and proof- 
reading of the book, but for any errors that remain I assume responsibility, and I will 
be grateful to hear of these. 

The author would like to express his appreciation to the editorial and production 
staff of Springer- Verlag New York, for their courteous cooperation in the production of 
this book and for the high standards of publishing, which they have set and maintained. 

Finally, but perhaps most importantly, I would like to thank my family for their 
forbearance, encouragement, and support in this endeavor. 



Dayton, Ohio 
November, 2003 



George M. Siouris 
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Introduction 



Rockets have been used as early as a.d. 1232, when the Chinese employed them as 
unguided missiles to repel the Mongol besiegers of the city of Pein-King (Peiping). 
Also, in the fifteenth century, Korea developed the sinkijon* (or Sin-Gi-Jeon ) rocket. 
Manufactured from the early fifteenth to mid-sixteenth century, the sinkijon was 
actively deployed in the northern frontiers, playing a pivotal role in fending off inva- 
sions on numerous occasions. Once out of the rocket launcher, the fire-arrows were 
set to detonate automatically near the target area. Also, the high-powered firearm was 
utilized in the southern provinces to thwart the Japanese marauders. The main body 
of the sinkijon ’ s rocket launcher was five to six meters long, the largest of its kind 
at that time** . A sinkijon was capable of firing as many as one hundred fire-arrows 
or explosive grenades. The fire-arrow contained a device equipped with gunpowder 
and shrapnel, timed to explode near the target. The introduction of gunpowder made 
possible the use of cannon and muskets that could fire projectiles great distances 
and with high velocities. It was desirable - in so far as the study of cannon fire is 
desirable - to learn the paths of these projectiles, their range, the heights they could 
reach, and the effect of muzzle velocity. Several years later, the sinkijon went through 
another significant upgrade, which enabled it to hurl a fire-arrow made up of small 
warheads and programmed to detonate and shower multiple explosions around the 
enemy. In 1451, King Munjong ordered a drastic upgrade of the hwacha (a rocket 
launcher on a cartwheel). This improvement allowed as many as one hundred sinki- 
jons to be mounted on the hwacha, boosting the overall firepower and mobility of the 
rocket. 

Since those early times and in one form or another, rockets have been used as 
weapons and machines of war, for amusement through their colorful aerial bursts, as 
life-saving equipment, and for communications or signals. The lack of suitable guid- 
ance and control systems may have accounted for the rocket’s slow improvement over 
the years. Strangely enough, it was the airplane rather than the rocket that stimulated 
the development of a guided missile as it is known today. 

* Sinkijon means “ghost-like arrow machine.” 

**The author would like to thank Dr. Jang Gyu Lee, Professor and Director of the Auto- 
matic Control Research Center, Seoul National University, Seoul, Korea, for providing the 
information on sinkijon. 
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In the twentieth century, the idea of using guided missiles came during World 
War I. Specifically, and as stated above, the use of the airplane as a military weapon 
gave rise to the idea of using remote-controlled aircraft to bomb targets. As early as 
1913, Rene Lorin, a French engineer, proposed and patented the idea for a ramjet 
powerplant. In 1924, funds were allocated in the United States to develop a missile 
using radio control. Many moderately successful flights were made during the 1920s 
with this control, but by 1932 the project was closed because of luck of funds. Radio- 
controlled target planes were the first airborne remote-controlled aircraft used by the 
Army and Navy. 

Dr. Robert H. Goddard was largely responsible for the interest in rockets back in 
the 1920s. Early in his experiments he found that solid-propellant rockets would not 
give him the high power or duration of power needed for a dependable supersonic 
motor capable of extreme altitudes. On March 16, 1926, Dr. Goddard successfully 
fired the first liquid-propellant rocket, which attained an altitude of 184 ft (56 m) and 
a speed of 60 mph (97 km/hr). Later, Dr. Goddard was the first to fire a rocket that 
reached a speed faster than the speed of sound. Moreover, he was the first to develop 
a gyroscopic steering apparatus for rockets, first to use vanes in the jet stream for 
rocket stabilization during the initial phase of a rocket in flight, and the first to patent 
the idea of step rockets. 

The first flight of a liquid-propellant rocket in Europe occurred in Germany 
on 14 March 1931. In 1932 Captain Walter Dornberger (later a general) of the 
German Army obtained the necessary approval to develop liquid-propellant rockets 
for military purposes [1]. Subsequently, by 1936 Germany decided to make research 
and development of guided missiles a major project, known as the “Peenemtinde 
Project,” at Peenemtinde, Germany. The German developments in the field of guided 
missiles during World War II were the most advanced of their time. Their most widely 
known missiles were the V-l and V-2 surface-to-air missiles (note that the designation 
VI and/or V2 is also found in the literature). As early as the spring of 1942, the original 
V-l had been developed and flight-tested at Peenemiinde. 

In essence, then, modern weapon (missile) guidance technology can be said 
to have originated during World War II in Germany with the development of the 
V-l and V-2 (German: A-4\ the A-4 stands for Aggregat-4, or fourth model in the 
development type series; the V stands for Vergeltungswajfe , or retaliation weapon, 
while some authors claim that initially, it stood for Versuchsmuster or experimental 
model) surface-to-surface missiles by a group of engineers and scientists at Peen- 
emiinde. It should be noted that static firing of rockets, notably the A-3, was per- 
formed as early as in the spring of 1936 at the Experimental Station, Kummersdorf 
West (about 17 miles south of Berlin). In the spring of 1942 the original V-l (also 
known by various names such as buzz bomb, robot bomb , flying bomb, air torpedo, 
or Fieseler Fi-103) had been developed and flight-tested at Peenemtinde. Thus, the 
V-l and V-2 ushered in a new type of warfare employing remote bombing by pilotless 
weapons launched over a hundred miles away through all kinds of weather, day and 
night [1], [3], 

The V-l was a small, midwing, pilotless monoplane, lacking ailerons but using 
conventional airframe and tail construction, having an overall length of 7.9 m (25.9 ft) 
and a wingspan of 5.3 m (17.3 ft). It weighed 2,180 kg (4,806 lb), including gasoline 
fuel and an 850 kg (1,874 lb) warhead. Powered by a pulsejet engine and launched 
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from an inclined concrete ramp 45.72 m (150 ft) long and 4.88 m (16 ft) above the 
ground at the highest end, the V-l flew a preset distance, and then switched on a release 
system, which deflected the elevators, diving the missile straight into the ground. The 
engine was capable of propelling the V-l 724 km/hr (450 mph). A speed of 322 km/hr 
(200 mph) had to be reached before the V-l propulsion unit could maintain the missile 
in flight. The range of the V-l was 370 km (230 miles). Guidance was accomplished 
by an autopilot along a preset path. Specifically, the plane’s (or missile’s) course 
stabilization was maintained by a magnetically controlled gyroscope that directed a 
tail rudder. When the predetermined distance was reached, as mentioned above, a 
servomechanism depressed the elevators, sending the plane into a steep dive. The V-l 
was not accurate, and it was susceptible to destruction by antiaircraft fire and aircraft. 
Several versions of the V-l were developed in Germany at that time. One version was 
designed for launch from the air. The missile could be carried under the left wing 
of a Heinkel He-111 aircraft. A manned V-l version was also developed, called the 
Reichenberg, flown first by Willy Fiedler, followed by Hanna Reitch. This version 
was planned for suicide missions. Three versions were built. 

The V-2 ( A-4 ) rocket was one of the most fearsome weapons of WWII. Successor 
to the V-l buzz bomb , the V-2 inflicted death, destruction, and psychological fear 
on the citizens of Great Britain. In essence, the V-2 was the first long-range rocket- 
propelled missile to be put into combat. Moreover, the V-2 was a liquid-propellant, 
14 m (45.9 ft) rocket that was developed between 1938 and 1942 under the tech- 
nical direction of Dr. Werner von Braun and Dr. Walter Dornberger, Commanding 
General of the Peenemiinde Rocket Research Institute. In addition to Great Britain, 
the V-2 was used to bomb other countries. However, although the first successful V-2 
test occurred on October 3, 1942, Adolf Hitler authorized full-scale development on 
July 27, 1943. The V-2 had movable vanes on the outer tips of its fins. These fins 
were used for guidance and control when the missile was in the atmosphere, which 
would be for most of its flight when used as a ballistic weapon. It also had movable 
solid carbon vanes projecting into the rocket blast for the same purpose when it was 
in rarified atmosphere. The first V-2, which landed in England in September 1944, 
was a supersonic rocket-propelled missile launched vertically and then automatically 
tilted to a 41°-47° angle a short time after launch. Furthermore, the V-2 had a liftoff 
weight of 12,873 kg (28,380 lb), developing a thrust of 27,125 kg (59,800 lb), a 
maximum acceleration of 6.4 g, reaching a maximum speed of about 5,705 km/h 
(3,545 mph), an effective range of about 354 km (220 miles), carrying a warhead of 
998 kg (2,201 lb). In addition, the powered flight lasted 70 sec, reaching a speed of 
about 6,000 ft/sec at burnout, with a burnout angle of about 45° measured from the 
horizontal. A flat-Earth model was assumed. Like the V-l, the V-2 was not known for 
its accuracy. For instance, the V-2 had a dispersion at the target of 10 miles (16 km) 
over a range of 200 miles (322 km). Active countermeasures against the V-2 were 
impossible at that time. Except for its initial programmed turn, it operated as a free 
projectile at extremely high velocity. The V-2 consisted of two main parts: (1) a 
directional reference made up of a gyroscopic assembly to control the attitude of the 
missile and a clock-driven pitch programmer, and (2) an integrating accelerometer in 
order to sense accelerations along the thrust axis of the missile, thereby determining 
velocity, and to cut off the engine upon reaching a predetermined velocity. In essence, 
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the V-2 system was the first primitive example of inertial guidance, making use of 
gyroscopes and accelerometers [3]. 

Several other German missiles were also highly developed during World War II 
and were in various stages of test. One of these, the Rheinbote (Rhein Messenger), was 
also a surface-to-surface missile. This rocket was a three-stage device with booster- 
assisted takeoff. Its range was 217 km (135 miles), with the third stage reaching over 
5,150 km/hr (3200 mph) in about 25 seconds after launch. The overall length of the 
rocket was about 1 1.3 m (37 ft). After having dropped a rearward section at the end 
of each of the first and second stages, it had a length of only 3.96 m (13 ft). The 
3.96 m (13 ft) section of the third stage carried a 40 kg (88 lb) high-explosive war- 
head. An antiaircraft or surface-to-air missile, the Wasserfall (Waterfall), was a remote 
radio-controlled supersonic rocket, similar to the V-2 in general principles of operation 
(e.g., both were launched vertically). When fully loaded, it had a weight of slightly 
less than 4,907 kg (5.4 tons). Its length was 7.62 m (25 ft). Designed for intercepting 
aircraft, the missile had specifications that called for a maximum altitude of 19,812 m 
(65,000 ft), a speed of 2,172 km/hr (1,350 mph), and a range of 48.3 km (30 miles). 
Its 90.7 kg (200 lb) warhead could be detonated by radio after the missile had been 
command-controlled to its target by radio signals. It also had an infrared proximity 
fuze and homing device for control on final approach to the target and for detonat- 
ing the warhead at the most advantageous point in the approach. Propulsion was 
to be obtained from a liquid-propellant power plant, with nitrogen-pressurized tanks. 
Another surface-to-air missile, the Schmetterling (Butterfly), designated IIS- 1 1 7, was 
still in the development stage at the close of the war. All metal in construction, it was 
3.96 m (13 ft) long and had a wingspan of 1.98 m (6.5 ft). Its effective range against 
low-altitude targets was 16 km (10 miles). It traveled at subsonic speed of about 
869 km/hr (540 mph) at altitudes up to 10,668 m (35,000 ft). A proximity fuze would 
set off its 24.95 kg (55 lb) warhead. Propulsion was obtained from a liquid-propellant 
rocket motor with additional help from two booster rockets during takeoff. Launching 
was to be accomplished from a platform, which could be inclined and rotated toward 
the target. The Schmetterling was developed at the Henschel Aircraft Works. 

The Enzian was another German surface-to-air missile (SAM). Designed to carry 
payloads of explosives up to 1000 pounds (453.6 kg), it was intended to be used against 
heavy-bomber formations. The Enzian was about 12 ft (3.657 m) long, had a wingspan 
of approximately 14 ft (4.267 m), and weighed a little over 2 tons (1,814.36 kg). 
Propelled by a liquid-propellant rocket, it was assisted during takeoff by four solid- 
propellant rocket boosters. The range of the Enzian was 16 miles (25.74 km), with 
a speed of 560 mph (901.21 km/hr), reaching an maximum altitude of 48,000 ft 
(14,630 m). In addition to the SAMs Germany had developed an air-to-air missile, 
designated the X-4. The X-4 was designed to be launched from fighter aircraft. Pro- 
pelled by a liquid-propellant rocket, it was stabilized by four fins placed symmetrically. 
Its length was about 6.5 ft (1.98 m) and span about 2.5 ft (0.762 m). Its range was 
slightly over 1.5 miles (2.414 km), and its speed was 560 mph (901.21 km/hr) at an 
altitude of 21,000 ft (6,401 m). Guidance was accomplished by electrical impulses 
transmitted through a pair of fine wires from the fighter aircraft. This missile was 
claimed to have been flown, but it was never used in combat. 
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The V-weapons, as mentioned earlier, were used to bombard London and 
southeastern England from launch sites near Calais, France, and the Netherlands. 
However, as the German armies were withdrawing from the Netherlands in March 
1945, the V-ls were launched from aircraft. Over 9,300 V-ls had been bred against 
England. By August 1944, approximately 1,500 V-ls had been shot down over 
England. Also, 4,300 V-2s had been launched in all, with about 1,500 against England 
and the remaining against Antwerp harbor and other targets. 

A project for developing missiles in the U.S.A. during World War II was started 
in 1941. In that year the Army Air Corps asked the National Defense Research Com- 
mittee to undertake a project for the development of a vertical, controllable bomb. 
The committee initiated a glide-bomb program, which resulted in standardization of 
a preset glide bomb attached to a 2,000 lb (907.2 kg) demolition bomb. The Azon, 
a vertical bomb controlled in azimuth only, went on the production line in 1943. 
Project Razon, a bomb controlled in both azimuth and range, was started in 1942. By 
1944, these glide bombs used remote television control. The Navy had a number of 
guided missile projects under development by the end of World War II. The Loon, a 
modification of the V-l, was to be used from ship to shore and to test guided-missile 
components. Another Navy missile, known as Gorgon IIC, used a ramjet engine with 
radar tracking and radio control. 

At the close of World War II the Americans obtained sufficient components to 
assemble two to three hundred V-2s from the underground factory, the Mittelwerk, near 
Nordhausen, Germany. The purpose of this was to use these V-2s as upper-atmosphere 
research vehicles carrying scientific experiments from JPL (Jet Propulsion Labora- 
tory), Johns Hopkins, and other organizations. 

In essence, the ballistic missile program in this country culminated with the 
development of the Atlas ICBM (intercontinental ballistic missile) (see Appendix F, 
Table F-l). In October 1953, and under a study contract from the U.S. Air Force, 
the Ramo-Woolridge Corporation (later Thomson-Ramo-Woolridge, or TRW) began 
work on a new ICBM. Within a year the program passed from top Air Force priority 
to top national priority. The first successful flight of a Series A Atlas ICBM took place 
on December 17, 1957, four months after the Soviet Union had announced that it had 
an ICBM. By the mid- 1959, more than eighty thousand engineers and technicians 
had participated in this program. 

Strictly speaking, missiles can be divided into two categories: (1) guided missiles 
(also called guided munitions), or tactical missiles, and (2) unguided missiles, or 
strategic missiles. Guided and unguided missiles can be defined as follows: 

Guided Missile: In the guided class of missiles belong the aerodynamic guided 
missiles. That is, those missiles that use aerodynamic lift to control its direction 
of flight. An aerodynamic guided missile can be defined as an aerospace vehicle, 
with varying guidance* capabilities, that is self-propelled through the atmosphere 
for the purpose of inflicting damage on a designated target. Stated another way, an 
aerodynamic guided missile is one that has a winged configuration and is usually 

* Guidance is defined here as the means by which a missile steers to, or is steered to, a target. 
In guided missiles, missile guidance occurs after launch. 
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fired in a direction approximately towards a designated target and subsequently 
receives steering commands from the ground guidance system (or its own, 
onboard guidance, system) to improve its accuracy. 

Guided missiles may either home to the target, or follow a nonhoming preset 
course . Homing missiles maybe active, semiactive, or passive. Nonhoming guided 
missiles are either inertially guided or preprogrammed [3]. (For more information, 
see Chapter 4.) 

Unguided Missiles: Unguided missiles, which includes ballistic missiles, follow the 
natural laws of motion under gravity to establish a ballistic trajectory. Examples of 
unguided missiles are Honest John, Little John, and many artillery-type rockets. 
Note that an unguided missile is usually called a rocket and is normally not a threat 
to airborne aircraft. (See also Chapter 6 for more details.) 

Typically, guided missiles are homing missiles, which include the following: (1) a 
propulsion system, (2) a warhead section, (3) a guidance system, and (4) one or more 
sensors (e.g., radar, sinfrared, electrooptical, lasers). Movable control surfaces are 
deflected by commands from the guidance system in order to direct the missile in 
flight; that is, the guidance system will place the missile on the proper trajectory to 
intercept the target. 

As stated above, homing guidance may be of the active, semiactive, or pas- 
sive type. Active guidance missiles are able to guide themselves independently after 
launch to the target. These missiles are of the so-called launch-and-leave class. For 
instance, air superiority fighters such as the F/A-22 Raptor that are designed with 
low-observable, advanced avionics and supercruise technologies are being developed 
to counter lethal threats posed by advanced surface-to-air missile systems (e.g., the 
U.S. HAWK MIM-23, Patriot MIM-104, Patriot Advanced Capability PAC-3, and the 
Russian SA-10 and SA-12 SAMs) and next-generation fighters equipped with launch- 
and-leave missiles. Therefore, an active guided missile carries the radiation source 
on board the missile. The radiation from the interceptor missile is radiated, strikes 
the target, and is reflected back to the missile. Thus, the missile guides itself on this 
reflected radiation. Consequently, a missile using active guidance will, as a rule, be 
heavier than semiactive or passive missiles. 

A semiactive missile uses a combination of active and passive guidance. A source 
of radiation is part of the system, but is not carried in the missile; that is, it is depen- 
dent on off-board equipment for guidance commands. More specifically, in semiactive 
missiles the source of radiation, which is usually at the launch point, radiates energy 
to the target, whereby the energy is reflected back to the missile. As a result, the mis- 
sile senses the reflected radiation and homes on it. A passive missile utilizes radiation 
originated by the target, or by some other source not part of the overall weapon system. 
Typically, this radiation is in the infrared region (e.g., Sidewinder-type missiles) 
or the visible region (e.g., Maverick), but may also occur in the microwave region 
(e.g., Shrike). Nonhoming guided missiles, as we shall presently discuss, are either 
inertially guided or preprogrammed. From the above discussion, we note that missile 
guidance can occur after launch. By guiding after launch, the effect of prelaunch aim- 
ing errors can be considerably minimized. Hence, the primary purpose of postlaunch 
guidance is to relax prelaunch aiming requirements. 
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Two common types of missiles that pose a threat to aircraft are the air-to-air 
{ AA ), or air- intercept, missile {AIM), and the surface-to-air missile {SAM) mentioned 
earlier. The AA and SAM missiles belong to the tactical and defense missile class, and 
are launched from interceptor fighter aircraft, employing various guidance techniques. 
Surface-to-air missiles can be launched from land- or sea-based platforms. They 
too have varying guidance and propulsion capabilities that influence their launch 
envelopes relative to the target. Furthermore, these missiles employ sophisticated 
electronic countermeasure (ECM) schemes to enhance their effectiveness. It should 
be pointed out that since weight is not much of a problem, these missiles are often 
larger than their air-to-air counterparts, and they can have larger warheads and longer 
ranges. 

In attempting to intercept a moving target with a missile, a desired trajectory will 
be needed in which the missile velocity leads the line of sight {LOS) by the proper 
angle so that for a constant-velocity target the missile flies a straight-line path to 
collision. In homing systems, for example, the target tracker is in the missile, and 
in such a case it is the relative movement of target and missile that is relevant. The 
two-dimensional end-game geometry of an ideal collision course will be discussed 
later in this book. Typically, an aerodynamic missile is controlled by an autopilot, 
which receives lateral acceleration commands from the guidance system and causes 
aerodynamic surfaces to move so as to attain these commanded accelerations. Since 
in general, there are two lateral missile coordinate axes, the general three-dimensional 
attack geometry can be resolved into these two directions. 

Ballistic missiles belong to the strategic missile class, and are characterized by 
their trajectory. A ballistic missile trajectory is composed of three parts (for more 
details, see Chapter 6). These are (1) the powered flight portion, which lasts from 
launch to thrust cutoff (or burnout); (2) the free -flight portion, which constitutes most 
of the trajectory, and (3) the reentry portion, which begins at some point (not defined 
precisely) where the atmospheric drag becomes a significant force in determining the 
missile’s path and lasts until impact on the surface of the Earth (i.e., a target). Typically, 
ballistic missiles rely on one or more boosters and an initial steering vector. Once in 
flight, they maintain this vector with the aid of gyroscopes. Therefore, a ballistic 
missile may be defined as a missile that is guided during the powered portion of the 
flight by deflecting the thrust vector, becoming a free-falling body after engine cutoff. 
However, as already noted, in ballistic missiles part of the guidance occurs before 
launch. Hence, prelaunch errors translate directly into miss distance. One important 
feature of these missiles is that they are roll stabilized, resulting in simplification of 
the analysis, since there is no coupling between the longitudinal and the lateral modes. 
Ballistic missiles are the type least likely to be intercepted. A ballistic missile can 
have surprising accuracy. Ballistic missiles can be classified according to their range. 
That is, short range (e.g., up to 300 nm (nautical miles) or 556 km), intermediate range 
(e.g., 2500 nm or 4632.5 km), and long range (over 2500 nm or 4632.5 km). Examples 
of these classes are as follows: (1) short range - Pershing , Sergeant, and Hawk class; 
(2) intermediate range - Thor, Jupiter, and Polaris/Poseidon/Trident, and (3) long 
range - Minuteman I-III, the MX, and Titan missiles. Note that ballistic missiles 
capable of attaining very long ranges (e.g., over 5000 nm) or intercontinental range, 
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are given the ICBM designator [2], [4]. Recently, the U.S. Air Force formulated plans 
for a new ICBM, likely to be named Minuteman TV. A possible start development 
date is for the year(s) 2004-2005. Among the enhancements being examined are 
communications upgrades, an additional postboost vehicle that could maneuver the 
warhead after separation from the missile, and a new rocket motor. 

In common use today are the following abbreviations, which use the term ballistic 
missile in the sense that the type of missile and its capacity are indicated (for a detailed 
list of acronyms, see Appendix C): 

IRBM : Intermediate Range Ballistic Missile 
ICBM : Intercontinental Ballistic Missile 
AICBM : Anti-Intercontinental Ballistic Missile 

SLBM : Submarine-Launched Ballistic Missile (or FBM - Fleet Ballistic Missile) 
ALBM : Air-Launched Ballistic Missile 
MMRBM : Mobile Mid-Range Ballistic Missile. 

The range has much to do with using this kind of missile designator, which like the 
point-to-point designator, is used with the vehicle’s popular name. It should be noted 
at this point that essentially, the difference between the ballistic and aerodynamic 
missiles lies in the fact that the former does not rely upon aerodynamic surfaces to 
produce lift and consequently follows a ballistic trajectory when thrust is terminated. 
Aerodynamic missiles, as stated earlier, have a winged configuration. 

Ballistic missiles use inertial guidance, sometimes aided with star trackers and/or 
with the Global Positioning System (GPS). More specifically, inertial guidance is used 
for a ballistic trajectory only during the very early part of the flight (i.e., up to fuel cut- 
off) in order to establish proper velocity for a hit by free fall. In ballistic missiles, the 
intent is to hit a given map reference, as opposed to aerodynamic missiles, whose intent 
is to intercept a moving and at times highly maneuverable target. Long-range inter- 
continental ballistic missiles are categorized as surface-to-surface. As stated above, 
ballistic missiles use inertial guidance to hit a target. The modern inertial naviga- 
tion and guidance system is the only self-contained single source of all navigation 
data. Self-contained inertial navigation depends on the integration of acceleration 
with respect to a Newtonian reference frame. That is, inertial navigation depends on 
integration of acceleration to obtain velocity and position. The inertial navigation 
system (INS) provides a reliable all-weather, worldwide navigation capability that is 
independent of ground-based navigation aids. The system develops navigational data 
from self-contained inertial sensors (i.e., gyroscopes and accelerometers), consisting 
of a vertical accelerometer, two horizontal accelerometers, and three single-degree- 
of-freedom gyroscopes (or 2 two-degree-of-freedom gyroscopes). In addition to the 
conventional mechanical gyroscopes, there is a new generation of inertial sensors such 
as the RLG (Ring Laser Gyro), the FOG (Fiber-Optic Gyro), and the MEMS (Micro 
Electro-Mechanical Sensor), which functions as both a gyro and an accelerometer. 
Note that the MEMS devices are fundamentally different from the RLG and FOG opti- 
cal sensors. The design of MEMS allows a single chip to function as both a gyro and an 
accelerometer. The sensing elements are mounted in a four-gimbal, gyro-stabilized 
inertial platform. The accelerometers are the primary source of information. They 
are maintained in a known reference frame by the gyroscopes. That is, the precision 
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gyro-stabilized platform is used for reference. Attitude and heading information is 
obtained from synchro devices mounted between the platform gimbals. Therefore, 
the heart of the inertial navigation system is the inertial platform. The platform has 
four gimbals for all-attitude operation, with the outermost gimbal being the outer roll, 
which has unlimited freedom. Proceeding inward, the next gimbal is pitch, which is 
normally limited to ±105° of freedom. The next inward gimbal is inner roll, which is 
redundant with the outer roll axis but is required in order to eliminate what is called 
gimbal lock and is limited to ±15° angular freedom. All inertial sensors are mounted 
on the azimuth gimbal, the innermost gimbal. The gyroscopes are mounted such that 
the vertical gyroscope is mounted with its spin axis parallel to the azimuth gimbal 
rotational axis and positioned to coincide with the local vertical when the platform 
is erected to X and Y (level) accelerometer nulls. The X and Y axis accelerome- 
ters, mounted on the azimuth structure, are aligned to sense horizontal accelerations 
along the gyro X and Y axes, respectively, while the Z, or vertical, accelerometer 
senses accelerations along the azimuth axis. After being supplied with initial position 
information, the INS is capable of continuously updating extremely accurate displays 
of position, ground speed, attitude, and heading. In addition, it provides guidance or 
steering information for autopilot and flight instruments (in the case of aircraft). 

Note that the above discussion was for gimbaled inertial navigation systems. There 
is also a class of strapdown INSs in which the inertial sensors are mounted directly on 
the host vehicle frame. In this way, the gimbal structure is eliminated. In the strapdown 
version of the INS, wherein sensors are mounted directly on the vehicle, the transfor- 
mation from the sensor to inertial reference is “computed” rather than mechanized. 
Specifically, the strapdown system differs from the gimbaled system in that the specific 
force is measured in the body frame, and the attitude transformation to the naviga- 
tion specific force is computed from the gyro data, because the strapdown sensors are 
fixed to the vehicle frame. Regardless of mechanization (i.e., gimbaled or strapdown), 
alignment of an inertial navigation system is of paramount importance. In alignment, 
the accelerometers must be leveled (i.e,, indicating zero output), and the platform 
must be oriented to true north. This process is normally called gyrocompassing. 

In ballistic missiles (in particular ICBMs ), rocket propulsion is employed to 
accelerate the missile to a position of high altitude and speed. This places it on a 
trajectory that meets certain guidance specifications in order to carry a warhead, or 
other payload, to a preselected target. An operational ballistic missile may acquire 
speeds up to 15,000 mph (24,140 km/hr) or better at heights of several hundred miles. 
After boost burnout ( BBO ), or engine shutoff, the missile payload travels along a 
free-fall trajectory to its destination; its motion follows, approximately, the laws of 
Keplerian motion. A special type of onboard navigation/guidance computer is used 
in ballistic missiles in which the platform (e.g., in gimbaled systems) maintains its 
alignment in space for the few minutes during which the inertial system is operating 
to launch the warhead. The computer is fed the velocity and position that the warhead 
ought to achieve when the motors are cut off. Consequently, the actual positions and 
velocities are recorded from the information taken from the inertial platform, and by 
comparing the two, a correction may be passed to the control system of the missile. 
Thus, the correction ensures that the motors are cut off when the warhead is traveling 
at a velocity and from a position that will enable it to hit the same target as if it had 
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followed exactly a planned (or programmed) flight path or trajectory. The planned 
path takes into account the change of gravity due to the forward movement of the 
missile, the change in the force of gravity due to upward movement of the missile, and 
the Earth’s tilt, rotation, and Coriolis acceleration. However, the planned path may 
involve a good deal of calculation, and as a result it may not be easy to alter the aiming 
point by more than a small amount without a completely new plan. It was mentioned 
earlier that part of the guidance of a ballistic missile occurs before launch. Moreover, 
during the powered portion of the flight, the objective of the guidance system is to 
place the missile on a trajectory with flight conditions that are appropriate for the 
desired target. This is equivalent to steering the missile to a burn-out point that is 
uniquely related to the velocity and flight-path angle for the specified target range. 

Another type of strategic missile is the now canceled USAF’s SRAM II missile. 
The SRAM (Short-Range Attack Missile) II was a standoff. air-launched, inertially 
guided strategic missile. As designed, the missile had the capability to cover a large 
target accessibility footprint when launched with a wide range of initial conditions. 
The missile was designed to be powered by a two-pulse solid-fuel rocket motor 
with a variable intervening coast time. The guidance algorithm was based on modern 
control linear quadratic regulator ( LQR ) theory, with the current missile state (a vector 
consisting of position, velocity, and other parameters) provided by a strapdown inertial 
navigation system. The SRAM II trajectory was dependent on the relative locations of 
the launch point and target, as well as the flight envelope characteristics of the carrier 
(i.e., aircraft). 

Still another class of strategic missiles is the nuclear ALCM (Air-Launched Cruise 
Missile) designated as AGM-86B. The ALCM uses an inertial navigation system 
together with terrain contour matching ( TERCOM ) for its guidance. A later version 
of the ALCM , known as the CALCM (Conventionally Armed Air-Launched Cruise 
Missile) and designated AGM-86C, uses an INS integrated with the GPS and/or 
TERCOM (for more information, see Chapter 7). 

It should be pointed out that there is still another class of missiles, namely, radia- 
tion missiles. In radiation missiles, radiation energy is transmitted as either particles 
or waves through space at the speed of light. Radiation is capable of inflicting damage 
when it is transmitted toward the target either in a continuous beam or as one or more 
high-intensity, short-duration pulses. Weapons utilizing radiation are referred to as 
directed high-energy weapons ( DHEW ). These are as follows: 

1. Coherent Electromagnetic Flux: The coherent electromagnetic flux is produced 
by a high-energy laser ( HEL ). The HEL generates and focuses electromagnetic 
energy into an intense concentration or beam of coherent waves that is pointed at 
the target. This beam of energy is then held on the target until the absorbed energy 
causes sufficient damage to the target, resulting in eventual destruction. On the 
other hand, radiation from a laser that is delivered in a very short period of time 
with a high intensity is referred to as a pulse-laser beam. (For more details on 
high-energy weapons see Section 6.9.) 

2. Noncoherent Electromagnetic Pulse ( EMP ): The noncoherent electromagnetic 
pulse consists of an intense electronic signal of very short duration that travels 
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through space just as a radio signal does. When an EMP strikes an aircraft, the 
electronic devices in the aircraft can be totally disabled or destroyed. 

3. Charged Nuclear Particles: The charged-particle-beam weapon is the newest of 
the developing threats that utilizes radiation in the form of accelerated subatomic 
particles. These particles, or bunches of particles, may be focused on the target 
by means of magnetic fields. Thus, considerable damage can result. This type of 
weapon has the advantage that it will propagate through visible moisture, which 
tends to absorb energy generated by the HEL. 

Regardless of the type of missile, a development cycle must be formulated that takes 
into account several phases of design and analysis. The missile development cycle 
commences with concept formulation, where one or more guidance methods are pos- 
tulated and examined for feasibility and compatibility with the total system objectives 
and constraints. Surviving candidates are then compared quantitatively, and a baseline 
concept is adopted. Specific subsystem and component requirements are generated 
via extensive tradeoff and parametric studies. Factors such as missile capability (e.g., 
acceleration and response time), sensor function (e.g., tracking, illumination), 
accuracy (signal to noise, waveforms), and weapons control (e.g., fire control logic, 
guidance software) are established by means of both analytical and simulation tech- 
niques. After iteration of the concept/requirements phase and attainment of a set of 
feasible system requirements, the analytical design is initiated. During this stage, the 
guidance law is refined and detailed, a missile autopilot and the accompanying con- 
trol actuator are designed, and an onboard sensor tracking and stabilization system is 
devised. This design phase entails the extensive use of feedback control theory and the 
analysis of nonlinear, nonstationary dynamic systems subjected to deterministic and 
random inputs. Finally, determination of the sources of error and their propagation 
through the system are of fundamental importance in setting design specifications 
and achieving a well-balanced design. 

From the above discussion, one can safely say that of vital interest in missile 
design is the development of advanced guidance and control concepts. For example, 
in the design of a guidance law for a homing missile, a continued effort should be 
the study of homing guidance and the means to optimize its performance in various 
intercept situations. The classical approach to missile guidance involves the use of a 
low-pass filter for estimating the line-of-sight angular rate along with a proportional 
guidance law. In addition to the classical methods, we will discuss the use of opti- 
mized digital guidance and control laws for highly dynamic engagements associated 
with air-to-air missiles, where the classical approaches often fail to achieve accept- 
able performance. Conventional proportional navigation systems, as will be discussed 
later in this book, have been improved with time-variable filtering, and the design pro- 
cess has been refined with automatic computer methods. Advanced guidance systems 
having superior performance have been designed with on-line Kalman estimation for 
filtering noisy radar data and with optimal control gains expressed in closed form. For 
instance, trajectory estimators are designed routinely using Kalman filtering theory 
and provide minimum variance estimates of key guidance variables based upon a 
linearized model of the trajectory. The guidance laws are commonly designed to 
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yield as small a miss distance as possible, consistent, of course, with the missile’s 
acceleration capability. This is accomplished by mathematically requiring the com- 
manded acceleration to minimize an appropriate performance index (or cost function) 
involving both the miss distance and the missile acceleration level. Today, the concept 
of optimized guidance laws is well understood in applications where information con- 
cerning the target range and line-of-sight angle is available. This is the case when the 
homing sensor is an active or semiactive radar (RF) or laser range finder. Moreover, 
considerable attention has been given to developing advanced guidance concepts for 
the situation in which direct measurements of range are unavailable, as with passive 
infrared or electro-optical sensors. 

Synthesis of sample data homing and command guidance systems is also of par- 
ticular importance, as will be discussed later. Classical servo theory has been used to 
design both hydraulic and electric seeker servos that are compatible with requirements 
for gyro-stabilization and fast response. Furthermore, pitch, yaw, and roll autopilots 
have been designed to meet such problems as Mach variation, altitude variation, 
induced roll moments, instrument lags, body-bending modes, guidance response, and 
guidance stability. Although classical theory is still applicable to autopilots, research 
efforts are continually made to apply modern control theory to conventional autopilot 
design and adaptive autopilot design. 

Optimal control and estimation theory is commonly used in the design of advanced 
guidance systems. Specifically, since the late 1960s and early 1970s, considerable 
research has been devoted to applying modern optimal control and estimation theory 
in the development of optimized advanced tactical and strategic missile guidance sys- 
tems. In particular, this technology has been used to develop tracking algorithms that 
extract the maximum amount of information about a target trajectory from homing 
sensor data and to derive guidance and control laws that optimize the use of this infor- 
mation in directing the missile toward the selected target. Performance improvements 
attainable with optimized systems over conventional guidance and control techniques 
are most significant against airborne maneuverable targets, where target acceleration 
information and rapid guidance system response time are required to achieve accept- 
able accuracy, in minimum time. Historically, surface-to-air missiles were among the 
first missiles to implement digital guidance systems. Such missiles may employ com- 
mand guidance whereby all digital computation is done on the ground with guidance 
commands telemetered to the missile. Today, the ease of availability of microproces- 
sors makes digital processing increasingly attractive for small, lightweight air-to-air 
missiles. Recently developed neural network algorithms and fuzzy logic theory serve 
as possible approaches to solving highly nonlinear flight control problems. Thus, the 
use of fuzzy logic control is motivated by the need to deal with nonlinear flight control 
and performance robustness problems. 

It was noted earlier that prior to beginning an engineering development program 
for a digital guidance and control system, it is desirable to perform a detailed computer- 
aided feasibility study within the context of a realistic missile-target engagement 
model. In order to accomplish these, guidance and control laws that have been 
developed and evaluated for simplified missile-target engagement scenarios must 
be extended and adapted to the air-to-air missile situation and then implemented in a 
complete three-dimensional engagement model. 
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Finally, microprocessor technology will allow future application of more 
sophisticated guidance and control laws that consider the effects of uncertain system 
parameters than have heretofore been considered for tactical missiles. System minia- 
turization is becoming more and more common in weapon systems. For example, a 
miniaturized system that can integrate GPS and inertial guidance to increase accuracy 
of Army and Navy artillery shells has already been developed. These systems can be 
placed on a circuit board and are small enough to fit into the nose of an artillery shell. 
Above all, a single processor placed on the board can be used to handle GPS and iner- 
tial data from MEMS. The Army’s XM-982 and the Navy’s Extended Range Guided 
Munition ( ERGM ) will use the GPS system (see also Appendix F). Missile guidance 
systems are advancing on several fronts as GPS spreads into old and new systems, 
automatic target recognition moves toward deployment, and ballistic missile defense 
programs improve the state of the art in data fusion and infrared sensors. Missile 
systems presently under research and development will evolve into smaller, more 
accurate missiles. 

A revolutionary new generation of miniature loitering smart weapons (or sub- 
munition) is the U.S. Air Force’s LOCAAS (Low-Cost Autonomous Attack System) 
missile that was designed and flight-tested in the 1990s as a gliding weapon for 
armored targets only. LOCAAS can be air launched singly or in a self-synchronizing 
swarm that will deconflict targets so only one LOCAAS pursues each target. This 
futuristic smart weapon has a mind of its own. Scanning the land below, these weapons 
can identify and destroy mobile launchers. The key here is that they can distinguish 
between different targets and then shape their warheads to inflict maximum damage. 
Nose to tail, these $40,000, 31-inch (0.787 meter) long air-to-surface weapons will 
be anything but small in performance. The current production version calls for a five- 
pound turbojet engine with thirty pounds of thrust to fly 100 m/sec (328 ft/sec) while 
hunting for fast-moving missile launchers over a large target area. The size of a soup 
bowl, the warhead uses a shaped charge to transform a copper plate into fragments, 
a shuttlecock-shaped slug, or a rod that can penetrate several inches of high-carbon 
steel. That is, its warhead can explode into fragments, a long-rod penetrator, or a 
slug, depending on the type of target it detects. Without designating a specific target, 
flight crews will leave the thinking to the missile’s three-dimensional imaging ladar 
(or laser radar) and use its target recognition system in its nose to continuously scan 
target areas. That is, the LOCAAS seeker uses advanced target recognition algorithms 
to detect, prioritize, reject, and select targets. As many as two hundred of these flying 
smart weapons can be swooping down on an enemy battlefield. 
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The Generalized Missile Equations of Motion 



2.1 Coordinate Systems 

2.1.1 Transformation Properties of Vectors 

In a rectangular system of coordinates, a vector can be completely specified by 
its components. These components depend, of course, upon the orientation of the 
coordinate system, and the same vector may be described by many different triplets 
of components, each of which refers to a particular system of axes. The three 
components that represent a vector in one set of axes, will be related to the com- 
ponents along another set of axes, as are the coordinates of a point in the two 
systems. In fact, the components of a vector may be regarded as the coordinates 
of the end of the vector drawn from the origin. This fact is expressed by saying 
that the scalar components of a vector transform as do the coordinates of a point. 
It is possible to concentrate attention entirely on the three components of a vector 
and to ignore its geometrical aspect. A vector would then be defined as a set of 
three numbers that transform as do the coordinates of a point when the system of 
axes is rotated. It is often convenient to designate the coordinate axes by numbers 
instead of letters x, y. z so that the components of a vector will be ci\ , ay, and ay. 
The designation for the whole vector is a,-, where it is understood that the sub- 
script i can take on the value 1, 2, or 3. A vector equation is then written in the 
form 



ai=bi. (2.1) 

This represents three equations, one for each value of the subscript i . The rotation 
of a system of coordinates about the origin may be represented by nine quantities 
Yiji, where ]/ f // is the cosine of the angle between the i-axis in one position of the 
coordinates and the /-axis in the other position. These nine quantities give the angles 
made by each of the axes in one position with each of the axes in the other. They are 
also the coefficients in the expression for the transformation of the coordinates of a 
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point. The cosines can be conveniently kept in order by writing them in the form of 
a matrix: 



Ynj' Y 12 ' Y 13' 
Y2V 722' 723' 
73 T 732' 733' 



( 2 . 2 ) 



Of the nine quantities, only three are independent, since there are six independent 
relations between them. Since Yij' can be considered as the component along the 
j '- axis in one coordinate system of a unit vector along the /-axis in the other, then 



Yfy + Ynj + 7,-y - £ y ij' ~ 1 • (2.3a) 

r 



This will be true for every value of / . Similarly, 



£4 = L ») 



The components of a vector, or the coordinates of a point, can be transformed from 
one system of coordinates to the other by 



at — Yivay + Yii'aj + Yiyay = Yij' a y. ( 2 . 4 ) 

Here a ,' represents the components of the vector a in one system of coordinates, and 
a i the components in the other. The summation sign is omitted in the last term, since 
it is to be understood that a sum is to be carried out over all three values of any index 
that is repeated. 

2.1.2 Linear Vector Functions 

If a vector is a function of a single scalar variable, such as time, each component 
of the vector is independently a function of this variable. If the vector is a linear 
function of time, then each component is proportional to the time. A vector may also 
be a function of another vector. In general, this implies that each component of the 
function depends on each component of the independent vector. Moreover, a vector 
is a linear function of another vector if each component of the first is a linear function 
of the three components of the second. This requires nine independent coefficients of 
proportionality. The statement that a is a linear function of b means that 

a i = C n bi + C12&2 + C13/73, 

«2 = C2I&I + C22&2 + C 23 ^ 3 > ( 2 . 5 ) 

«3 = C31&1 + C32&2 + <233/73. 

Using the summation convention as in ( 2 . 4 ), this becomes 



at — Cijbj. 



(2.6) 
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A relationship such as that in (2.6) must be independent of the coordinate system 
in spite of the fact that the notation is clearly based on specific coordinates. The com- 
ponents a i and b, are with reference to a particular coordinate system. The constants 
Cjj also have reference to specific axes, but they must so transform with a rotation of 
axes that a given vector b always leads to the same vector a. 

If the coordinate system is rotated about the origin, the vector components will 
change so that 



a i = Yij'a y = Ci j Y jk ’ h’ ■ (2.7) 

If both sides of this equation are multiplied by yi’i and the equations for the three 
values of i are added, the result is 



Yl’i Yij' a j 1 = a l’ = (Yl’iQjYjk')bk'- 



(2.8) 



If the quantity Yl’iCijYjk’ * s called Cpk’, then 



av = Ci'k'bk'. (2.9) 

This relationship between the components in this system of coordinates is the 
same vector relationship as was expressed by the Ca in the original system of 
coordinates. 



2.1.3 Tensors 

Tensor is a general name given to quantities that transform in prescribed ways when 
the coordinate system is rotated. A scalar is a tensor of rank 0, for it is independent 
of the coordinate system. A vector is a tensor of rank 1 . Its components transform as 
do the coordinates of a point. A tensor of rank 2 has components that transform as do 
the quantities Cij . Put another way, a scalar is a quantity whose specification (in any 
coordinate system) requires just one number. On the other hand, a vector (originally 
defined as a directed line segment) is a quantity whose specification requires three 
numbers, namely, its components with respect to some basis. In essence, scalars and 
vectors are both special cases of a more general object called a tensor of order n, 
whose specification in any given coordinate system require 3" numbers, again called 
the components of the tensor. In fact, 

scalars are tensors of order 0, with 3° = 1 components, 
vectors are tensors of order 1, with 3 1 = 3 components. 

Tensors can be added or subtracted by adding or subtracting their corresponding 
components. They can also be multiplied in various ways by multiplying components 
in various combinations. These and other possible operations with tensors will not be 
described here. 

A tensor of the second rank is said to be symmetric if Cij = Cj ; and to be antisym- 
metric if Cij = —Cjj. An antisymmetric tensor has its diagonal components equal to 
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zero. Any tensor may be regarded as the sum of a symmetric and an antisymmetric 
part for 

Cij = \[Cij + Cji] + \[Cij - Cji] (2.10a) 

and 

2 [Cij + Cji ] - Sij \ [ Cij - Cji ] = Aij , (2.10b) 

where Sij is symmetric and A// is antisymmetric. Numerous physical quantities have 
the properties of tensors of the second rank, so that the inertial properties of a rigid body 
can be described by the symmetric tensor of inertia. By way of illustration, consider 
that we are given two vectors A and B . There are nine products of a component of A 
with a component of B. Thus, 



AiBi(i,k= 1,2,3). 

Suppose we transform to a new coordinate system K ' , in which A and B have compo- 
nents A' and B'. . Then the transformation of a coordinate system can be expressed as 

A/ = oijfAk, 

where A^, A( are the components of the vector in the old and new coordinate systems 
K and K' , respectively, and a,-^ is the cosine of the angle between the ;th axis of K' 
and the kth axis of K. Thus, 



Aj — (%i'k Af , B/. — Otk'm , 



and hence 



Aj Bj. — Oti'lCtlc'm A / B m . 



Therefore, A, Bj is a second-order tensor. 



2.1.4 Coordinate Transformations 

There are three commonly used methods of expressing the orientation of one three- 
axis coordinate system with respect to another. The three methods are ( 1 ) Euler angles, 
(2) direction cosines, and (3) quaternions. The Euler angle method, which is the con- 
ventional designation relating a moving-axis system to a fixed-axis system, is used 
frequently in missile and aircraft mechanizations and/or simulations. The common 
designations of the Euler angles are roll (</>), pitch (9), and yaw (\[r). Its strengths lie 
in a relatively simple mechanization in digital computer simulation of vehicle (i.e., 
missile or aircraft) dynamics. Another beneficial aspect of this technique is that the 
Euler angle rates and the Euler angles have an easily interpreted physical signifi- 
cance. The negative attribute to the Euler angle coordinate transformation method is 
the mathematical singularity that exists when the pitch angle 9 approaches 90°. The 
direction cosine method yields the direction cosine matrix (DCM), which defines the 
transformation between a fixed frame, say frame a, and a rotating frame, say frame b. 
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such as the vehicle body axes. Specifically, the DCM is an array of direction cosines 
expressed in the form 



c 



C] 1 C12 C 13 
C21 C22 C23 
C31 C32 C33 



where Cjk is the direction cosine between the / th axis in the a frame and the fcth axis 
in the b frame. Since each axis system has three unit vectors, there are nine direction 
cosines. Direction cosines have the advantage of being free of any singularities such 
as arise in the Euler angle formulation at 90° pitch angle. The main disadvantage 
of this method is the number of equations that must be solved due to the constraint 
equations. (Note that by constraint equations we mean cn = C22C33 — C23C32, C21 = 
C13C32 -C12C33, etc.) 

In order to resolve the ambiguity resulting from the singularity in the Euler angle 
representation of rotations about the three axes, a four-parameter system was first 
developed by Euler in 1776. Subsequently, Hamilton modified it in 1843, and he 
named this system the quaternion system. Therefore, a quaternion [Q] is a quadruple 
of real numbers, which can be written as a three-dimensional vector. Hamilton adopted 
a vector notation in the form 



[£?] = qo + i«7i +j <?2 + kc /3 = ( qo , 7t> 72, 73 ) = (qo, q), (2.1 1) 

where 70 , q 1 , 72 , 73 are real numbers and the set { i, j , k ) forms a basis for a quaternion 
vector space. From the orthogonality property of quaternions, we have 

7o + 7f + 72 + 73 = 1- (2.12) 

In terms of the Euler angles <p, 0, \j/, we have 

70 = cos(Vf/ 2 ) cos( 0 / 2 ) cos( 0 / 2 ) — sin(i/f/ 2 ) sin( 0 / 2 ) sin(</>/ 2 ), 
q l — sin(d/ 2 ) sin((/)/ 2 ) cos(i/f/ 2 ) + sin(i/f/ 2 ) cos( 0 / 2 ) cos(</>/ 2 ), 
q 2 — sin( 0 / 2 ) cos(i/f/ 2 ) cos(</>/ 2 ) — sin(i/r/ 2 ) sin((/>/ 2 ) cos( 0 / 2 ), 

<73 = sin(</>/ 2 ) cos(i/r/ 2 ) cos( 0 / 2 ) + sin(i/r/ 2 ) sin( 0 / 2 ) cos(</>/ 2 ). 



Suppose now that we wish to transform any vector, say V, from body coordinates 
Y h into the navigational coordinates V". This transformation can be expressed as 
follows: 

\ n = C n b V 6 , 

where C'l is the direction cosine matrix, or equivalently, using quaternions, 

Y" = qV h q*, 



where q* is the conjugate of q. Then [7] 



C n — 
L b ~ 



'%■ 



- qi - q \- q \ 2(7172-7073) 

2(7172 + 7073) ql ~ q \ 

2(7173-7072) 2(7273 + 7071) 75 - q \ 

For more details on the quaternion and its properties, the reader is referred to [7], 



■ q ] 



■qj 



2(7173- 

2(7273 



- 7072) 

-7071) 



7o ' 



■7 2 2 - 



■73. 
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The coordinate system that will be adopted in the present discussion is a 
right-handed system with the positive .r-axis along the missile’s longitudinal axis, the 
y-axis positive to the right (or aircraft right wing), and the 4 -axis positive down (i.e., the 
Z-axis is defined by the cross product of the x- and y-axis). This coordinate system 
is also known as north-east-down (NED) in reference to the inertial north-east-down 
sign convention [5], [7]. It should be noted here that the coordinate system used in the 
present development is the same one used in aircraft. Four orthogonal-axes systems 
are usually defined to develop the appropriate equations of vehicle (aircraft or missile) 
motion. They are as follows: 

1 . The inertial frame, which is fixed in space, and for which Newton’s Laws of Motion 
are valid. 

2. An Earth-centered frame that rotates with the Earth. 

3. An Earth-surface frame that is parallel to the Earth’s surface, and whose origin is 
at the vehicle’s center of gravity (eg) defined in north, east, and down directions. 

4. The conventional body axes are selected to represent the vehicle. The center of 
this frame is at the eg of the vehicle, and its components are forward, out of the 
right wing, and down. 

In ballistic missiles, two other common coordinate systems are used. These coordinate 
systems are 

1. Launch Centered Inertial: This system is inertially fixed and is centered at launch 
site at the instant of launch. In this system, the x-axis is commonly taken to be in 
the horizontal plane and in the direction of launch, the positive z-axis vertical, and 
the y-axis completing the right-handed coordinate system. 

2. Launch Centered Earth-Fixed: This is an Earth- fixed coordinate system, having the 
same orientation as the inertial coordinate system (1). This system is advantageous 
in gimbaled inertial platforms in that it is not necessary to remove the Earth rate 
torquing signal from the gyroscopes at launch. 

Figure 2. 1 illustrates two posible methods for defining the missile body axes with 
respect to the Earth and/or inertial reference axes. These coordinate frames will be 
used to define the missile’s position and angular orientation in space. 

Referring to Figure 2. 1, we will denote the Earth-fixed coordinate system by (X e , 
Y e , Z e ). In this right-handed coordinate system, the X e — Y e lie in the horizontal plane, 
and the Z^-axis points down vertically in the direction of gravity. (Note that the posi- 
tion of the missile’s center of gravity at any instant of time is given in this coordinate 
system). The second coordinate system, the body axis system, denoted by (Aft, Yb, Zb), 
is fixed with respect to the missile, and thus moves with the missile. This is the mis- 
sile body coordinate system. The positive A/, -axis coincides with the missile’s center 
line (or longitudinal axis) or forward direction. The positive K/,-axis is to the right of 
the Ai,-axis in the horizontal plane and is designated as the pitch axis. The positive 
Zft-axis is the yaw axis and points down. This coordinate system is similar to the NED 
system. The Euler angles (f, 9, <p) are commonly used to define the missile’s attitude 
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Fig. 2.1. Orientation of the missile axes with respect to the Earth-fixed axes. 



with respect to the Earth-fixed axes. These Euler angles are illustrated in Figure 2.1, 
whereby the order of rotation of the missile axes is yaw, pitch, and roll. This figure 
also illustrates the angular rates of the Euler angles. The transformation CJ ( i from the 
Earth-fixed axes coordinate system to the missile body-axes frame is achieved by a 
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Inertial axes (X t , Y r Z ; ) 




Fig. 2.2. Representation of the inertial coordinate system (inertial. Earth, and body coordinate 
systems). 



yaw, pitch, and roll rotation about the longitudinal, lateral, and normal (i.e., vertical) 
axes, respectively. The resultant transformation matrix C’j’ is [2], [7] 



1 


0 


0 


cos 9 0 — sin 6 


cos ip 


sin xp 0 


0 


cos cp 


sin (p 


0 1 0 


— sin xp 


cos xp 0 


0 


— sin </> cos 0 I 


sind 0 cos 9 


0 


0 1 



cos 9 cos xp cos 9 sin xp — sin 9 

sin 4> sin 9 cos xp — cos (p sin t p sin </> sin 6 sin xp + cos (p cos i fr sin (p cos 9 

cos (p sin 9 cos i p + sin (p sin t p cos (p sin 6? sin t/z — sin (p cos i p cos (p cos 9 



It should be noted here that ambiguities (or singularities) can result from using the 
above transformation (i.e., as 6 , <p, xp -> 90°). Therefore, in order to avoid these 
ambiguities, the ranges of the Euler angles (<p, 6, xp) are limited as follows: 



—u <<p <n or 0 <<p < 2 jt, 

— 7 r <Xp <7T, 

—tt/2<9<tc/2 or 0<\p<2jr. 

The inertial coordinate system described above is shown in Figure 2.2. 



2.2 Rigid-Body Equations of Motion 

In this section we will consider a typical missile and derive the equations of motion 
according to Newton’s laws. In deriving the rigid-body equations of motion, the 
following assumptions will be made: 

1. Rigid Body: A rigid body is an idealized system of particles. Furthermore, it 
will be assumed that the body does not undergo any change in size or shape. 
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Translation of the body results in that every line in the body remains parallel 
to its original position at all times. Consequently, the rigid body can be treated 
as a particle whose mass is that of the body and is concentrated at the center 
of mass. In assuming a rigid body, the aeroelastic effects are not included in the 
equations. With this assumption, the forces acting between individual elements of 
mass are eliminated. Furthermore, it allows the airframe motion to be described 
completely by a translation of the center of gravity and by a rotation about this 
point. In addition, the airframe is assumed to have a plane of symmetry coinciding 
with the vertical plane of reference. The vertical plane of reference is the plane 
defined by the missile Xb- and Z/,- axes as shown in Figure 2.1. The K/,-axis, 
which is perpendicular to this plane of symmetry, is the principal axis, and the 
products of inertia Ixy and Iyz vanish. 

2. Aerodynamic Symmetry in Roll: The aerodynamic forces and moments acting 
on the vehicle are assumed to be invariant with the roll position of the missile 
relative to the free-stream velocity vector. Consequently, this assumption greatly 
simplifies the equations of motion by eliminating the aerodynamic cross-coupling 
terms between the roll motion and the pitch and yaw motions. In addition, a 
different set of aerodynamic characteristics for the pitch and yaw is not required. 

3. Mass: A constant mass will be assumed, that is, dm /dt = 0. 

In addition, the following assumptions are commonly made: 

4. The missile equations of motion are written in the body-axes coordinate frame. 

5. A spherical Earth rotating at a constant angular velocity is assumed. 

6. The vehicle aerodynamics are nonlinear. 

7. The undisturbed atmosphere rotates with the Earth. 

8. The winds are defined with respect to the Earth. 

9. An inverse-square gravitational law is used for the spherical Earth model. 

10. The gradients of the low-frequency winds are small enough to be neglected. 

Furthermore, in the present development, it will be assumed that the missile has 
six degrees of freedom (6-DOF). The six degrees of freedom consist of (1) three 
translations, and (2) three rotations, along and about the missile ( Xb , 17,, Z/,) axes. 
These motions are illustrated in Figure 2.3, the translations being (u, v, w ) and the 
rotations (P, Q, R). In compact form, the traslation and rotation of a rigid body may 
be expressed mathematically by the following equations: 

Translation : F = ma, 

E d 

r — — (rx mV) 
dt 

where is the net torque on the system. 

Aerodynamic forces and moments are assumed to be functions of the Mach* 
number ( M ) and nonlinear with flow incidence angle. Furthermore, the introduction 

*The Mach number is expressed as M = Vm/V s , where Vm is the velocity of the missile 
and V7 is the local velocity of sound, a piecewise linear function of the missile’s altitude. 



(2.13) 

(2.14) 
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Fig. 2.3. Representation of the missile’s six degrees of freedom. 



of surface winds in a trajectory during launch can create flow incidence angles that are 
very large, on the order of 90°. Nonlinear aerodynamic characteristics with respect 
to flow incidence angle must be assumed to simulate the launch motion under the 
effect of wind. Since Mach number varies considerably in a missile trajectory, it is 
necessary to assume that the aerodynamic characteristics vary with Mach number. 

The linear velocity of the missile V can be broken up into components n, v, and w 
along the missile (X/,, K/,, Z/,) body axes, respectively. Mathematically, we can write 
the missile vector velocity, Vm, in terms of the components as 

Vm = ni + i>j T u>k, 

where (i, j, k) are the unit vectors along the respective missile body axes. The mag- 
nitude of the missile velocity is given by 



|VmI — Vm — (u 2 + v 2 + w 2 ) 1 / 2 . 

These components are illustrated in Figure 2.3. 

In a similar manner, the missile’s angular velocity vector w can be broken up into 
the components P, Q, and R about the ( Xb , K/, , Z/,) axes, respectively, as follows: 

a) = Pi + Qj + Rk, 

where P is the roll rate, Q is the pitch rate, and R is the yaw rate. Note that some 
authors use lowercase letters for roll, pitch, and yaw rates instead of uppercase letters. 
Therefore, these linear and rotational velocity components constitute the 6 -DOF of 
the missile. As stated in the beginning of this section, the rigid-body equations of 
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motion are obtained from Newton’s second law, which states that the summation 
of all external forces acting on a body is equal to the time rate of the momentum of 
the body, and the summation of the external moments acting on the body is equal to 
the time rate of change of moment of momentum (angular momentum). Specifically, 
Newton’s laws of motion were formulated for a single particle. Assuming that the 
mass m of the particle is multiplied by its velocity V, then the product 

p = mV (2.15) 

is called the linear momentum. Thus, the linear momentum is a vector quantity having 
the same direction and sense as V. For a system of n particles, the linear momentum 
is the summation of the linear momenta of all particles in the system. Thus [8], 

n 

p = y^(m,V i ) = miVi + m 2 V 2 4 f m„V„, (2.16) 

(=t 

where i denotes the i th particle, and n denotes the number of particles in the system. 
Note that the time rates of change of linear and angular momentum are referred to 
an absolute or inertial reference frame. For many problems of interest in airplane and 
missile dynamics, an axis system fixed to the Earth can be used as an inertial reference 
frame (see Figure 2.1). Mathematically, Newton’s second law can be expressed in 
terms of conservation of both linear and angular momentum by the following vector 
equations [1], [8], [11]: 



E F = 



diniS m) 

dt _ j 



(2.17a) 



E M 



d H" 

dt _ j 



(2.17b) 



where m is the mass, H the angular momentum, and the symbol ]/ indicates the time 
rate of change of the vector with respect to inertial space. Note that (2.17a) is simply 



or 



F = 



dp 

dt ' 



(2.18a) 



fdS\ 

F = m — 1 = ma. 

V dt J 



(2.18b) 



Equations (2.17a) and (2.17b) can be rewritten in scalar form, consisting of three 
force equations and three moment equations as follows: 



F x = 



d(mu) 

dt 



Fy = 



d(mv) 

dt 



F z = 



d(mw) 

dt 



(2.19) 
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where F x , F y , F- and u. v, w are the components of the force and velocity along the 
missile’s Xb, Yb, and Z/, axes, respectively. Normally, these force components are 
composed of contributions due to (1) aerodynamic, (2) propulsive, and (3) gravita- 
tional forces acting on the missile. In a similar manner, the moment equations can be 
expressed as follows [6]: 



dH x dH y dH- 

L = — -, M — y ~, N= -, 

dt dt dt 



( 2 . 20 ) 



where L, M, N are the roll moment, pitch moment, and yaw moment, respectively, 
and H x , FI y , IF are the components of the moment of momentum along the body 
X, Y, and Z axes, respectively. 

At this point, let us summarize the various forces, moments, and axes used in 
developing the missile 6-DOF equations of motion. 

Force: 

F = F x i + F y \ + F,k, 

where F x , F y , F- are the (x, y, z.) components of the force. 

Velocity: 

V = ui + uj + u>k, 

where u, v, w are the velocity components along the (x, y, z) axes, respectively. 
Moment of External Forces: 



^2 AM — Li + Mj + Aik, 

where L is the rolling moment, M is the pitching moment, and N is the yawing 
moment. 

Angular Momentum: 

H = H x \ + Hyj + H z k, 

where H x , H y , H- are the components of the angular momentum along the x,y,z 
axes, respectively. 

Angular Velocity: 

w = co x i + coy j + co z k = Pi + Qj + Rk, 

where P is the roll rate, Q is the pitch rate, and R is the yaw rate, (i = unit vector along 
the x-axis, j = unit vector along the y-axis, and k = unit vector along the "-axis). 

We now wish to develop an expression for the time rate of change of the velocity 
vector with respect to the Earth. Before we do this, we note that in general, a vector 
A can be transformed from a fixed (e.g., inertial) to a rotating coordinate system by 
the relation [6], [7] 



d A\ 


d A 


/ fixed(X'.y'.Z') 


. dt . 



rot.(X.F.Z) 



+ to x A, 



(2.21a) 
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Table 2.1. Axis and Moment Nomenclature 



(a) Axis Definition 



Axis 


Direction 


Name 


Linear 

Velocity 


Angular 

Displacement 


Angular 

Rates 


OX 


Forward 


Roll 


u 


<P 


P 


OY 


Right Wing 


Pitch 


V 


9 


Q 


OZ 


Downward 


Yaw 


w 




R 


(b) Moment Designation 




Moment of Product of 






Axis 


Inertia 


Inertia 


Force 


Moment 


OX 


lx 


*xy 


= 0 


Fx 


L 


OY 


h 


*yx 


= 0 


Fy 


M 


OZ 


h 


hx 


7^0 


F z 


N 



or 



/ d\ ' 




f r/Al 


V dt , 


| 

' inertial 


1 

^3 



■ to x Vm, 



body 



(2.21b) 



where w is the angular velocity of the missile body coordinate system (X, Y, Z) 
relative to the fixed (inertial) system (X', Y Z '), and x denotes the vector cross 
product. Normally, the missile’s linear velocity Vm is expressed in the Earth-fixed 
axis system, so that (2.21a) can be written in the form 



dV m \ 

dt ) E \ dt ) rot. coord. 



■ (0 X Y M, 



( 2 . 22 ) 



where co is the total angular velocity vector of the missile with respect to the Earth. 
In terms of the body axes, we can write the force equation in the form 



F = m 



dV m 
dt 



+ /«(&) x Vm)- 



body 



The first part on the right-hand side of (2.22) can be written as 
( dV m\ (du\, ( dv\ . ddw\ 

v dt j rot . coord . Ur + Ur + Ur 



(2.23) 



(2.24) 



where 



( du/dt ) — forward (or longitudinal) acceleration, 
( dv/dt ) = right wing (or lateral) acceleration, 

( dw/dt ) = downward (or vertical) acceleration, 
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O 



(b) 

Fig. 2.4. General rigid body with angular velocity vector to about its center of mass. 

and the vector cross product as 

" i j k" 

wxVm= PQR = (wQ — vR)i+ (uR — wP)j + (vP — uQ)k. (2.25) 
U V w 

Next, from (2.17a) we can write the sum of the forces as 

£af=E AF,i+^AF v j + ^AF s k. (2.26) 

Equating the components of (2.24), (2.25), and (2.26) yields the missile’s linear 
equations of motion. Thus, for a missile with an X i, — Z[, plane of symmetry (see 
rigid-body assumption #1) we have 



~Y\ A F x =m(u + wQ — vR), 


(2.27a) 


Y, AF y = m{v + uR — wP), 


(2.27b) 


Y, A F z =m(w + vP — uQ ). 


(2.27c) 



From (2.17b) we can obtain in a similar manner the equations of angular motion. 
However, before we develop these equations, an expression for H is needed. To this 
end, consider Figure 2.4. 
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Now let dm be an element of mass of the missile, Y the velocity of the elemental 
mass relative to the inertial frame, and <5 F the resulting force acting on the elemental 
mass. 

First of all, and with reference to Figure 2.4(a), the position vector of any particle 
of the rigid body in a Newtonian frame of reference is the vector sum of the position 
vector of the center of the mass and the position vector of the particle with respect to 
the center of mass. Mathematically, 

fp —^cm "b 7* p/cm - 

where 



r p — the position vector of the particle, 
r cm — position vector of the center of mass of the particle, 
r p / cm — position vector of the particle with respect to the 
center of mass. 



Note that if this equation is differentiated, we obtain 

dr p dr cm dr p /cm 

dt dt dt 

Also, from Figure 2.4(a) we can write the velocity of the point p in the form 



\7 — cm) 

Y p — — hWXt p jcm i 



or 



V p — Y (;m ~h Y p/cm- 

Then, from Newton’s second law we have 



<5F = dm 




(2.28) 



The total external force acting on the missile is found by summing all the elements 
of the missile. Therefore, 

^<$F = F. (2.29) 

The velocity of the differential mass dm is 

V = v„„ + g), (230, 



where \ cm is the velocity of the center of mass (cm ) of the missile, and dr/dt is 
the velocity of the element relative to the center of mass. Substituting (2.30) for the 
velocity into (2.29) results in 




+ 




VSF = F 



dm. 



(2.31) 
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Assuming that the mass of the missile is constant, (2.31) can be written in the form 



or 




(2.32a) 



(2.32b) 



Since r is measured from the center of the mass, the summation is equal to 0. 

Thus, the force equation becomes simply 



F = m 




(2.33) 



which relates the external force on the missile to the motion of the vehicle’s center 
of mass. Similarly, we can develop the moment equation referred to a moving center 
of mass. For the differential element of mass, dm, the moment equation can then be 
written as 



8M = d 





(r x Y)dm. 



(2.34) 



The velocity of the mass element can be expressed in terms of the velocity of the 
center of mass and the relative velocity of the mass element to the center of mass. 
Therefore, 



V P =V C 



f dr pj C m 

V dt 



= Y cm + w X r, 



(2.35) 



where co is the angular velocity vector of the vehicle and r is the position of the mass 
element measured from the center of mass (see Figure (2.4a)). In relation to (2.35) 
and Figure 2.4(a), we can write the equation 



/Jr\ 


' d r" 


V dt / inertial 


_ dt _ 



+ co x r. 



Jrel. to coord. 



The reader will note that this is the well-known Coriolis equation, which is important 
in dynamics where body axes are used. Furthermore, it will be noted that the term 
<0 x r occurs in addition to the vector change relative to the coordinate system, so 
that the total derivative relative to the inertial axes is expressed by this equation. The 
rigid-body assumption implies that dr p / cm /dt — 0. Therefore, we can write the linear 
velocity of the point p in the simple form 



Y p — co x v p f cm 

In general, the moment about an arbitrary point O of the momentum p = mV (2.15) 
of a particle is 



H = r x mV = mr x (w x r). 
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Referring to Figure (2.4(b)), it will be observed that this vector is perpendicular to 
both r and V. Furthermore, this vector lies along MP and is directed toward M. The 
moment of momentum (or angular momentum) of the entire body about O is therefore 

H = r x m V = mr x (to x r) = ^2 m [w(r • r) — r(r • w)] 

(note that this result was obtained using the formula a x (b x c) = (a • c)b - 
(a • b)c). This equation can also be written as 



H = mr ~ 



2 ) to — mr ( r ■ w). 

From Figure (2.4(a)), the total moment of momentum can be written as 

H = <$H = ^(r x V cm )dm + ?> x (co x r )]dm. (2.36) 



Note that the velocity V cm is constant with respect to the summation and can be taken 
outside the summation sign. Thus [1], [3] 

h =£ r dm x Vcm + ^[rx(uxr)](im. (2.37a) 

As stated above, the first term on the right-hand side of (2.37a) is zero. Therefore, we 
have simply 



and 



m=£[r x (w x r )]elm 




(w x r )dm. 



Performing the vector operations in (2.37c) and noting that 



(2.37b) 



(2.37c) 



we have 



Finally, 



w = co x i + (o y j + w-k = Pi + Qj + Rk, 
r = xi + yj + zk. 



wxr = 



i j k 

P Q R 



lx y z J 

= (zQ- yR ) i + (xR -zP ) j + (yP - xQ) k. 



(2.38a) 



r x (to x r) 



i j k 

x y z 

_(zQ-yR ) (. xR-zP ) (yP-xQ)_ 

i[(v 2 + z 2 )P — xyQ — xzR] + j[(r + x 2 ) Q - yzR - xyP] 

+ k[(x 2 + y 2 )R-xzP-yzQ]. (2.38b) 
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Substituting (2.38b) into (2.37c), we have 






i[(v 2 +z 2 )P -xyQ -xzR]dm + J j [(z 2 + x 2 )Q — yzR - xyP]dm 



7 



k[(x 2 + y 2 )/? — xzP — yzQ]dm, 



(2.38c) 



where the f ( y 2 + z 2 ) is defined as the moment of inertia, I x , and f xydm is defined 
as the product of inertia, I xy . The remaining integrals in (2.38c) are similarly defined. 
By proper positioning of the body axis system, one can make the products of inertia 
l xy — I yz equal to 0. This will be true if we can assume that the x-y plane is a plane of 
symmetry of the missile. Consequently, (2.38c) can be rewritten in component form 
as follows: 



H x = 
Hy = 
H z = 



7 



7 



(y 2 + z 2 )dm — R I xzdm — PI X — RI X 



Q J (x 2 + z 2 )dm = QI y , 
R J(x 2 + y 2 )dm - P J 



xzdm — RL — P I x 



(2.39a) 

(2.39b) 

(2.39c) 



From (2.17b), we note that the time rate of H is required. Now, since H can change 
in magnitude and direction, (2.17b) can be written as [ 1] 



AM=1„ 




+ w x H. 



Next, the components of 1 n(dH/dt) assume the form 



(2.40) 



dH x _(dP\ ( dR \ 

~dt~ ~ \dT ) Ix ~ \~dt~ J 

dt \ dt ) v ’ 
dH z _/dR\ ( dP \ 
dt \ dt ) \dt ) 



Ixz- 



(2.41a) 

(2.41b) 

(2.41c) 



Since initially we assumed a rigid body with constant mass, the time rates of change 
of the moments and products of inertia are zero. The vector cross product in (2.40) is 



(iixH = 



i j k 

P Q R 

H X Hy Hy 



= (QH Z - 



RH y )i+(RH x - PH Z ) j + (PH y - QH X ) k. (2.42) 
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Similar to (2.26), we can write an equation for the summation of all moments in 
the form 



AM = i^ AL+jJ] AM + k^ AN. (2.43) 

Equating the components of (2.41), (2.42), and (2.43) and substituting for H x , H y , 
and H z from (2.39) yields the angular momentum equations. Thus [1], [5], 



PI * - Khz + QR(h - Iy) - PQIxz • 


(2.44a) 


J2 A M = Qh + PR(1 * - 7 z> + {P 2 - Rl )hz, 


(2.44b) 


N = SIz- Phz + PQ(h - h) + QRI XZ , 


(2.44c) 


J2 AL = p h + (h - h)QR~ (R + P Q)hz, 


(2.44d) 


J2 A M = Qh + (l * ~ h)PR + ( P 2 ~ R 2 )h z , 


(2.44e) 


J2 AN = Pl z + (h - h)PQ- (P - QR)hz - 


(2.44f) 



where dP /dt is the roll acceleration, dQ/dt is the pitch acceleration, and dR/dt 
is the yaw acceleration. The set of equations (2.27a)-(2.27c) and (2.44d)-(2.44f) or 
(2.44a)-(2.44c) represents the complete 6- DO I- missile equations of motion. Specifi- 
cally, equations (2.27) describe the translation, and equations (2.44) describe the 
rotation of a body. The set of equations (2.27) and (2.44) are six simultaneous non- 
linear equations of motion, with six variables u, v, w . P . Q, and R. which com- 
pletely describe the behavior of a rigid body. Moreover, these equations can be 
solved with a digital computer using numerical integration techniques. An analytical 
solution of sufficient accuracy can be obtained by linearizing these equations. These 
equations are also known as Euler’s equations. Note that I x , I y , I xz are constant 
for a given rigid body because of our choice of coordinate axes. Due to the usual 
symmetry of the aircraft (or missile) about the x-y plane, the products of inertia 
that involve y are usually omitted, and the moment equations may be rewritten as 
follows (note that for cruciform missiles with rotational symmetry, I y — I z and 
hz = 0): 

A L = PI x + QR(I z -I y ), (2.45a) 

A M= QI y + (I x — I Z )PR, (2.45b) 

A N = RI Z + (I y - / x)PQ. (2.45c) 

It should be noted that the L and N equations indicate that a rolling or yawing moment 
excites angular velocities about all three axes. Therefore, except for certain cases. 
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Fig. 2.5. Rotational dynamics of a rigid body. 



these equations cannot be decoupled. Solving (2.45a)-(2.45c) for d P /dt, d(Q/dt , 
and dR/dt, we obtain the rotation accelerations as follows: 



— = QR[(Iy-I z )/I x ] + (L/I x ), 
dt 


(2.46a) 


^ = PR[(I Z - I X )/Iy] + (M/Iy), 


(2.46b) 


d ^ = PQWx-Iy)/I Z ] + (N/I z ). 


(2.46c) 


The relationship of the three coordinate systems discussed in 


Section 2.1 can be 


described in terms of the body dynamics. Figure 2.5 illustrates the manner in which 
these three methods are integrated into computational sequence of representing the 


vehicle dynamics. 




The equations for the angular velocities (< i\jr / dt , dtj) / dt , dG /dt) in terms of the 
Euler angles {\jr, cf>, 0) and the rates ( P , Q. R) can be written from Figure 2.1 as 


follows [1]: 




d\[/ 

— = (Q sin cp + R cos (p) / cos 9 , 
dt 


(2.47a) 


d l =P+ (i 'iW 

dt \ dt ) 


(2.47b) 


dd 

— = Q cos (p — R sin (p , 
dt 


(2.47c) 



where P is the roll rate, Q is the pitch rate, and R is the yaw rate. The values of 
(t/r, <j > , 6) can be obtained by integrating (2.47a)-(2.47c). Thus, 

(?)*• 



(2.48a) 
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</> = 0O+ [ 

Jo 


{ dtp ^ 
V dfy 


) dt, 


(2.48b) 


e = o 0 + f 

Jo 


(dd\ 

V fifty 


) dt. 


(2.48c) 



From the transformation matrix C b e of Section 2.1, the components of the missile 
velocity dX e /dt, dY e /dt, dZ e /dt in the Earth-fixed coordinate system ( X e , Y e , Z e ) 
in terms of (u, v, w ) and (t/r, cp, 0) are given as follows: 



dX e 

dt 

dY e 

dt 

dZ e 

dt 

or in matrix form. 



(cos 0 cos \[r )u + (cos t jr sin 0 sin 0 — sin xfr cos 0)t> 
+ (cos cos 0 sin 6 + sin sin 0) w , 

(cos 6 sin t/f )u + (sin xp sin 0 sin 9 + cos i p cos (p)v 
+ (sin i p cos (p sin 0 — cos t/r sin0)ui, 

— (sin0)t< + (sirup cos6)v + (cos 0 cos <p)w, 





— c h 


u 

V 






w 



From (2.49) we can obtain the equations for (X e , Y e , Z e ) in the form 




and the altitude is 



(2.49) 



(2.50a) 

(2.50b) 

(2.50c) 



h — Z e . 



(2.50d) 



In the foregoing discussion, only the missile velocities relative to the ground or inertial 
velocities have been mentioned. If wind is being considered, the missile velocities 
relative to the wind must be computed, since these velocities are needed in computing 
the aerodynamic forces and moments (see Chapter 3). 

It should be noted here that stability and control for fixed-wing aircraft are assessed 
through six rigid-body degrees of freedom models. Rotorcraft models provide three 
more degrees of freedom for the main flapping plus a rotational degree of freedom. 
Additional degrees of freedom for structural modes and other dynamic components 
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Fig. 2.6. Coordinate system 



(e.g., transmissions) can be added as necessary. The aircraft models are based on 
aerodynamic coefficient representations of the major aircraft components, including 
the wing, fuselage, vertical tail, and horizontal tail. Mass distribution is represented 
by the center-of-gravity location and mass moments of inertia for the aircraft. A 
stability analysis is performed by trimming the forces and moments on the air- 
craft model for each flight condition. Force and moment derivatives are obtained 
through perturbations from trim in the state and control variables. These derivatives 
are used to represent the rigid-body motion of the aircraft as a set of linear first-order 
differential equations. The matrix representation of the aircraft motion is then used 
in the linear analysis package MATLAB* to assess stability and to investigate feed- 
back control design. Aircraft dynamics and control system conceptual designs are 
typically analyzed with respect to dynamic performance, stability, and pilot/vehicle 
interface. 

Example 1. In this example, we will consider an aircraft whose equations of 
motion can be represented as a point mass, based on five variables (i.e., 5 -DOF). 
The coordinate system for this example is illustrated in Figure 2.6. 

Furthermore, the variables are defined as follows: 

Let a ' , y,z = position variables, 

Y = velocity vector, 
a = angle-of-attack (AOA), 
ft = velocity heading angle, 
y = velocity elevation angle, 

p = orientation angle of the aircraft body axes relative 
to the velocity vector. 

From the definition of the above variables, the orientation of the velocity vector Y is 
through the angles ft and y, while the orientation of the aircraft body axes relative 

* MATLAB is a commercially available software package for use on a personal computer. 
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to the velocity vector is through the angle p and the AOA a in the pitch plane. The 
yaw of the aircraft about the velocity vector, V, is assumed to be zero (i.e., sufficient 
control power exists that all maneuvers are coordinated). For this 5 -DOF (x, y, z, 
a, /a) point mass model, the equations of motion are as follows: 

x — V cos y cos /3, 
y — V cos y sin/3, 
z — V sin y, 

1 

V — —[T cos a — D] — g sin y , 
m 

1 

/3 = [T sina + L](sin/x/ cos y), 

m V 

1 

y = [T sina + L] cos p — (g/F)cosy, 

m V 

m = m(M, z, n), 

h = n(a, Vi as ). 
a = a (a, V rAS ), 

A = A(«. Vias ), 

where 

M — Mach number, 
g — acceleration of gravity, 
m = mass, 
n = throttle setting, 

Vj as — indicated airspeed, 

T — Thrusts T (M, z, n), 

D — drag= \pV 2 SCo (see Section 3.1), 

L — lift = 5 pV 2 SCl (see Section 3.1), 
p — atmospheric density = p (z) 

(i.e., a function of altitude), 

S — aerodynamic reference area, 

Cd — coefficient of drag. 

Cl = coefficient of lift. 

Several approximations can be made in the above model. These are: 

1. The dp/dt equation of motion becomes undefined for vertical (i.e., y = ± 90°) 
flight. 

2. Thedn/dt, da/dt . dp jdt equations are at best first-order approximations to actual 
aircraft control response. 
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From the above results and discussion, a 6 -DOF model can be implemented that 
approximates an actual 6 -DOF control response with a standard transfer function/filter 
whose input constants can be selected by the designer to more accurately match actual 
aircraft/missile control response. The roll, pitch, and yaw transfer functions are then 
as follows: 



Roll: 

Pitch: 

Yaw: 



^stab lb ) 1 

Atab cmd (^) ~~ TS + 1 
n z (s ) co 2 (rs + l) 

' h - cmd ^) s 2 + 2$ojs+w 2 
n y (s ) _ co 2 

n y cmd(^) s 2 + 2i;ws + co 2 



where r is the time constant, .v is the Laplace operator, and o> is the frequency. 

Under 6 -DOF modeling, the dyt/dt, dy /dt , and dft/dt kinematic relationships 



are 



where 



dyt 

dt 

dy 

dt 

dp 

dt 



P + tan y(Q sin/r + R cos yt), 

Q cos /a — R sin /z, 

sec y ( Q sin yt + R cos yt ) , 



P — body axes roll rate, 

Q = body axes pitch rate, 

R — body axes yaw rate. 

Next, in order to eliminate the dfi/dt equations anomaly at y = ± 90°, the quaternion 
system of coordinates will be used; the kinematic rate equations are [7] 



cle i 
i it 
de 2 

dt 
de 3 
i it 
de 4 
dt 



(~e\P — e-$Q — eiR)/2, 
(~e 3 P -e 4 Q-eiR)/2, 
{-eiP + e\Q - e$R)/2, 
{-e\P - e 2<2 + e^R) /2, 
and the Earth-to-body direction cosine matrix is, as before. 



C 



b 

e 



Cll C\2 Ci 3 
C21 C22 C23 
C31 C32 C33 
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where 

Cn = e\-e\-e\ + e\, 

C 12 = 2(e3eq + e\e2), 

C 13 = 2(e2<?4 — et<?3), 

C21 = 2(^364 — eie2), 

n 2 2,22 

C22 — e l e 2 + e 3 ~ e 4’ 

C23 = 2(^263 + e\e 4 ), 

C31 = 2(ei?3 + 62^4), 

C32 = 2(^2e3 — ei<?4), 

^ 2,222 

C33 = e } + e 2 — e 3 — e 4 , 

and 

P = tan- 1 (Ci 2 /Ci 1 ), 

K = - sin -1 (C13), 

H = tan -1 (C23/ C33). 

Finally, we note that the same 6 -DOF equations of motion can be used to model both 
aircraft and missiles. 

Example 2. Based on the discussion thus far, let us now consider in this exam- 
ple a 6-DOF aerodynamic model. Furthermore, let us assume an NED coordinate 
system, in which all units are metric. This model is designed for a generic aircraft. 
A quaternion fast-processing technique will be employed to simulate aircraft navi- 
gation. This technique avoids not only time-consuming trigonometric computations 
in the fast-rate direction cosine updating, but also singularities in aircraft attitude 
determination [7], 

6 -DOF Initialization 

Before processing begins, these initialization functions must be performed. Compute 
the initialized Earth reference velocity: 

U e — V * cos (6) * cos (i/O, 

V e — V * cos (<p) * sin(i /r), 

W e = —V * sin(0), 

where 

U e = Earth X- velocity, 

V e = Earth Y -velocity, 

W e = Earth Z-velocity, 

6 = pitch angle, 

(f> — roll angle, 
ijr — yaw angle, 

V — airspeed. 
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Compute Initialized Quaternions 

A = [sin(i//-/2) * sin(0/2) * cos(</>/2)] — [cos(i/r/2) * cos(0/2) * sin(</>/2)], 

B = — 1 * [cos(i/r/2) * sin(0/2) * cos(<£/2)] — [sin(i/r/2) * cos(0/2) * sin(0 /2)], 

C = —1 * [sin ( i/z/2) * cos(0/2) * cos(0/2)] + [cos(t/r/2) * sin(0/2) * sin(</>/2)], 

Z) = — 1 * [cos(t/f/2) * cos(0/2) * cos(0/2)] — [sin(t/r/2) * sin(0/2) * sin(</>/2)], 

where A, B, C, D — quaternion parameters of the direction cosine matrix. 

Now compute the initialized direction cosine matrix: 

Cm( 1, 1) = A 2 — B 2 - C 2 + D 2 , 

Ctn( 1 , 2) = 2 * (A * 5 — C * Z)), 

Cm(l, 3) = 2 * (A * C + Z? * Z)), 

Cm(2, 1) = 2 * (A * Z? + C * Z)), 

Cm( 2,2) = -1*A 2 + B 2 -C 2 + D 2 , 

Cot (2, 3) = 2*(Z?*C — A* Z)), 

Cm (3, 1) = 2 *{A*C-B*D), 

Cm( 3, 2) = 2*(fi*C + A* Z)), 

Cot ( 3, 3) = - 1 * A 2 - B 2 + C 2 + D 2 , 
where Cot direction cosine matrix. 

Compute the initial body velocity 

U b = Cot( 1, 1) * Z/ e + Cot ( 2, 1) * V e + Cot ( 3, 1) * W e , 

Vb = Cm(l, 2) * U e + Cot (2, 2) * V, + Cot ( 3, 2) * W e> 

= Cm(l, 3) * U e + Cot ( 2, 3) * V c + Cot (3, 3) * 

where 

Ub — body A- velocity, 

Vb — body K-velocity, 

Wb — body Z- velocity. 

6 -DOF Processing 

The following computations are performed at every simulation cycle. 

Compute the dynamic pressure 

q — \pV 2 , 

where 

q — dynamic pressure, 
p — pressure in standard atmosphere, 

V — airspeed. 
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Compute the wing lift 



where 



L = Cl*ci* S, 



L = lift, 

Cl = coefficient of lift, 

S = surface area of the wing. 

Compute the wing drag 

D = Cr>*q*S, 

where Co — coefficient of drag. 

Compute the lift acceleration 

L a = L l w , 

where w = weight of the airplane. 

Compute the drag acceleration 

D a = D/w. 

Compute the thrust acceleration 

T a — T/m , 

where 

m — mass of the airplane, 

T = thrust. 



Compute the body accelerations 

Xba — T a * L a * sin(a) — D a * cos(a) + Cm( 3, 1) * g + R * Vb — Q * Wb, 
Yba — Cm( 3, 2 ) * g — R * Ub + P * Wb , 

Zba — — 1 * L a * cos(a) — D a * sin(a) + Cm( 3, 3) * g + Q * Xb — P * Vj, 
where 

Xba — 2f-axis body acceleration, 

Yb a = Y -axis body acceleration, 

Zba — Z - axis body acceleration, 

P — roll rate, 

Q — pitch rate, 

R — yaw rate, 

g — acceleration due to gravity, 
a — angle of attack. 
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Compute the Earth accelerations 

X ea = Cm{ 1, 1) * Xba + Cm( 1, 2) * Y ba + Cm( 1, 3) * Z^, 

= Cw (2, 1) * + Cot (2, 2) * Y ba + Cm{ 2, 3) * Z ba , 

Z ea = Cot (3, 1) * + Cm (3, 2) * + Cot (3, 3) * Z ba , 

where 

X ra = X-axis Earth acceleration, 

Y ea = Y -axis Earth acceleration, 

Z ea = Z-axis Earth acceleration. 

Compute the angular deltas 

A 0 — Q*ti, 

A <p — P *tj, 

At// = R*t(, 

where 

A0 = pitch delta, 

A cp = roll delta, 

At j/ = yaw delta, 

ti = simulation cycle time. 

Compute Cn and Sn 

Cn = 1.0 — (A0 2 + A</> 2 + Ai/r 2 )/8 + (A0 4 + A</> 4 + At/f 4 )/384, 
Sn = 0.5 - (A 6 2 + A(f> 2 + Ai// 2 )/48, 

where 

Cn = nth-order Maclaurin Series of cos(A0/2), 

Sn = nth-order Maclaurin Series of sin(A0/2), 

A 6 = total body angle increment in f,- . 

Compute the Quaternions 

A' — A* Cn + B * Sn * Ai/r + C * — 1 * Sn * A 0 + D * Sn * A </>, 
B' — A * — 1 * Sn * A\[r + B * Cn + C * Sn * A </> + D * Sn* Ad, 
C' — A * Sn * A 6 + B * — 1 * Sn * A <p + C * Cn + D * Sn * At//, 
D' — A * — 1 * Sn * A </> B * — 1 * Sn * A 9 + C * — 1 * Sn * Ai j/, 
A = A', 

B = B' , 

C = C', 

D = D' . 
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Renormalize the Quaternions 

Normalizer = 0.5 * (3.0 — A 2 — B 2 — C 2 — D 2 ), 

A = A* Normalizer, 

B = B * Normalizer, 

C = C * Normalizer, 

D — D * Normalizer. 

Compute the direction cosine matrix 

Cm(l, 1) = A 2 - B 2 - C 2 + D 2 , 

Ctn( 1 , 2) = 2 * (A * B — C * D ), 

Cm{ 1,3) = 2*(A*C + 5* D ), 

Cwi(2 , 1) = 2 * (A * B + C * D ), 

Cm (2, 2) = — 1 * A 2 + B 2 — C 2 + O 2 , 

Cm (2, 3) = 2 * (B * C — A* D), 

Cm( 3, 1) = 2 *(A*C-B*£>), 

Cm( 3,2) = 2*(B*C + A*£)), 

Cm (3, 3) = — 1 * A 2 — B 2 + C 2 + D 2 . 

Compute the Earth velocities 

U e — U c + * tj , 

Ve — Ve + Y ea * tj , 

W e — W e + Z ea * ti . 

Compute the body velocities 

U b = Cm{ 1, 1) * U e + Cm (2, 1) * y e + Cm(3, 1) * lV e , 

= Cm(l, 2) * C e + Cm (2, 2) * y e + Cm(3, 2) * 1V C , 

= Cm(l, 3) * £/ e + Cm (2, 3) * y e + Cm(3, 3) * W,,. 

Compute the airspeed 

V = (U 2 +V 2 +W 2 ) l/2 . 

Compute the Earth referenced position 

U — U e * / i , 

v e = y e * t,- , 

= W e *r,'. 
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Compute the angle of attack 

a — Atan2(Wb/Ub). 

Compute the attitudes 

6 = A sin(— 1 * Cm(3, 1)), 
f = Atan2(Cm(3,2)/Cm(3, 3)), 
t jf = A tan2(C/«(2, 1)/Cw(l, 1)). 

Compute the sideslip angle 

P = A tan2(V b /U b ), 

where f is the sideslip angle. 

Compute the flightpath angle 

Y = A tan2((— 1 * W e )/{U 2 e + V, 2 ) 1 / 2 ), 

where y is the flightpath angle. 

Earlier in this section, the equations of motion for a missile were discussed assum- 
ing the missile to be a rigid body. However, all materials exhibit deformation under 
the action of forces: elasticity when a given force produces a definite deformation, 
which vanishes if the force is removed; plasticity if the removal of the force leaves 
permanent deformation; flow if the deformation continually increases without limit 
under the action of forces, however small. 

A "fluid” is material that flows. Actual fluids fall into two categories, namely, 
gases and liquids. A “gas” will ultimately fill any closed space to which it has access 
and is therefore classified as a (highly) compressive fluid. A “liquid” at constant 
temperature and pressure has a definite volume and when placed in an open vessel 
will take under the action of gravity the form of the lower part of the vessel and 
will be bounded above by a horizontal free surface. All known liquids are to some 
extent compressible. For most purposes it is, however, sufficient to regard liquids as 
incompressible fluids. It should be pointed out that for speeds that are not comparable 
with that of sound, the effect of compressibility on atmospheric air can be neglected, 
and in many experiments that are carried out in wind tunnels the air is considered to 
be a liquid, in the above sense, which may conveniently be called incompressible air. 

All liquids (and gases) in common with solids exhibit viscosity arising from 
internal friction in the substance. For those readers interested in pursuing more thor- 
oughly the area of incompressible air and/or fluid flow, the Navier-Stokes equation 
is a good start. The Euler and Navier-Stokes equations describe the motion of a fluid 
in R" (n — 2 or 3). These equations are to be solved for an unknown velocity vector 
u(jc, t ) = (u,- (jc, t))i<i< n e R" and pressure p(x, t) e R. defined for position x e R" 
and time t > 0. We restrict attention here to incompressible fluids filling all of R". 
The Navier-Stokes equations are then given by 

n 

(d/dt)uj + Uj(dUj/dXj ) = vAui — ( dp/dxi ) + f\ (x , t ) (reR'',/> 0), 
j = 1 



(2.51) 
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n 

divu = J2(d»i/dxi) = 0 (jc e R", t > 0), (2.52) 

i=i 



with initial conditions 



u(jc, 0) = u°(jc), jc e R" , 



(2.53) 



where 

A = ^(9 2 /9*, 2 ) 

i=i 

is the Laplacian in the space variables. 

Here, h°(jc) is a given C°° divergence-free vector field on R”, f ) (x. t) are the 
components of a given externally applied force (e.g., gravity), and v is a positive 
coefficient (the viscosity). Equation (2.51) is just Newton’s law / = ma for a fluid 
element subject to the external force/ = ( f , (;c . f))i<i<« and to the force arising from 
pressure and friction. Equation (2.52) just says that the fluid is incompressible. For 
physically reasonable solutions, we want to make sure that «( jc, t) does not grow large 
as |jc| —> oo. Moreover, we accept a solution of equations (2.5 1)— (2.53) as physically 
reasonable only if it satisfies 



and 



p,ue C°°(R" x [0, oo)) 




|h(x, t)\ 2 dx < C 



Vf > 0 (bounded energy). 



2.3 D’Alembert’s Principle 

In Section 2.2 we discussed the rigid-body equations of motion. Specifically, we 
discussed Newton’s second law as given by (2.17b) and (2.18b). In the fundamental 
equation (2.18b), F — m a, the quantity m(— a) is called the reversed effective force 
or inertia force. D’Alembert’s principle is based on Newton’s second and third laws 
of motion and states that ‘the inertia force is in equilibrium with the external applied 
force,’ or 



F + m(— a) = 0. (2.54) 

This principle has the effect of reducing a dynamical problem to a problem in statics 
and may thus make it easier to solve. Based on the principle of virtual work, * which 
was established for the case of static equilibrium, we can proceed as follows: Let p be 
the momentum of a particle in the system, and separate the forces acting on it into an 

* Consider a particle acted upon by several forces. If the particle is in equilibrium, the resul- 
tant R of the forces must vanish, and the work done by the forces is a virtual displacement 
c$r is zero. Thus, R • <5r = 0. 
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applied force F and a constraint force f. Then the equation of motion of the particle 
can be written as [9] 



F + f- 




= 0 . 



The quantity (dp/dt) is usually referred to as the reverse effective force discussed 
above. Note that the virtual work of the constraint force is zero, since f and <5r are 
mutually perpendicular. The virtual work of the forces acting on the particle is 




• Sn — 0 (i = l,2, 



and for a system of N particles. 



N 



E 



F, 




• Sri — 0 (i= 1,2, 



Another way of writing this equation is 



N 



E 





. Srt = 0 0 = 1,2,..., N), 



where r, is the position vector of the particle. The term —mj(d 2 rj /dt 2 ) has the 
dimensions of force and is known as the inertia force acting on the ith particle (see 
also discussion above). This is the Lagrangian form of d’Alembert’s principle and is 
one of the most important equations of classical dynamics. 



2.4 Lagrange’s Equations for Rotating Coordinate Systems 

The missiles considered thus far were assumed to obey the laws of rigid bodies. 
However, in analyzing the dynamics of flexible missiles, such as intermediate-range 
ballistic missiles or intercontinental ballistic missiles, it is convenient to use a set of 
coordinates moving with the missile. In this case, the missile can be considered as 
a system of particles whose position relative to the moving axes can be defined by 
generalized coordinates qt . Specifically, we will consider the motion of a holonomic* 
system with n degrees of freedom. Let {q \ , qi , . . . , q n ) be the coordinates that specify 
the configuration of the system at time t. Furthermore, we will consider a mechanical 
system of n particles whose coordinates are (xi, yi, zi, X 2 , yi, Z 2 , ■ ■ ■ , x n , y„, z n )■ 
The motion of the system is known when the value of every coordinate is known as 
a function of time. Suppose now that the system moves from a certain configuration 
given by (X| , . . . , z' n ) at time t\ to another configuration given by ( x ", y”, . . . , z") at 

*A dynamical system for which a displacement represented by arbitrary infinitesimal 
changes in the coordinates is, in general, a possible displacement is said to be holonomic. 
When this condition is not satisfied, the system is said to be nonholonomic. 
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time t 2 . During all of the motion between these two configurations, the Newtonian 
equations of motion will be followed, and the acceleration of each particle will be 
given by the total force acting on it. Moreover, this motion can be described by 
expressing each coordinate as a function of time. 

There are then 3 n dependent variables depending on the one independent variable 
t. These functions can be written in the form 

X\=x\(t), yi — y\ (f), . . . , Zn=Zn(t). (2.55) 

In deriving the equations of motion, it is common practice to start the derivation 
using the concepts of kinetic and potential energies of the system using Lagrange’s 
equation. As will be noted, however, these equations differ from the usual Lagrange 
equations for fixed coordinates. Now consider some other way in which the system 
might have moved from the initial configuration to the final configuration in the 
same amount of time, t 2 — t\ . This new motion is to be one that satisfies the geometric 
conditions, or the constraints of the problem. If this new motion is just slightly different 
from the original motion, the coordinates, as functions of time, can be written as 
follows: 

xi ( t ) + 8xi ( t ), yi ( t ) + <5yi (f), . . . , z n (?) + $z n (?)• 

The variation of a coordinate x is a function of time and is the difference between 
the x coordinate of the comparison path and that of the true path. It is also 
assumed that the true path is a continuous function with continuous first derivatives 
satisfying Newton’s equations. The same ideas apply to the comparison path. 
Therefore, 



8x\ (?i) = <5*t (ti) = <5yi (t\ ) = 8y\ (t 2 ) = ■■■= 8z n fe) = 0. (2.56) 



The true path was originally defined in terms of the Newtonian equations of motion. 
For the true path there are 3 n equations of the form 



mi 




(2.57) 



where m, typifies the mass of one of the particles of the system. The quantity A, may be 
a function of the coordinates, of the time explicitly, or of both. It may be considered, 
however, as a function of time only, since the dependence on the coordinates is a 
dependence upon the positions of the particles, and these are uniquely determined 
by the time along any path that may be considered. In general, the coordinates of 
the individual particles (referenced to some fixed set of rectangular coordinates) are 
known functions of the coordinates (q\ , c/ 2 ,..., q n ) of the system, and possibly of t 
also. Let this dependence be expressed by the equations [11] 



Xi = fi(qi,q 2 , • • • , qn , 1 ), 
yi = giiquqi — , q„,t), 

Zi = hj(q \ , q 2 , q„,t). 



(2.58) 
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Furthermore, let (Xj, Yj, Zj ) be the components of the total force (external) acting on 
the particle m,-. Then, the equations of motion of this particle are 



d z Xj 

~di 2 



= X U 



d 2 yj 

, i it 2 



= Yi 



d~z j 
~dt 2 



= Z h 



(2.59) 



If each component of the force Xj is now multiplied by the variation of path in the 
direction of the force and all the resulting equations are added together, the result is 
[4], [11] 



SU = J2 (X > Sx i + YiSyi + ZjSzi ) 

i 

= E 



' d 2 Xj d 2 Vj 

mi i — —8xj + 



dt 



dt 



d 2 Zi 
dt 2 



2 to + 






— (xjSxj + yjSyi + ij&Zi) - XjSxj - y,<5y; - ij&Zi 
dt 



(2.60) 



where the symbol ^ denotes summation over all the particles of the system; this can 
be either an integration (if the particles are united into rigid bodies) or a summation 
over a discrete aggregate of particles. The quantity 8U is defined by the first equality 
in (2.60). It is the work done by the forces of the system during the infinitesimal 
displacement (8xj , . . . , 8z n ) and is a function of the time and the independent coor- 
dinates of the system. If the forces do not depend explicitly on the time, SU can be 
expressed as a function of the coordinates only. The last part of (2.60) represents the 
variation of the kinetic energy ST. Hence the equation can be written as [3] 



8T- 






m j — (xj 8xj + yj 8y t + Zj Szj). 
dt 



(2.61) 



It should be noted that in the above expressions t is the independent variable. Now, if 
both sides of (2.61) are integrated with respect to this independent variable between 
the limits t\ and t 2 , the result is 



f 

Jh 



(8T + 8U)dt = 



r'2 rt: 

= 8 Tdt + I 

J tt Jti 



8Udt — 0. 



(2.62) 



In this equation, we note that the right-hand side is zero because all of the variations 
are zero at both limits. Therefore, (2.62) is a property of the path that satisfies the 
equations of motion, and this property furnishes a way of defining the true path of the 
system. In the special case in which the forces are conservative, that is, when they 
can be derived from the potential energy, 8U is the negative of the variation of the 
potential energy. Consequently, we have 



ft 2 ft 2 

8 / (T — U)dt — 8 / Ldt — 0, 

Jt\ Jt\ 



(2.63) 



where T is the kinetic energy and U is the potential energy. 
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Commonly, the quantity (T — JJ) is denoted by L and is called the Lagrangian 
function or the kinetic potential of the system, or L — T — U. The function L, 
defined as the excess of kinetic energy over potential energy, is the most fundamen- 
tal quantity in the mathematical analysis of mechanical problems. The Lagrangian 
function can be expressed in any convenient coordinate system, and the variation 
principle will still apply. Thus, if we introduce a new function L of the variables 
(qi, qi, ■ ■ ■ , qn, qi, qi, • • • , qn, t), defined by the equation 

L=T — U, 



then Lagrange’s equations can be written in the form [ 1 1] 



d /9L\ dL 
dt \ dqi ) dqt 



i = 1, 2, . . . , n, 
t\ <t<t 2 . 



(2.64) 



Hamilton’s principle for the motion of a mechanical system states that 

r> 2 

8 L(qi,q 2 ,...,q n ,qi,q 2 ,---,qn,t)dt. (2.65) 

Jti 

In (2.65) the q’s represent the coordinates necessary to specify the configuration of the 
system. Note that the time appears explicitly in the Lagrangian function only in case 
the forces are explicit functions of time, or the coordinates used are in motion. In the 
simple conservative cases the Lagrangian function depends upon the coordinates and 
their first derivatives only. If, as has been assumed, the coordinates are all independent, 
then the path can be described by the set of differential equations (2.64). 

The Euler-Lagrange equations for Hamilton’s principle (2.64) are usually called 
simply Lagrange’s equations. They contain nothing more than was contained in the 
Newtonian equations, but they have the decided advantage that the coordinates may 
be of any kind whatever. It is necessary only to write the potential and kinetic energies 
in the desired coordinates to obtain the equations of motion by simple differentiation. 
This is usually much simpler than transforming the differential equations themselves. 
Finally, we note here that Lagrange’s equations and Newton’s equations are entirely 
equivalent. 

Now, if the ( x-y ) plane is rotated by an angle 6 , the coordinate axes in plane 
motion will have three degrees of freedom, namely, xo , yo, and 6, which can be 
varied independently. In this case, the Lagrange equations can be written in the form 
[4], [9] 



d dT ■ dT 
dt dxo dyo 

d dT . dT 

— + 0 — 

dt dyo dxo 




d dT . dT . dT ^ 

dt dQ dyo dx 0 ^ 



(2.66a) 

(2.66b) 



(2.66c) 
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where the equation for the q ,• remains unaltered. In accounting for the terms in these 
equations, the partials 3 T /Sxq and 3 T /3yo will be recognized as the general momenta, 
and dT/dd as the generalized angular momentum. The linear momentum can then 
be represented by the expression 

P — (3r/3i 0 )i+ (3r/3yo)j. (2.67) 

Since the force equation is the rate of change of the linear momentum, we have 

F — [dp/dt] + to xp. (2.68) 



The terms of (2.66a) and (2.66b) are immediately accounted for. Moreover, the terms 
of (2.66c) can be identified from the expression 



M 0 = 




(2.69) 



The term 3 T /3 6 is the angular momentum ho, and the remaining two terms are equal 
to (dRo/dt) x m i (dr,/dt), where dRo/dt — (dxo/dt)i + (dy 0 /dt)j. 

Example 3. A typical example illustrating the above principles will now be given. 
Specifically, we will work out Problem 2, p. 118, of reference [4], Consider a particle 
moving in a plane attracted toward the origin of coordinates with a force inversely 
proportional to the square of the distance from it. In plane polar coordinates (r, 0) 
one has 




From these expressions, we form the Lagrangian function as follows: 




Furthermore, using (2.64), we obtain the derivatives as 




which give for this equation of motion 




For the other equation in the variable 6, we have 




d6 

dt 



and — = 0, 
dd 




References 
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so that this equation of motion is 




Note that since 0 is not explicitly present in L, the derivative of L with respect to 
dO/dt is a constant. 
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Aerodynamic Forces and Coefficients 



3.1 Aerodynamic Forces Relative to the Wind Axis System 

In this section we will discuss briefly the aerodynamic forces acting on a missile. 
In particular, and unless otherwise specified, we will assume a skid-to-tum missile 
because this technique is used in the majority of both surface-to-air and air-to-air mis- 
sile applications (for more details see Section 3.3.2). However, the reader should be 
cautioned that both the aerodynamics and rigid-body dynamics are highly nonlinear. 
For a more in-depth discussion of these forces the reader is referred to [2], [6], and 
[8]. Generally, the magnitude of the forces and moments that act on an air vehicle 
depend on the combined effects of many different variables. Briefly, the parameters 
that govern the magnitude of aerodynamic forces and moments include the follow- 
ing: (1) configuration geometry, (2) angle of attack, (3) vehicle size, (4) free-stream 
velocity, (5) density of the undisturbed air, (6) Reynolds number (i.e., as it relates to 
viscous effects), and (7) Mach number (i.e., as it relates to compressibility effects). In 
order to correlate the data for various stream conditions and configurations, the mea- 
surements are usually presented in dimensionless form. In practice, however, flow 
phenomena such as boundary-layer separation, shock-wave/boundary-layer inter- 
action, and compressibility effects limit the range of flow conditions over which 
the dimensionless force and moment coefficients remain constant. In essence, the 
motion of the air around an aircraft or missile produces pressure and velocity 
variations, which produce the aerodynamic forces and moments. As discussed in 
Section 2.2, the forces acting on a missile in flight consist of aerodynamic, propul- 
sive (i.e., thrust), and gravitational forces. These forces can be resolved along 
the missile’s body-axis system (Xb, 17,. Z/,) and fixed to the missile’s center of 
gravity (eg). The reference axis system standardized in guided weapons is cen- 
tered on the eg and fixed in the body. Thus, any set of axes fixed in a rigid 
body is a body-fixed reference frame. Before we proceed with the present dis- 
cussion, some of the fundamental concepts and definitions of aerodynamics will 
be reviewed. These definitions and nomenclature will be given with reference 
to Figure 3.1. 
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Relative 

wind 



Lift F 




Angle of Attack (a): The angle between the relative wind and the chord line. 

Aerodynamic Center: The point on the chord of an airfoil about which the moment 
coefficient is practically constant for all angles of attack. Moreover, the aerody- 
namic center is that point along the chord where all changes in lift effectively 
takes place. Since the moment about the aerodynamic center is the product of a 
force (the lift acts at the center of pressure) and a lever arm (the distance from the 
aerodynamic center to the center of pressure), the center of pressure must move 
toward the aerodynamic center as the lift increases. 

Center of Gravity: The forces due to gravity are always present in an aircraft (or 
missile) and act at the center of gravity (eg). Since the centers of mass and gravity 
in an aircraft practically coincide, there is no external moment produced by gravity 
about the eg. The gravitational force acting upon an aircraft is commonly expressed 
in terms of the Earth axes (see also Section 3.2.1). 

Center of Pressure:* The point on the chord of an airfoil through which all of the 
aerodynamic forces act. The center of pressure (cp) in general will not be located 
at the center of gravity of the airfoil; thus a moment will be produced. 

Dynamic Pressure: The aerodynamic pressure appears frequently in the derivation 
of aerodynamic formulas. Dynamic pressure, denoted by the symbol q, is given 
by the expression q — \ pV 2 , where p is the air density, and V is the free-stream 
velocity. 

Center of Mass: The origin of the body axes is usually the mass center (cm). 

Relative Wind: Refers to the motion of air relative to an airfoil and is equal and 
opposite to the forward velocity of the air vehicle. 

Resultant Aerodynamic Force: The vector summation of all of the aerodynamic forces 
acting on the airfoil. Its point of application is at the center of pressure. 



*Note that in aircraft design, aerodynamicists call the center of pressure (cp) the aerody- 
namic center ( ac ). Therefore, cp and ac will be assumed here to denote the same thing. 
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Fig. 3.2. Aerodynamic forces and thrust acting on a missile. 



It is conventional in aerodynamics to resolve the sum of the normal (or pressure) forces 
and the tangential (or viscous shear) forces that act on the surface due to the fluid 
motion around a vehicle into three components along axes parallel and perpendicular 
to the free-stream direction. These forces are lift (L), drag (D), and side force (Y). 
The relation of the lift and the drag forces to the free-stream velocity is shown in 
Figure 3.2. It should be noted from this figure that if an angle of attack is generated, 
the lift vector acting at the center of pressure ( cp ) has a destabilizing effect and must 
be controlled. 

We will now define these forces in some detail. 

Lift-Lift is the component of the resultant aerodynamic force that is perpendicular 
(i.e., upward) to the relative wind (direction of flight) or to the undisturbed free- 
stream velocity. The aerodynamic lift is produced primarily by the pressure forces 
acting on the vehicle surface. Also, the lift force is perpendicular to the missile’s 
velocity vector in the vertical plane. 

Drag-Drag is the component of the resultant aerodynamic force that is parallel to the 
relative wind. In other words, it is net aerodynamic force acting in the same direc- 
tion as the undisturbed free-stream velocity. The aerodynamic drag is produced 
by the pressure forces and by skin friction forces that act on the surface. The drag 
force is measured along the velocity vector, but in the opposite direction. 

Side Force-Side force is the component of force in a direction perpendicular to both 
the lift and the drag and is measured in the horizontal plane. The side force is 
positive when acting toward the starboard wing, provided that the bank angle is 
zero. If the bank angle is not zero, L and Y will be rotated by a negative angle 
about the velocity vector. 

The definitions of the aerodynamic forces, moments, and velocity components in the 
body-fixed coordinate system, which will be used in this chapter, are summarized 
in Table 3.1. 
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Table 3.1. Missile Aerodynamic Moments, Coordinates, and Velocity Components 





Roll Body 

Axis (Xb) 


Pitch Body 
Axis (Y b ) 


Yaw Body 
Axis (Zb) 


Angular Rates 


P 


Q 


R 


Velocity Components 


u 


V 


w 


Aerodynamic Force 


Fx 


Fy 


F Z 


Components 


Aerodynamic Force 


Cd 


Cy 


Cl 


Coefficients 


Aerodynamic Moment 


Ci 


Cm 


C n 


Coefficients 



The basic aerodynamic forces are commonly defined in terms of dimensionless 
coefficients, the flight dynamic pressure, and a reference area. For missiles that skid 
(yaw) to turn (see Section 3.3.2 for more details), the basic aerodynamic forces are 
illustrated in Figure 3.2 and are calculated as follows [2], [6]: 

Drag: D = CoqS, 

Lift: L — CiqS, 

Side Force: Fy — CyqS, 

— Coefficient of drag in the wind axis system, 

= Coefficient of lift in the wind axis system, 

= Side force coefficient, 

= Free-stream dynamic pressure at a point far from the airfoil = y pV 2 , 

= Reference area, usually the area of one of the airfoils, 

= Free-stream velocity, 

n ^ — 4 

= Atmospheric density = 2.3769 x 10 lb-sec~-ft 
at sea level (see also Appendix D). 

For missiles that roll to turn, drag is the same as in (3.1), but the lift and side force 
are as follows: 



where 

Cd 

Cl 

Cy 

q 

s 

v 

p 



(3.1) 

(3.2) 

(3.3) 



Lift= Clt ( cos (p)qS, 
Side Force= C^tsin cp)qS , 



(3.4) 

(3.5) 



Clt — Total lift coefficient in the maneuver plane = (Cj + Cy ) l 22 , 
cp — Roll angle. 



where 
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For the purposes of the discussion in this book, the three most important aerodynamic 
force coefficients are commonly defined as* 

C L = L/qS, 

Cd — D/qS , 

Cm — M/qSd, 

where M is the moment and d is the mean missile diameter from a body cross-section. 
Figure 3.3 illustrates the aerodynamic forces relative to the wind-axis system. 

The aerodynamic forces may also be expressed in the form 

Axial Force (Drag): Fx = q(V, h)SCo(V, h, a, ft ), (3.6) 

Side Force: Fy = q(V, h)SCy(V , /;, a, /3), (3.7) 

Normal Force (Lift): Fz = q(V, Ii)SCl(V, h, a , /3), (3.8) 

showing the dependence on the angle of attack (a), sideslip angle (/?), and altitude 
h. On occasion, it may be convenient to measure the aerodynamic forces in the 
body axis coordinate system. In this case, we have the normal force (Fz) along the 
Zfc-axis, side force (Fy ) along the K/,-axis, and the axial force (Fx) along the A/, -axis. 
The specification of forces in the body-axis system is similar to that in the wind-axis 
system. That is, 



Fxb — qSCo , 

Fyb — qSCy, 

Fzb — qSCL- 

The aerodynamic force coefficients Cl, Cd, and Cy are commonly expressed in the 
wind-axis system oriented relative to the free-stream. Since the aerodynamic force 
components of the equations of motion (see Section 2.2) are required to be in the body- 
fixed coordinate system, one must express these coefficients in terms of the angle of 
attack and sideslip angle. The aerodynamic force coefficients can be determined in the 
wind tunnel in the body-fixed axis system, designated as Cxb , Cyb, and Czb- Thus, 

Cxb — —Cd cosa cos f — Cy cosa sin f + Cl sina, (3.9) 

Cyb — — Cd sin/3 + Cy cos /l, (3.10) 

Czb — —Cd sin a cos /3 — Cy sin a sin f — Cl cos a. (3.11) 

Figure 3.4 illustrates these coefficients. 

For a simple point mass case, relative to the airstream, the aerodynamic force 
coefficients Cd, Cl, and Cy will be assumed to be functions of one or more of 
the following: 

*Note that here we assume that the specification of forces in the body-axis system is similar 
to that in the wind-axis system. 
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Earth axis system 



Thrust (7) 




(1) Angle of attack and sideslip angle. 

(2) Lift and/or side force. 

(3) Mach number and/or Reynolds number (plays a role only in the drag force). 

(4) Center of gravity location. 

(5) Altitude. 

As mentioned above, all the aerodynamic force coefficients are, in general, functions 
of the state variables and the control variables, so that one can write, for example, 



Cd = Cd(u, P, M,q, 



3.1 Aerodynamic Forces Relative to the Wind Axis System 



59 




Fig. 3.4. Wind tunnel representation of the aerodynamic coefficients. 



By taking the partial derivatives, we have 

clC d = (9C o /da)da + (dC d / df)df) (dCo/dM)dM (dCo/dq)dq + ■ ■ ■ . 

Therefore, the aerodynamic force coefficients Cd , Cl, and Cy may be expressed in 
terms of the aerodynamic derivatives as follows: 

Cd = Cdq + Cdu\u\ + C Da icr + Cop \P\ + C D pifl 2 + Cdup\u\\P\, 

Cl = Clo + CLa |a| + C La icr + Clp\P\ + C Ljg 2 ^ 2 + CL a p\u\\P\, 

Cy = Cyo + Cy a \a\ + C Yol 20 T + Cyp\fi\ + Cypp 2 + Cy a p \a \ \fi \ , 

where Cdq — (9CD/9°Ola=o 0- e -> evaluated at a — 0), Cr> a = dCo/da, etc. 
purposes, the functional dependence of the aerodynamic force coefficients 
assumed to take the simpler form as follows [6], [8]: 

Drag Coefficient (Cd) 

Cd — Cdo + CoaOt, 

where 

Cdo — total drag coefficient evaluated at a — 0 (or close to it) = (9Co/9tf )|„ = o, 
Coct — total drag coefficient variation with angle of attack = DC o /'da, 
a — angle of attack (in radians). 

The derivatives are evaluated at constant Mach number and Reynold’s number. The 
drag polar is written in the form [4] 

Cd — Cdo + KC 2 , 



(3.12) 

(3.13) 

(3.14) 

For our 
will be 



(3.15a) 



(3.15b) 
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where 

Cdo — zero lift drag coefficient, 

K — drag due to lift factor (also called the separation drag due to lift factor) 

= dC d / dC~^, 

Cl — lift coefficient. 

Equation (3.15b) states that the total drag may be written as the sum of (1) the drag 
that exists when the configuration generates zero lift (Cdo), and (2) the induced drag 
associated with lift (KCj). 

Lift Coefficient Cl 

Cl — Clo + CLad, (3.16a) 

where 



Clo — total lift coefficient evaluated at a — 0 
= OC L /9«) | a=0 , 

Cloi — total lift-curve slope. 

The derivatives here are evaluated at constant Mach number. The lift coefficient can 
also be written as 

CL = (dCL/da)\ a= oa + C,a 2 , (3.16b) 

where C\ is a nonlinear factor. 



Side Force Coefficient (Cy) 

The functional dependence of the side force coefficient on sideslip angle, ft, aileron 
angle, Sa, etc., is expressed as 

Cy = Cy 0 + Cypfi + Cy$&A, (3.17) 

where 



Cy 0 = side force coefficient for zero sideslip and zero control deflection 
= (dCy/df)\p = o, 

Cyp = change in side force coefficient due to a unit sideslip angle 

= dCy/dffi 

f> — sideslip angle (in radians). 



The derivatives here are evaluated at constant Mach number and constant angle of 
attack. 

The components of the normalized instantaneous accelerations in the wind-axis 
system are calculated as follows (see Figure 3.3): 



Ax = (T cos(a + £) cos f — Coq S)/ W, 



(3.18a) 
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Ay = (T cos(a + £) sin/1 + CygSVW, (3.18b) 

A z = (T sinter + ft) + C L qS)/ W, (3.18c) 



where 

W — the missile weight, 

T — the thrust, 

a — angle of attack in the pitch plane, 
ft — sideslip angle, 

0 = missile pitch reference angle = a + y, 
xfr = missile yaw reference angle = ft + r, 
y = flight path angle in the vertical plane, 

ft = thrust inclination relative to the missile body axis in the pitch plane, 
r = flight path angle in the horizontal (X p. Y p) plane. 



The instantaneous accelerations in the Earth-axis system (see Figure 3.3) are obtained 
through the following transformation [6] : 



d 2 X p 
dt 2 



— (Ax cos y cos r — Ay sin r — Az sin y cos r)g, 



d 2 Y E 
dt 2 



= (Ax cos y sin r + Ay cos r — Az sin y sin r)g, 



d 2 Zp 
dt 2 



= ({-A x sin y - A z cos y)g. 



(3.19a) 

(3.19b) 

(3.19c) 



where g is the gravitational acceleration and y is the flight path angle in the verti- 
cal plane. Integration of (3. 19a)— (3. 19c) yield the velocities dX E /dt , dY E /dt , and 
dZ E /dt. The velocities can then be integrated to obtain the missile position coordi- 
nates X E , Yp, and Z E . The angle of attack (a) and sideslip (ft ) can be defined in terms 
of the velocity components as shown in Figure 3.3. Mathematically, the equations for 
these angles are given in the form 

a = tan ~ 1 (w/u), (3.20a) 

ft — sin -1 (u/ V m), (3.20b) 



where Vm = (« 2 + v 2 + w 2 ) 1 ^ 2 . If the angle of attack and sideslip are small, say, 
< 15°, then (3.20a) and (3.20b) assume the simpler form 



a = w/u, (3.21a) 

ft = v/u, (3.21b) 



where a and ft are given in radians. Angle of attack and sideslip completely define 
the attitude of the vehicle with respect to the velocity vector. These angles can also be 
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expressed in terms of the other angles indicated in Figure 3.5 (see also Figure 2.1(b)) 
as follows: 



tan a = tan 0 cos tf> , 


(3.22a) 


tan y = tan 0 sin tp, 


(3.22b) 


sin /l = sin 0 sin 4> , 


(3.22c) 


cos 0 = cos a cos y6 = (tan 2 a + tan 2 yS) 1 / 2 , 


(3.22d) 


tan </> = cot a tan y — tan /3 /sin a. 


(3.22e) 



3.2 Aerodynamic Moment Representation 

In a similar manner to the aerodynamic forces of the previous section, the moments on 
the missile can be divided into moments created by the aerodynamic load distribution 
and the thrust force not acting through the center of gravity. Specifically, the moment 
due to the resultant force acting at a distance from the origin may be divided into 
three components, referring to the missile’s body reference axes. The three moment 
components are the pitching moment , the rolling moment, and the yawing moment. 
These moments will now be defined more closely: 

Pitching Moment: The pitching moment is the moment about the missile’s lateral 
axis (i.e., the F/,-axisj. The pitching moment is the result of the lift and the drag 
forces acting on the vehicle. A positive moment is in the nose-up direction. 
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Rolling Moment: This is the moment about the longitudinal axis of the missile 
(i.e., the A/, -axis). A rolling moment is often created by a differential lift, generated 
by some type of aileron. A positive rolling moment causes the right or starboard 
wingtip to move downward. 

Yawing Moment: The moment about the vertical axis of the missile (i.e,, the Z*- axis) 
is the yawing moment. A positive yawing moment tends to rotate the nose to the 
right. 

It should be pointed out here that the calculation of the aerodynamic forces and 
moments acting on a vehicle often requires that the engineer be able to relate data 
obtained at other flow conditions to the conditions of interest. For example, the 
engineer often uses data from the wind tunnel, where scale models are exposed to flow 
conditions that simulate the design environment or data from flight tests at other flow 
conditions. In order that one can correlate the data for various stream conditions and 
configuration scales, the measurements are usually presented in dimensionless form. 
The procedure used to nondimensionalize the moments created by the aerodynamic 
forces is similar to that used to nondimensionalize the lift. Thus, the components of 
the aerodynamic moment can be expressed in terms of dimensionless coefficients, 
flight dynamic pressure, reference area, and a characteristic length as follows [6]: 

Rolling Moment (L): L — CiqSl, (3.23a) 

Pitching Moment (M): M = C m qSl, (3.23b) 

Yawing Moment {N)\ N = C n qSl, (3.23c) 

where C/, C m , and C n are the aerodynamic moment coefficients in roll, pitch, and 
yaw, respectively. Note that in missiles, the reference area S is usually taken as the 
maximum cross-sectional area, and the characteristic length / is taken as the mean 
diameter, whereas in aircraft [2], [6], 

L = C,qSb , (3.24a) 

M — C m qSc, (3.24b) 

N — C n qSb, (3.24c) 

where b is the wingspan, c is the aerodynamic chord, and S is the wing planform area 
used to nondimensionalize the aerodynamic forces. In general, and as stated earlier, 
the standard six-degree-of-freedom aerodynamic coefficients Cl, Cd, Cy , C/, C m , 
and C n are primarily a function of center-of-gravity location, altitude, Mach number, 
Reynolds* number, angle of attack (a), and sideslip angle (fi), and are 
secondary functions of the time rate of change of angle of attack and sideslip, and the 
angular velocity of the missile. (The pitching moment coefficient C m is independent 
of the Reynolds number). We will now develop the aerodynamic moments and their 
associated coefficients in terms of measured quantities. 



*The Reynolds number is a nondimensional number defined as R = (pVl)/ p = (Vl)/v, 
where p is the density of the fluid, p is the coefficient of absolute viscosity of the fluid, v is 
the velocity, l is the characteristic length, and v is the kinematic viscosity (v = p/p). 
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(a) Aerodynamic Rolling Moment 

We begin the development of the moments by noting that the rolling moment L is a 
function of a p , a y , and S a , where 

oip — angle of attack in the pitch plane; in the present development, the angle of 
attack will be taken to be the acute angle between the A/, -axis of the missile 
and the line of relative airflow or the missile velocity; 
a y = angle of attack in the yaw plane (this angle is identified as the sideslip 
8 a = deflection angle of the aileron (this angle controls the roll of the missile). 

The moments L , M, and N can be linearized using the Maclaurin series for three- 
variable functions. The Maclaurin series is a type of Taylor series, and so a function 
in three variables can be approximated by the sum of the first three terms of the series. 
As a result, the rolling moment L can be approximated as follows: 

L = [(dL/doip) \a p =o\Up + [( dL/da y )\ ay =o\ciy + [(3L/3^ fl )U Q , = o]5 a = 0. (3.25) 

Let us now define the following coefficients: 

[(3L/3a p )| ap= o] = L Up = the rate of change of the rolling moment due to a 

change in angle of attack in pitch, 

[(3L/3a y )| a>I= o] = L a y — the rate of change of the rolling moment due to a 

change in angle of attack in yaw, 

[(3L/3a y )|^ = o] = L$ a — the rate of change of the rolling moment due to a 
change in aileron deflection angle. 

In order to simplify equations and calculations, the above coefficients can be non- 
dimensionalized so that they become normalized moments: 

Lot p /qSl re f — c, ap — rolling moment coefficient due to angle of attack in pitch, 

L ay /q SI re f — C j = rolling moment coefficient due to angle of attack in yaw, 

L sJqSl ref — C ha — rolling moment coefficient due to aileron deflection angle, 

where 



q = free-stream dynamic pressure at a point far from the airfoil, 
S — reference area, usually the area of one of the airfoils, 

I re f — reference length, usually the mean missile diameter from 
a body cross-section. 

Figure 3.6 illustrates the angles under discussion. 
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Y aw plane 



Fig. 3.6. Missile angle representation. 



The rolling moment L can now be rewritten using the above-introduced coeffi- 
cients and definitions: 

L — [(dL/da p )\ ap= o]a p + [(dL /da y )\ ay -{)]oiy + [(3L/95 a )|fl 0= o]5 a 
— (Lct p )a p + ( L ay )oiy + (L$ a )8 a 

= qSl ref {C lap )u p + qSl ref (Ci ay )a y +qSl ref (C ha )a p 
= qSl ref m ap )a p + ( Ci ay )a y + (Ci s J8 a ]. (3.26) 

Solving (2.29a) for dP /dt, we can now write the rotational acceleration equation for 
roll in the form 

Cl £ = QR[(Iy - I z )/I x ] + 0 qSlref/I x )m ap )a p + (C lay )a y + (C ka )8 a ]. (3.27) 

Before we proceed with the derivation of the aerodynamic pitching and yawing 
moments, it will be necessary to develop the nondimensionalized aerodynamic nor- 
mal force coefficients, as well as the nondimensionalized pitching and yawing moment 
coefficients. We begin by noting that when an inclined surface moves through the air, 
there is a force perpendicular to that surface, caused by the deflecting stream. This 
force is called the aerodynamic normal force, and is normal (perpendicular) to the 
Xfe-body axis. Now let us define the following quantities: 

Fy = component of the aerodynamic normal force (side) along the T/,- body axis, 
being positive from the origin in the direction of the negative K/,-axis, 

Fn = component of the aerodynamic normal force (normal) along the Z/,-body axis, 
positive from the origin in the direction of the negative Z/,-axis. 

These two components of the aerodynamic normal force can be non-dimensionalized 
as follows: 



Cy = Fy/qS, 

Cn = Fn /q S , 



(3.28a) 

(3.28b) 
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where Cy is called the lateral force coefficient, and Cm is called the normal force 
coefficient. The components Cy and Cn are usually defined as positive for positive 
angles of attack and zero control surface deflection (see Figure 3.7). These aerody- 
namic normal force coefficients are extremely nonlinear and usually cannot be accu- 
rately linearized, as could the rolling moment coefficient. The aerodynamic normal 
force coefficients are functions of Mach number (M), a p or a y , and 8 p or 8 y , and are 
commonly written in the form 

Cy(M, a y , <5 V ), 

Cn(M , a p , 8 p ), 

where 8 p is the control surface deflection in pitch, and <$ v is the control surface 
deflection in yaw. 

(b) Aerodynamic Pitching Moment 

The aerodynamic pitching moment M is a function of the pitch rate Q as well as a p 
and S p . Define now the components of the pitching moment as 

M = M„ + M q , (3.29) 

where 

M 0 = moment contribution from the AO A (angle of attack) in pitch (a p ) 
and pitch control surface deflection ( 8 p ), 

M q — pitching moment rate. 

These two components are commonly described in the following manner: 

M„ = C m qSl re f, (3.30a) 

M q = (C M QqS(lref) 2 )/ 2Vm, (3.30b) 

where 

C m = pitching moment coefficient, which is a function of Mach 
number, u p . and 8 p , 

Cm — coefficient of moment due to pitch rate, or rate of change 
of pitching moment, 

Vm — total velocity of the missile. 



(c) Aerodynamic Yawing Moment 

The aerodynamic yawing moment A is a function of yaw rate R as well as a y and 
8 y . Now, define the components of the yawing moment as follows: 

N = N a + N r , 



(3.31) 
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8 = Positive deflection of control surface, measured with respect to X^-axis; 
positive for anti-clockwise sense of rotation about the hinge looking into 
the F^-axis towards the missile body. 

Tail - Positive deflection as shown is producing: 

1) Negative Pitching Moment (DM) about the eg. 

2) Positive Normal Force (+F N ). 

Canard - Positive deflection as shown producing: 

1) Positive Pitching Moment (+M) about the eg. 

2) Positive Normal force (+F N ). 



(a) X-Z (pitch) plane. 




8 = Positive for clockwise sense of rotation about the hinge looking down 
the Y b - axis away from the missile body. 

Tail - Positive deflection as shown is producing: 

1) Negative Yawing Moment (DAO about the eg. 

2) Negative Lateral Force (DFy). 

Canard - Positive deflection as shown producing: 

1) Positive Yawing Moment (+N) about the eg. 

2) Negative Lateral force (DFy). 

(b) X-Y (yaw) plane. 

Fig. 3.7. Control surface deflections. 



where 



N 0 = the moment contribution from the AO A in yaw (a y ) 
and yaw plane control surface deflection (<5 V ), 

N r = yawing moment rate. 
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F = mass x vertical 
acceleration 




Free-stream 

velocity 



Fig. 3.8. Forces acting on a missile in the pitch plane. 



These two components are usually described in the following manner: 

N 0 = (■// (j Sl re f , 

N r = (C N RqS(l re f) 2 )/ 2V M , 

where 



(3.32a) 

(3.32b) 



C„ = yawing moment coefficient, which is a function of Mach number, 
a y , and<5y, 

C,v = coefficient of moment due to yaw rate. 

In missiles, the center of gravity (eg) normally shifts due to the burning off of fuel. 
Consequently, a new moment is produced. Therefore, M a and N 0 are evaluated at a 
reference center of gravity at launch. Also, the actual center of gravity moves when 
the missile begins to move. However, the referenced center of gravity remains in the 
original position. Thus, applying this concept of a referenced center of gravity, the 
moment contributions due to the movement of the actual center of gravity can be 
evaluated. A simple free-body diagram that shows the forces and moments acting on 
the missile in the pitch plane is illustrated in Figure 3.8. 

When the moment coefficients are evaluated, the vector normal force is 
assumed to act through the referenced center of gravity (cg re f). Moreover, the missile 
vertical acceleration vector (F /in) acts through the same point, which corresponds to 
the momentary center of gravity at launch. As the real center of gravity moves, the 
acceleration vector stays at the same location, but an additional moment arises due to 
the motion of the center of gravity. Since the force is taken to be applied through the 
reference center of gravity, the induced moment is clockwise about the actual center 
of gravity, is negative, and is equal to 

(— )M C g = F]\/(dcg—ref — d C g), 



(3.33) 
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where M cg is the moment of about the actual center of gravity and ( d cg - re f — d cg ) 
is the distance between the actual center of gravity and eg (note that d is usually 
measured from the nose of the missile). Equation (3.33) can also be written in the 
form 



M cg — Fn ( d cg d cg —ref ) — ( \ t/ S(d l g d cg —ref ) • (3.34) 

It should be pointed out that the force F would produce no moment about the center 
of gravity, since it acts through that point. An induced moment will similarly be 
produced in the yaw plane with Fy and lateral acceleration. A positive aerodynamic 
force will produce a positive moment: 

N C g — Fy(d cg -ref d cg ) , (3.35a) 

or 

N cg = Fridcg d cg ) — CyqSidcg d cg —ref )• (3.35b) 

The above equations for the pitching moment and yawing moment must be rewritten 
in order to include the change in the center of gravity position. Thus, 

PitchingMoment : M — M a + M q + M cg , 

YawingMoment : N — N 0 + N r + N cg . 

From these expressions, we can now obtain the pitching moment (M) and yawing 
moment (N) by substitution as follows: 

M — C m qSl re f + CMQqSQref) )/2 Vm CyqSidcg d cg —ref ) 

= q SI ref {C m T" C M ((Qlref )/2 V'v/ ) + Cm[ (dcg d C g—ref)/ lref]}> (3.36) 
N = C n q Slyef + Cn Rq S(l re f )~ /2 Vm Cy q S (d cg dcg—ref ) 

= qSlref{C n + CN((Rlref) /2 Vm) - Cy\(d cg - dcg-ref) / Iref]}- (3.37) 

Substituting (3.25) and (3.26) into (2.30b) and (2.30c), we obtain 

= (Uz ~ h)/Iy)P R + (<7 Sl re f / 1 yy){C M + Cl\[{(d cg — d C g—ref) / Iref) 

+ C M (Qlref/ 2 Vm)} (3.38) 

— ((/x - Iy)/h)PQ + (q Slref /Iyy){C N ~ Cy{{d cg ~ dcg-ref) / Wef) 

+ C N (Rlref/ 2 Vm)}. (3.39) 

Equations (3.38) and (3.39) are the rotational accelerations for pitch and yaw, 
respectively. The following definitions from second-order differential equations 
should also be noted: 

PitchDampingMoment : Cm — 9Cm /<HQ l>-ef / 2 Vm), 

YawDampingMoment: Cn — dC n / d(Rl re f / 2Vm) ■ 
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These damping moments are important to the missile designer, in that they are 
necessary to keep the missile from oscillating. The rotational rates P, Q. and R 
can be obtained by integrating the equations for d P fell , dQ/dt , and dR/dt. Thus, 




(3.40a) 

(3.40b) 

(3.40c) 



In terms of the Euler angles cp (roll), (/(pitch ), and V/fyaw), the rotational rates P, Q, 
and R (assuming that the missile A/, -axis is along the longitudinal axis, the K/,-axis 
out to the right, and the Z*- axis down) can be expressed in the form [1] 



P = 



dtp 

dt 




(3 .4 Id) 



Q = 




COS 4 > + 




cos 9 sirup. 



(3.41e) 



R — cos 9 cos tp — sin tp , (3.41f) 

or in matrix form. 



p~ 




'1 0 


— sinO 


0 


Q 


= 


0 cos cp 


cos 6 sin cp 


9 


R 




0 — sin (p 


cos 6 cos (p 


\p 



(d) Derivation of the Translation Equations (A* -Body Axis) 

The translational equations of motion can be derived from (2.11). Therefore, rear- 
ranging (2.1 1), we have 



— Rv - Qw + ( F x /m ), 
dt 

— = Pw — Ru + ( Fy/m ), 
dt 

dw 

— = Qu - Pv + (F-Jm). 
dt 



(3.42a) 

(3.42b) 

(3.42c) 



Next, we must determine the forces acting along each of the body axes. Assume now 
that a flat surface moves through a mass of air at some angle of attack. Friction acting 
between the mass of air and the inclined surface, resulting in a force that tends to move 
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the surface parallel to itself. There is also a force perpendicular to the surface, caused 
by the pressure difference across the plate. These combined forces give a resultant 
force backward in the direction of the moving air mass. Resolving this resultant force 
into components, we have one force with its vector perpendicular to the moving mass 
of air, called the lift, and another force parallel to the mass of air, called the drag. 
Figure 3.9 illustrates these forces. 

This figure refers to one panel of the missile airfoil.* However, the concept can 
be extended to illustrate the effects on a full airfoil. At each panel there is a pressure 
force acting normal to the panel and a friction force acting parallel to the panel. By 
summing these forces and taking the component parallel to the velocity vector and 
component perpendicular to the velocity vector, the lift and drag can be obtained. Drag 
acts along the total velocity vector, but the largest component is along the A/, -axis. The 
A/, -component of drag is referred to as the axial force Fa, and as is customary in 
aerodynamics, this force is described in terms of a nondimensionalized coefficient 
Ca • Thus, 



F A =qSC A , (3.43) 

where Ca is the axial force coefficient, which is a function of Mach number, S p , S y , 
and altitude (with respect to the earth) of the missile body. The other primary force 
acting along the A/, -body axis is the thrust. Thrust is defined as the forward force 
produced by the propulsion system to sustain the aircraft in flight. Sometimes, thrust 
is expressed in the same manner as drag and other aerodynamic forces. However, it is 
not usually convenient to do so, since the thrust is often constant or is some unknown 
function of altitude, whereas the other aerodynamic forces are not. Typically, a missile 
has two modes of thrusting: (1) the boost phase, and (2) the sustain phase (a boost- 
sustain mode is also used. See Sections 3.3.1 and 4.5 (Table 4.5)). The boost phase is 
the first stage of missile flight. During this phase, the missile is propelled from rest to 
slightly supersonic speeds by a high-powered rocket engine. The sustain phase is the 



* An airfoil is a streamlined body that when set at a suitable angle of attack, produces more 
lift than drag. 
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second stage of missile propulsion, which begins when the missile is traveling close to 
the speed of sound (i.e., Mach 1 ). The missile booster is an auxiliary propulsion system 
that imparts thrust to the missile during the initial phase of flight. It generally consists 
of a solid rocket motor and an attaching device. A solid rocket motor is best suited for 
this purpose because of its simple construction and operation and its ability to develop 
a high thrust in a short time. The missile sustainer, on the other hand, is usually a 
ramjet engine, which allows a greater missile range than would a rocket engine. Often, 
the booster and sustainer burn simultaneously for a short time prior to jettison of the 
booster, in order to ensure a reasonably smooth transition to the sustain phase. 

Prior to the sustain phase, the missile sustainer is open, allowing airflow that will 
exert a negative thrust (i.e., an axial drag force of (—C A qS), since F A = C A qS). 
Typically, the total thrust is expressed as 

T — 7b oos t + r sus t , (3.44) 

where T is the total thrust, 7j, 00st is the thrust from the booster, and 7’ sust is the thrust 
from the sustainer. The remaining force acting along the A/, -ax is of the missile is 
the A-component of weight due to gravity. Referring to Figure 3.10, we note that 
the missile body and inertial Earth axes will show that the weight component of the 
missile along the A/, -axis is 



x mg = mg sin0, 



(3.45) 



where 

x mg = weight component along the A/, -axis, 

Z mg = weight component along the Z/,- axis, 
g = acceleration due to gravity, 

9 — angle between the body axes and Earth axes (pitch angle), defined as 
positive for the counterclockwise sense of rotation in the A — Z plane. 



Summing all the force components along the A/, -axis, we have 



Fx — Tboost T" Tsiist x mg 

— Iboost + 7’sust - mg sin 9 - C A qS. 



(3.46) 



(e) Derivation of Translation Equations (F/,-Body Axis) 

The forces acting along the K/,-body axis of a missile consist of the K/, -components 
of weight and aerodynamic normal force. Recall that the K/, -component of the aero- 
dynamic normal force is 

Fy = CyqS. 

Figure 3.11 illustrates the weight component of the missile along the K/,-axis in three 
dimensions with respect to the missile-body axes and inertial Earth axes. 
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Horizon 



Fig. 3.10. Missile body and Earth-fixed axes. 




Fig. 3.11. Three-dimensional view of the missile weight components. 



The weight component of the missile along the T^-axis is 



where 



y mg = mg sin cp cos 9, 



y mg — weight component along the JVaxis, 



(3.47) 



4> — angle between the projection of the weight vector onto the F/, — Z/, plane 
and the Zb- axis; this is the roll angle and is defined to be positive for the 
counterclockwise sense of rotation in the Yb — Zb- 

Note that the weight vector is directed along the Ze- axis. If we project this vector 
onto the Yb — Z], plane by dropping a perpendicular, we obtain the vector mg cos 6. 
Next, we drop a perpendicular from this vector onto the Zb- axis, thus obtaining the 
Z/, -component of the weight: 



Zmg — mg cos 9 cos cp . 



(3.48) 
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The third side of this vector triangle is parallel to the i^-axis, and so is the 
K/, -component of the weight. This vector is positive because of its orientation with 
respect to the K/,-axis. Now we can sum the force components along the K*-axis: 

F y — y mg — Fy — mg sin f cos 6 — CyqS. (3.49) 



(f) Derivation of Translation Equations (Z/,-Bodv Axis) 

The forces acting along the Z/,-body axis of a missile consist of the Z/, -component 
of weight and the Z/, -component of the aerodynamic normal force. We have noted 
earlier that the Z/, -component of the aerodynamic normal force is 

Fn — CnqS 

and that the Z/, -component of the weight is 

Zmg — m g cos 4> cos 

Summing the force components along the Z/,- axis, we have 

F z = z m g ~ Fn = mg cos </> cos 6 — CNqS. 

Substituting the above results into (3.42) results in the following equations: 

C ^~ = Rv - Qw + (F x /m) = Rv — Qw + [(Tboost - T sust )/m] 
at 

—g sin0 — ( CAqS/m ), 
civ 

— = Pw — Ru + (Fy/m) — Pw — Ru + g cos 6 sin (p — ( CyqS/m ), 
dt 

dw 

— = Qu — Pv + (F z /m) = Qu — Pv + geos </> cos 6 — ( CNqS/m ). 
dt 



(3.50) 

(3.51a) 

(3.51b) 

(3.51c) 



Finally, in order to obtain the translational velocity components n, v, and w, (3.51) 
must be integrated. Thus, 

Longitudinal Component of Velocity 




Lateral Component of Velocity 




Vertical Component of Velocity 




(3.52a) 



(3.52b) 



(3.52c) 



We will now summarize the translational and rotational equations of motion of a 
missile in terms of the body axes: 




Translation Equations 

(1) Longitudinal Acceleration 
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— = Rv - Qw + (Tb 00st - r sust )/m -g sin 0 - ( C A qS)/m . 
at 



(2) Lateral Acceleration 

dv 

— — Pw — Ru + g cos 6 sin <p — ( CyqS)/m . 
dt 

(3) Vertical Acceleration 

dw 

— = Qu — Pv + geos 9 cos0 — ( CuqS)/m . 
dt 

Rotation Equations 



(3.53a) 



(3.53b) 



(3.53c) 



(4) Roll Acceleration 

^ = (qSl ref /I x )[C L (8 a ) + C L (a p ) + C L (a y )] + QR[(I y - I z )/I x ]. 

(3.53d) 

(5) Pitch Acceleration 

dQ 

— — (qSl ref /I y )[C M + C N ((d cg - d cg - ref )/l ref ) + C M (Qlref/2V M )] 

+ PR[(I Z -I x )/I y l (3.53e) 

(6) Yaw Acceleration 

— = (qSl ref /I z )[C N — CY((d cg — d cg - re f)/ l re f) + Csj{Rl re f )2V m)\ 

+ PQ[(I x -Iy)/hl (3.53f) 

Example 1. In Section 2.1, the transformation matrix from the Earth to body axes was 
developed. Consider now the free-flight dynamic model of a missile. The mathemati- 
cal model describing the missile motion consists of six rigid-body degrees of freedom 
(i.e,, three body inertial position coordinates and three Euler-angle body attitudes). 
In this example, we will use the Earth’s surface (or ground) as the inertial refer- 
ence frame. The body frame is defined in the conventional manner, and the dynamic 
equations are written with respect to this coordinate system [9]. The missile’s 
translation and rotation kinematic and dynamic equations are given by 



X 




CQC-Is S(f)SQC\ls C(j)S\Js CtySQCy}; S(j)SiJs 
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where 
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c — cos, 
s — sin, 
t — tan, 

m — mass of the missile, 

L — rolling moment, 

M — pitching moment, 

N — yawing moment, 

P — roll rate, 

Q = pitch rate, 

R — yaw rate, 
cp — roll angle, 

6 — pitch angle, 
x[r — yaw angle, 

u,v,w — components of velocity of the center of mass 
relative to the atmosphere. 



( 2 ) 

( 3 ) 

( 4 ) 



The total applied force is composed of the weight W and body aerodynamic force A 
terms. The weight portion of the external loads is given by 



x w 




~sg 


Yw 


= mg 


S<j>Cg 


_Zir_ 




_ C(j>Cg _ 



where g is the gravitational acceleration. The aerodynamic force contribution is 
given by 
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( 6 ) 
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where 



a = tan _1 (w/i>), — n < a < n, (7a) 

ji = tan _1 (i>/«), — 7r < a < jt, (7b) 

q a = \p(u 2 + v 2 + w 2 )tt D 2 , (8) 

where D is the drag. The right-hand side of the rotation kinematic equations contains 
the externally applied moments. Note that the external moment components contain 
contributions from steady and unsteady body aerodynamics. The steady-body aero- 
dynamic moment is computed by a cross product between the distance vector from 
the center of gravity to the center of pressure and the steady-body aerodynamic force 
vector given above. As in the case of the aerodynamic coefficients, the center of 
pressure location is dependent on the local Mach number and is computed by linear 
interpolation. 

In Section 3.1 we discussed briefly the role of the airfoil. The airfoil used in 
modern airplanes has a profile of the “fish” type, as illustrated in Figure 3.1. Such an 
airfoil has a blunt leading edge and a sharp trailing edge. The projection of the profile 
on the double tangent, as shown in the diagram, is the chord. The ratio of the span 
to the chord is the aspect ratio. The camber line of a profile is the locus of the point 
midway between the points in which an ordinate perpendicular to the chord meets 
the profile. The camber is the ratio of the maximum ordinate of the camber line to the 
chord. From the theory of the flow around such an airfoil, the following assumptions 
apply: 

(1) That the air behaves as an incompressible inviscid fluid (see also Section 2.2). 

(2) That the airfoil is a cylinder whose cross-section is a curve of the above type. 

(3) That the flow is two-dimensional irrotational cyclic motion. 

The above assumptions are, of course, only approximations to the actual state of 
affairs, but by making these simplifications, it is possible to arrive at a general under- 
standing of the principles involved. 



3.2.1 Airframe Characteristics and Criteria 

In this section we will examine the general airframe features and stability from the 
point of view of the guidance designer. Commonly, the airframe is symmetrically 
cruciform, with four fixed wings and four movable control surfaces. The cruciform 
configuration permits lateral maneuvering in any direction without first rolling (i.e., 
banking to turn, as required of an airplane). In a wing-controlled airframe with 
movable wings slightly forward of the center of gravity (eg) and fixed stabilizing 
tail surfaces, variable downwash from the control surfaces impinges on the fixed 
surfaces and may induce undesirable roll moments, etc. On the other hand, in this 
type of airframe, the normal force on the control surface is in the direction of the 
desired maneuver, so that this feature aids the overall response of the guidance system. 




78 



3 Aerodynamic Forces and Coefficients 



• Wing control (e.g.. Sparrow III AIM-7F): 
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• Tail control (e.g., Phoenix AIM-54A): 




• Canard* (e.g.. Sidewinder AIM-9): 
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p 1 y 









Fig. 3.12. Airframe types. 



Similar generalities hold true for a canard airframe with control surfaces far forward 
of the eg. There are three typical types of airframes used in aerodynamically guided 
missiles (or munitions). Figure 3.12 illustrates these types. 

Each of these types has its own advantages and disadvantages, and the missile 
system designer must exploit these advantages and disadvantages for the given threat 
and operating environment. Also, the selected airframe must be able to deliver the 
required aerodynamic performance (see also Section 3.3.1). 

From Figure 3.12, we note that the typical Sparrow missile uses wing control 
with tail fins for stability and aerodynamics. The Phoenix (or Falcon ) missile uses 
tail control, via fins with fixed wings. The typical Sidewinder uses fixed tail wings, 
with movable nose fins (or canards). We note here that the type and size of airframe is 
strongly dependent on the guidance characteristics, motor size, and warhead size. It 
should be noted that situations do arise, however, in which some of these choices are 
dictated: for example, to use a motor existing in the inventory, or to use an existing 
airframe. Since typically these items were designed for some other threat, and with 
other system considerations, the system design problem focuses on obtaining the most 
out of these designs. 

The present discussion will be restricted to the tail-controlled configuration, which 
has no downwash interference from the control surfaces. If the autopilot pitch and 
yaw axes are each 45° from the planes of adjacent control surfaces, then all four 

*Canards are forward control surfaces, placed far forward of the center of gravity. An 
advantage of using canards is that the downwash from canards onto the main lifting body 
can, in certain configurations, nullify any attempt to control the missile in roll. 
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control surfaces are deflected equally by the pitch (or yaw) autopilot. However, in 
some applications it is preferable to put the autopilot axes in the planes of the control 
surfaces, so that only two surfaces are deflected by the pitch autopilot and two by 
the yaw autopilot. It is apparent that the tail-controlled airframe has a tail normal 
force opposite to the direction of the desired maneuver acceleration, which causes 
a small initial airframe acceleration in the wrong direction. Analytically, this effect 
manifests as a right-half-plane zero in the transfer function from the control surface 
deflection, <5, to a lateral (i.e., normal to the missile’s longitudinal axis) acceleration 
a n at the missile eg, thus tending to limit the speed of response of the guidance 
system. However, it is very desirable to have the roll autopilot command all four 
surfaces, so that the induced roll moments from “shadowing” of the control surfaces 
at high angles of attack will be minimized. In most cases, the control surfaces exert 
only aerodynamic forces, although it is possible to augment these at high altitudes by 
reaction jets embedded in the control surfaces. 

In Section 3.2 the various forces and moments acting on a missile were devel- 
oped. In this section, we will discuss the stability conditions in terms of the airframe 
dynamics. Specifically, in an aerodynamically maneuvering missile, the function of 
the control surfaces is to exert a moment so that the missile can develop an angle 
of attack (ACM), and thereby achieve lift from the body and wings (if any). (Note 
that in some applications, a wingless airframe relying only on body-lift is preferred.) 
As will be discussed further in the autopilot section, the pitching of the airframe 
causes the seeker to develop a spurious component of the measured line-of-sight 
(LOS) angular rate dX/dt , which results in a parasitic attitude loop* An impor- 
tant measure of the necessary pitching of the airframe is the alpha over gamma dot 
(a/y) time constant, r, of the linearized airframe response, defined by the following 
relation: 



where 



r — a/y — 



2 M 
pV m S 



dC L 

da 



9 C,„ 

HT 

/ 9 c,„ \ 9C L /dc„,\ 

V 9S ) dS V 9a ) _ 



M = mass of the missile, 

V m — missile speed, 

S — missile reference area, 
a — angle of attack, 
y — flight path angle, 
p — air (or atmospheric) density, 



*A parasitic attitude loop is defined as a control loop that interferes with the guidance 
stability, resulting in a larger miss distance. 
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(d = diameter ot the missile) 

Fig. 3.13. Aerodynamic missile stability conditions. 



and the aerodynamic derivatives are conventionally defined. Here, it is assumed that 
thrust and drag are negligible. The (ajy) term in the above equation is the angle of 
attack required to achieve a given flight path rate, and therefore represents the gimbal 
angle change in achieving a given missile maneuver. At high altitudes, it is necessary 
for the airframe to have a relatively large angle of attack in order to develop 1 g 
of lateral acceleration. The ratio a/y can be reduced, if it becomes large, by using 
movable wings near the eg, so that very little body pitching is required to develop 
lateral acceleration. However, a practical difficulty with this scheme is that the body 
eg moves appreciably as the booster (and sustainer, if used) motor bums. 

Based on the discussion of this section, the stability condition for an aerodynamic 
missile is illustrated in Figure 3.13 (see also Sections 3.1 and 3.3.2 for a discussion 
of the center of pressure ( cp )). 

The quantity M a is a measure of responsiveness of the airframe to pitch-moment 
changes with alpha. A small positive M a and slightly unstable airframe can be toler- 
ated, but a highly negative M a is preferable for stability of the autopilot and attitude 
loop. It appears that a zero or slightly negative M a (fairly stable airframe) is the best 
total compromise that the guidance designer could ask for. In words, the conditions 
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depicted in Figure 3.13 can be stated as follows: (a) If the center of pressure ( cp ) 
is ahead of the center of gravity (eg), the missile is said to be statically unstable; 
(b) if the cp is behind the eg, the missile is said to be statically stable; and (c) if 
the cp and eg coincide, then the missile is said to be neutrally stable. Note that the 
distance between the cp and eg is called the “static margin.” A summary of some of 
the important airframe characteristics is given below: 



G pr — pitch rate 9 m /surface angle 8 
= K pr (l + rs)/(l+b n s + b n s 2 ) 

= M s [ 1 + (1 /s)]/[s 2 + (b n /bn)s - M a ], 



1 2 Sd 

Ms — ~pV m ~j C ms \ cg , 

z 1 yy 



1 2 Sd 

M a =-pVl—C, 



Stable if M a < 0, 



/ cg Unstable if M a > 0, 

yy 



t = o‘/y\s.s 



Gi a = lateral accel. /(//surface angled 

= K] a ( 1 + nips + ci\2s~) / (1 + fens + bns 2 ) 



1 2 Sd 

K & = -pV m —Cu 



i XX 



Ms_ 

K S 



I XX C ws 
Iyy C h 



eg 



where 



Ms — Surface pitch effectiveness, 

K f , — Surface roll effectiveness, 
s.s — steady-state (subscript). 

Note that instabilities can be eliminated by means of a suitable feedback control 
system. Finally, it is noted that the majority of tactical missiles have fixed main lifting 
surfaces (often called wings) with their cp somewhere near the missile eg, and rear 
control surfaces. In subsonic missiles, it may be more efficient to use control surfaces 
as flaps immediately behind the wings, since the flaps control the circulation over the 
entire surface. In supersonic flow, the control surface cannot affect the flow ahead 
of itself, and therefore it is placed as far to the rear as possible in order to exert the 
maximum moment on the missile. Rear control surfaces often make for a convenient 
placement and/or arrangement of components [2], 

For those readers interested in aerodynamics, it should be mentioned that an 
experimental flexible-wing jet made its first flight in November 2002, from NASA’s 
Dryden Flight Research Center at Edwards AFB, California. During an hour-long 
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test, the modified F/A-18 took off, climbed to 30,000 ft., and flew over a test range 
northeast of Edwards AFB. The first flight followed a three-year period of modifica- 
tion and ground tests at the NASA facility. Officials from AF Research Laboratory, 
Boeing’s Phantom Works, and NASA Dryden collaborated on the research effort, 
called the Active Aeroelastic Wing ( AAW) program. The project intends to demon- 
strate improved aircraft roll control through aerodynamically induced wing twist on 
a full-scale manned supersonic aircraft. The new wing technology is important to the 
AF because it represents a new approach to designing wings that are more, structurally 
and aerodynamically from a control effectiveness standpoint. With AAW, leading and 
trailing edge control surfaces are deflected, which causes a change in the aerodynamic 
pressure distribution on the wing’s surface, causing it to warp or twist. The surfaces 
are deflected such that the wing twists into a shape that helps the wing perform better 
than it would if it did not twist at all. 

Moreover, AAW is applicable to a wide variety of future air vehicle concepts that 
are under study and not only to supersonic flight. While the technology was conceived 
during a supersonic fighter aircraft design study, aircraft that fly subsonically can also 
experience a high degree of wing deformation, and therefore could benefit from the 
AAW design approach. Since AAW exploits wing flexibility, it also is viewed as a first 
step toward future “morphing” wings that can sense their environment and adapt their 
shape to perform optimally in a wide range of flight conditions. 

The test-bed aircraft was modified with additional actuators, a split leading edge 
flap actuation system, and thinner skins on a portion of the upper wing surface that 
will allow the outer wing panels to twist up to 5°. With the first flight completed, 
the modified aircraft will undergo 30 to 40 flights over a three to four month period. 
NASA expects the second phase of the research flights with new control software to 
begin in mid to late 2003. 



Example 2. In this example we will apply the most important results of Sections 3.1 
and 3.2. Consider an air-to-air missile represented in Figure 3.14. 

The problem will be formulated as follows: 



Missile Dynamics 



dq 

dt 



— — {[— C(a, M, h, t) + T (t)]/m(t)} — qw 
dt 

= {— N[a, M, 8 z )/m(t)} + qu , 

dt / c u 

[M y (a, M, 8 Z , X cg -X cp ,t]-q l -Jp 






where 



u = longitudinal velocity component, 
C — axial force, 

M — Mach number = V m / y s , 

V m = missile speed, 

V s = speed of sound, 
a = angle of attack = tan _1 (w/u;), 



( 1 ) 

( 2 ) 

(3) 
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P 





Fig. 3.14. Representation of the missile in the local reference frame XYZ. 



h = flight altitude, 
t = time, 

m = mass of missile, 
w = normal velocity component, 
q = pitch rate, 

T = rocket motor thrust, 

N = normal force (yawing moment), 

S z = thrust deflection (or tail fin) angle, 

M — pitching moment, 

X cg , X cp = center of gravity and center of aerodynamic pressure, respectively, 

I y = missile’s moment of inertia about the 7-axis. 

The aerodynamic forces with respect to the local reference frame XYZ are as follows: 
Axial Force 



C — qC c (oi, M, h, t), 



(4) 
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C c = F fr (M){[(h - 6.096)/3.048] + F afc (a, M) + F b (M , t)[l - (A e /A h )]. (5) 
Normal Force 



N — qC n (a, M, S z ), 


(6) 


C n = F n f c (a, M)a - F tnf (M)8 z , 


(7) 



where 

F n f c = normal force coefficient, 

F tn f = trim normal force effectiveness, 
F b — base drag coefficient, 

F f r = friction drag coefficient, 

F a f c = axial force coefficient. 

Aerodynamic Couple 



M y — qL r C m (a , M, 8 Z , q , f), 

C Ml — C ma & “1“ Cm8z^z ”1” Cmqtf ? 

C ma = F n f c (u, M){[X cg (t) - X cp (a , M)]/L r }, 
C,„Sz = Ftpm(M) - F tn f(M){[X cg (t) - X cg b\/L r ], 
Cmq — -Fpdm (M)]500L r /V m , 

where 



( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 



q — dynamic force = S r q, 

S r — reference area, 
q = dynamic pressure = \p a (h) V^, 
p 0 = volumic air mass, 

A b = base area, 

A e — nozzle exit area, 

L r = reference length, 

X cg b = center of gravity location at burnout. 

The aerodynamic forces as given in the above equations are functions of aerodynamic 
coefficients, which depend on the angle of attack a and on the Mach number M. 
Values of these coefficients for various a and M are available in tabular form, which 
are obtained from wind tunnel experiments [2], [8]. Note also that numerical values 
for p„(h ), V s , T ( t ), and X cp (t) are tabulated along with the analytical expressions of 
the functions m(t), X cg {t), and I y (t). 
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3.3 System Design and Missile Mathematical Model 

3.3.1 System Design 

In general, a missile can be defined as an aerospace vehicle with varying guidance 
capabilities that is self-propelled through space for the purpose of inflicting damage 
on a designated target (see also Chapter 1). These vehicles are fabricated for air-to-air, 
surface-to-air, or surface-to-surface roles. They contain a propulsion system, warhead 
section, guidance systems, and control surfaces, although hypervelocity missiles do 
not use warheads or control surfaces. The guidance capabilities of the different mis- 
siles vary from self-guided to complete dependence on the launch equipment for 
guidance signals. 

Specifically, a guided missile is typically divided into four subsystems: (1) the 
airframe, (2) guidance, (3) motor (or propulsion), and (4) warhead. These subsystems 
will now be examined in more detail: 

1 . Airframe: The type and size of airframe is strongly dependent on guidance char- 
acteristics, motor size, and warhead size (see also Section 3.2.1). 

2. Guidance: The type of guidance that can be used is also dependent on the motor, 
warhead, and threat. More specifically, the type of guidance chosen is dependent 
on the overall weapon system in which the missile will be used, on the type of 
threat the missile will be used against, the characteristics of the threat target, and 
other factors. Guidance, as we have seen earlier, is the means by which a missile 
steers or is steered to a target. 

3. Motor: The motor characteristics are dependent on guidance requirements, the 
threat, and the airframe characteristics. 

4. Warhead: The warhead is dependent on the threat and type of guidance. 
Commonly, the procedure is to size the guidance requirements (e.g., accuracy, 
response time, range capability) from the threat, select an airframe that can 
deliver the required aerodynamic performance, size the motor based on threat 
and airframe considerations, and size the warhead from guidance and airframe 
considerations. 

Figure 3.15 illustrates the above missile characteristics for an aerodynamic air- 
to-air missile. 

In addition to the above considerations, there are other basic factors that affect the 
design of any weapon system. These are ( 1 ) the threat, (2) the operating environment, 
(3) cost, and (4) state of the art. Typically, the threat and operating environment are 
known or are given. Also, the state of the art is known. Therefore, the design effort 
of any missile centers on meeting the threat in the environment with the state of the 
art, at minimum cost. Consequently, three of these four factors are specified, with the 
fourth being either minimized (i.e., cost, state of the art) or maximized. 

In Section 3.2 the boost and sustain motor thrusting methods were discussed. 
However, in certain applications, missiles are designed to be of the boost-sustain 
type. Consequently, a major consideration in any missile design is the motor (or 
propulsion system) type selection. The important factors in selecting a motor type are 
(1) aerodynamic heating due to the incremental missile velocity, (2) aerodynamic 
drag, which decreases missile velocity, (3) maximum altitude at which the missile 
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Canards 




1 Radome 

2 Planar array active radar 
antenna 

3 Proximity fuze antenna 



8 Moving control fins 

9 Nozzle 

10 Rear detection antenna 

1 1 Hydraulic power unit 




must perform, and (4) maximum and minimum intercept ranges required. For the 
reader’s convenience, these missile motor types are summarized below: 

All-Boost: An all-boost motor typically will make the missile accelerate rapidly, 
causing high peak velocities. However, this causes high missile drag, high aero- 
dynamic heating, and short time of flight, for a given range. This motor is suitable 
for a rear hemisphere, tail chase encounter. 



Thrust 




All boost 



Time 
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All-Sustain: The all-sustain motor has low missile acceleration, resulting in lower 
aerodynamic drag and longer time of flight, for a given range. Since the motor 
burns for a long period of time, the motor can be used to overcome gravity 
in a look-up engagement, and to provide sufficient velocity for maneuvering 
at high altitude. This motor is suitable for head-on engagements, or in look-up 
engagements at high altitudes. 



Thrust 



V 

All sustain 



Time 



Boost-Sustain: The boost-sustain motor represents an attempt to combine the 
best features of the all-boost and the all-sustain designs. 



Thrust 




The warhead is typically an input to the system design. Since almost all warheads 
for missiles are GFE (government furnished equipment), we see that the role of the 
system designer is to design a system to deliver a given warhead to a given point in 
space, within a given accuracy, and to fuze it at the appropriate time. 

The threat typically describes what the target can do, in terms of performance, 
which translates into what the missile has to do, in terms of performance. These 
important threat factors are: 

1. The rate of closure: This is the combination of interceptor velocity and target 
velocity. 

2. Engagement altitude: This is the altitude regime over which the target can be 
expected. 

3. The engagement range: This is the limit over which the missile can be launched. 

All of these factors reflect missile requirements or constraints. 
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In addition to solid and/or liquid fuel missile engines (or motor), there is another 
class, namely the ramjet. The ramjet is the simplest (in terms of moving parts) of all 
air-breathing engines, in which the compression is produced by the ram, or forward 
velocity pressure in the intake; combustion downstream of the intake thus gives the 
high-velocity jet. Thus, the ram effect can be defined as the pressure increase obtained 
in an air-breathing engine intake by virtue of its forward motion. 

Ramjet propulsion is a key ingredient of new antiradar missile technology. With 
a Mach 3 or more capability, ram rocket propulsion will have increased range and 
reduced time to target. Research is also being conducted into the so-called scramjet 
engines. These engines would have the power to hurtle through the air at Mach 
6 or more. The engines initially will be installed in cruise missiles, making them 
hypersonic, meaning that they move faster than five times the speed of sound, and 
more deadly. Cruise missiles, which are launched from aircraft, now travel slower 
than the speed of sound. The scramjet engine should enable cruise missiles to strike 
targets such as transportable Scud missile launchers or fixed targets that must be taken 
out quickly. It also could be used to hit deeply buried targets. A 1996 study conducted 
at Maxwell Air Force Base in Alabama showed that hypersonic cruise missiles could 
have a range of more than 1,000 miles, allowing the aircraft from which they are 
launched to stay out of harm’s way. 

Unlike conventional jet engines, the scramjet does not require rotating fans and 
compressors. Instead, it relies on the forward motion of the vehicle to compress the 
air. Once inside the engine, the air mixes with the injected fuel and is burned. The hot 
gas exits the rear of the engine and provides thrust by pushing against a nozzle-like 
surface. The biggest hurdle for researchers is dealing with the high speed of the wind 
as it enters the engine. The main issue is to get the fuel and the air mixed well enough 
so that fuel and air can burn in the very limited time available. 

A European development of a new beyond-visual-range air-to-air missile 
(BVRAAM) would employ a Mach 3 ram-rocket propulsion system for increased 
range and reduced time to target. The ram-rocket motor would feature four inlets in 
the center of the missile body and high boron content in the sustainer propellant for 
high specific impulse with low volume. After being boosted to the required operating 
speed, the air-breathing ramjet sustainer would take over for the rest of the flight, mix- 
ing fuel-rich gas from a boron gas generator. The Raytheon Company also is offering 
a next-generation version of the AIM- 120 AMRAAM, dubbed the future medium- 
range air-to-air missile ( FMRAAM ). Raytheon’s FMRAAM design employs a liquid 
fuel ramjet developed by Aerospatiale Missiles. Still another European concept is to 
employ a solid fuel, variable flow ducted ramjet developed by Germany’s DASA sub- 
sidiary Bayern-Chemie. The Russian Kh-31, which has an active or passive RF seeker 
for antiship or antiradiation missions; is one of the few operational ramjet missiles; 
it flies at Mach 2.7 while sea-skimming. The Defense Advanced Research Projects 
Agency’s ( DARPA ) Affordable Rapid Response Missile Demonstration ( ARRMD ) is 
examining two different concepts: one from the Air Force’s HyTech program and the 
other from the Office of Naval Research’s Hypersonic Weapon Technology Program. 
The Aerojet Corporation, which builds the dual combustion ramjet ( DCR ), proposes 
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Fig. 3.16. Missile development phases. 



a central subsonic ramjet combustor that feeds fuel-rich exhaust into the surrounding 
supersonic scramjet stream, where combustion is completed. Before we leave this 
topic, it should be pointed out that there is another type of propellant, namely, the 
hypergolic* propellant. This is a rocket propellant that ignites spontaneously upon 
contact with an oxidizer. 

In the development of missile systems, a broad spectrum of engineering dis- 
ciplines is employed with primary emphasis placed on the guidance and control 
subsystems. A comprehensive systems-oriented approach is applied throughout the 
system development. Specifically, there are four basic phases in the development of 
a missile. These are: 

1. Concept formulation and/or definition; 

2. Requirements; 

3. Design; 

4. Evaluation. 

In diagram form, these four phases are shown in Figure 3.16. 

Consequently, the development cycle for a missile system commences with the 
concept formulation, where one or more guidance methods are postulated and exam- 
ined for feasibility and compatibility with the total system objectives and constraints. 
Surviving candidates are compared quantitatively and a baseline concept adopted. 
Specific subsystem and component requirements are generated via extensive tradeoff 
and parametric studies. In particular, such factors as missile capability (e.g., acceler- 
ation and response time), sensor function (e.g., tracking, illumination), accuracy (i.e, 

* Hypergolic is a coined word, the element golic being obtained from a German code word 
Gols, used to refer to a series of rocket propellants containing methylaniline, organic amine, 
and certain other compounds. 
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SNR and waveforms), and weapons control (e.g., fire control logic, guidance software) 
are established by means of both analytical and simulation techniques. After iteration 
of the concept/requirements phases, the analytical design is initiated. The guidance 
law is refined and detailed, a missile autopilot and the accompanying control actuator 
are designed, and an onboard sensor tracking and stabilization system devised. This 
design phase entails the extensive use of feedback control theory and the analysis 
of nonlinear, nonstationary dynamic systems subjected to deterministic and random 
inputs. Determination of the sources of error and their propagation through the system 
are of fundamental importance in setting design specifications and achieving a well- 
balanced design. Finally, overall evaluation is conducted using numerous simulation 
techniques and test facilities to verify the design, predict performance, define zones 
of effectiveness, and analyze post-flight results. 

As stated above, classical servomechanism theory has been used extensively to 
design both hydraulic and electric seeker servos that are compatible with require- 
ments for gyro-stabilization and fast response. Pitch, yaw, and roll autopilots have 
been designed to meet problems of Mach variation, altitude variation, induced roll 
moments, instrumentation lags, body-bending modes, guidance response, and stabil- 
ity. Although classical theory presently is most applicable to autopilots and airframe 
stability, research efforts are continually made to apply modern control and estima- 
tion theory to conventional as well as adaptive autopilot design. Modern advanced 
guidance and control systems having superior performance have been designed with 
on-line Kalman filter estimation for filtering noisy radar data and with optimal con- 
trol gains expressed in closed form. Synthesis of sampled data homing and command 
guidance systems are being used extensively today. For example, a vital point in 
the development of advanced homing guidance and control systems is to optimize 
the performance of the missile under design for various intercept situations and tar- 
get maneuvers. Furthermore, trends in operational requirements indicate that future 
air-to-air or air-to ground missiles will have to have a high probability of kill under 
total sphere launch engagement conditions and a launch and leave capability (such as 
the AGM-154 Joint Standoff Weapon (JSOW)) when employed against a wide vari- 
ety of highly maneuverable intelligent targets. In order to satisfy these requirements, 
future air-to-air missiles will require complex guidance algorithms. Additionally, 
in order to implement these complex guidance algorithms, more information about 
the missile and target dynamic states will have to be accurately measured or esti- 
mated on board the missile. The very nature of this problem lends itself to the use of 
modern control theory to derive the advanced guidance laws and modern estimation 
theory to develop techniques to process the available information and estimate the 
unavailable information. The key to a successful problem formulation is a translation 
of the mission requirements into a mathematical performance index (P.I.). No matter 
what theoretical techniques are used to develop the optimal control strategy, they will 
always be based on the minimization (or maximization) of some performance index. 
Hence, the optimal control strategy will be no better than our selection of the critical 
performance drivers and our translation of these into concise mathematical terms. In 
addition to the performance index, two other formulations will impact the optimal 
control strategy, namely, the mathematical model of the system and the additional 
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equality and inequality constraints to be placed on the system. In general, a more 
detailed system model results in a more accurate control strategy, but this is achieved 
at the expense of additional complexity in both the derivation and resulting algorithms. 
The selection of appropriate equality/inequality constraints can be based upon either 
actual system parameters or trajectory properties that we want the optimal solution to 
possess. Some of the modern control techniques that have been investigated and/or 
applied are (1) reachable set theory, (2) singular perturbation theory, (3) differential 
game theory, and (4) adaptive control theory. 

The performance index (also known as cost functional) study is a fundamental but 
extremely complex problem. To be sure, the mission objectives influence the choice 
of the parameters/states that are included in the performance index. However, there 
are many other factors that must be considered in the construction of the performance 
index. These additional factors are due to the interrelationships of the steps involved in 
the modern control problem formulation. Specifically, for every different performance 
index or cost functional there is a different optimal guidance law. In essence, the 
measure of performance for any guidance law will be the ability of the missile to hit 
the target (e.g., minimum terminal miss distance) in various engagement scenarios. 
Additional measures of performance involve (1) launch envelope (full 360° launch 
aspect angle, minimum inner launch boundary), (2) fuel considerations, (3) time 
considerations, and (4) maneuver capability. 

3.3.2 The Missile Mathematical Model 

Guided tactical missiles are sometimes referred to according to their airspeed relative 
to the speed of sound and their type of propulsion system. Generally, the highest rate 
of airspeed that can be reached safely and still ensure correct operation is considered 
as that missile’s classification. In essence, the general means of classification of a mis- 
sile’s airspeed is related to the speed of sound (or Mach 7), which varies with respect 
to the ambient temperature. Commonly, there are four groups that are considered in 
classifying a particular missile. These are: 

1. Subsonic: Airspeeds less than Mach 1. 

2. Sonic: Airspeeds equal to Mach 1. 

3. Supersonic: Airspeeds ranging between Mach 1 and Mach 5. 

4. Hypersonic: Airspeeds exceeding Mach 5. 

Practically all AIM and SAM missiles can be placed in the supersonic classification, 
since modern military aircraft are capable of attaining Mach 1 speed. 

A commonly used missile mathematical model in the analysis and design of 
surface-to-air and air-to-air weapon systems is the skid-to-turn ( STD missile, in which 
both the pitch and yaw plane systems have identical response behavior. As discussed 
in Section 3.1, the main reason for using the skid-to-turn design is that the inertial cross 
coupling between roll, pitch, and yaw is negligible. For this reason, we will assume that 
our missile mathematical model is of the skid-to-turn type. However, as noted earlier, 
in this technique both aerodynamics (the aerodynamics cannot be described in closed 
form, but are available in look-up table form) and rigid-body dynamics are highly 




92 



3 Aerodynamic Forces and Coefficients 



Table 3.2. Comparison of Weapon System Characteristics 



Feature 


Sparrow-type weapon 


Bank-to-turn weapon 


Guidance Mode 


Skid-to-Turn 


Bank-to-Turn 


Control Surfaces 


Wing Control 


Tail Control 


Autopilot Sensors 


Accelerometers/Rate Gyros 


Accelerometers/Rate Gyros 


Maximum Acceleration 


45g’s (32 g’s per axis) 


100 g’s in pitch, 10 g’s in yaw 


Guidance Delay 


0.75 seconds 


0.40 seconds 


Launch Speed 


0. 5-2.0 Mach 


0.5-2.0 Mach 


Speed Range 


0. 5-3.0 Mach 


0.5-4.0 Mach 


Maximum Roll Rate 


Not Applicable 


±600°/sec 



nonlinear (see Chapter 2). The highly nonlinear aerodynamic effects occur at high 
Mach numbers. Also, the skid-to-turn missile will experience difficulty when attack- 
ing high-g targets. The Sparrow weapon is a skid-to-turn type missile. At this point it 
should be mentioned that another missile configuration is the bank-to-turn ( BTT ) type. 

This configuration is suitable for highly maneuverable ramjet missiles. However, 
the asymmetrical airframe of this missile design requires rolling the missile to main- 
tain target motion in the missile pitch plane. Another drawback of this design is 
that the bank angle causes a coupling between the pitch and yaw channel dynamics, 
which can vary considerably even for a short period of time. Table 3.2 summarizes 
the skid-to-turn and bank-to-turn missile characteristics. 

Specifically, the pitch/yaw plane rotational responses behave like a spring- 
mass damper system. Mathematically, this system response can be expressed in the 
form 



+2£tt> + 0 J 2 y — u> 2 u[t). (3.54) 

Equation (3.54) can also be written in the usual frequency domain as follows: 

y(s)/u(s) =ar/(s 2 + 2 t;cos + or), (3.55) 



where 



y(s) = output, 
u(s ) = input, 

£ = damping ratio (dimensionless), 
u> — frequency (rad/sec), 
s = Laplace operator (rad/sec). 

The above continuous system can be constructed in a simulation as a feedback net- 
work that represents a load factor command system in both pitch/yaw planes (see also 
Section 3.5). Figure 3.17 represents a typical pitch/yaw network. 
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Proportional + integral path M a q 




Fig. 3.17. Pitch-yaw feedback network. 



Let us now digress for a moment and return to (3.1). In (3.1) the drag was given 
simply as a function of the coefficient of drag, dynamic pressure, and reference area. 
However, drag is a nonlinear function of velocity. For the purposes of control design, 
drag can be modeled by the parabolic drag form 

D = qSC D + (KL 2 )/qS, (3.56) 

where K is defined as in (3.15b). In the present analysis, we will consider lift as the 
control. Lift is chosen subject to the constraint 

L<Wg m ( v), (3.57) 

where W is the weight and g m (v) represents the load factor limit, which may arise 
due to a structural limit, control surface actuator, or autopilot stability considerations. 
In general, lift is a function of missile speed. From the above discussion, the load 
factor is simply expressed by the equation 

g m (v) = r 1 = L/W = jpV 2 SC L /W. (3.58) 

The dynamics for the angle of attack ( AOA ), a, as well as dot/dt, load factor n z , and 
pitch rate, are commonly modeled after the short-period approximations of longitu- 
dinal motion. The short-period approximation for the angle of attack is given by the 
following transfer function: 

JHS L= + 1 ) ( 3 . 59 .) 

acmd(s) S z + 2^ cos + co l 

where 

T a = AOA time constant (sec), 

£ = short-period damping ratio 
(dimensionless), 

co = short-period frequency (rad/sec), 
s — Laplace transform operator (rad/sec). 
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M Alt Thrust 




Fig. 3.18. Load factor command system. 



This transfer function is expressed in the frequency domain and represents the lower- 
order equivalent system of the full closed-loop system. We note that the load factor 
and angle-of-attack transfer functions are identical in form. Specifically, the dynamics 
for the load factor in the pitch plane, n z , can be modeled by the following transfer 
function: 



n z (s) _ w 2 (T a s+ 1) 
WzcmdO) s 2 + 2^ cos + co 2 ’ 



(3.59b) 



where 



n zc = commanded load factor in the pitch plane, 

T a = AOA time constant (sec), 

and &>, C, and 5 are defined as in (3.59a). The parameters co, £, and T a can be found 
by linear analysis of the entire closed-loop system. The above transfer function is 
valid, provided that the load factor being modeled is located at the center of pres- 
sure ( cp ), that is, the point ahead of the center of gravity (eg) where the effect 
of pitch acceleration and horizontal tail force cancel. Moreover, load factor mea- 
sured at the center of pressure will reflect forces mostly due to angle of attack, 
which is why it has the same transfer function form. This assumption eliminates 
having to deal with a pair of complex second-order zeros in the numerator for an n z 
accelerometer located away from the cp. Figure 3.18 illustrates the load factor 
command system. 

We have seen that in terms of the continuous system (3.55), a feedback network 
for the purposes of simulation that represents the load factor command system in both 
the pitch/yaw planes was shown in Figure 3.17. The load factor control system can 
be reduced to have an identical response to that of (3.55) with y = n and u — n cm d- 
Specifically, this is done by properly computing the inner-loop parameters. This can 
be done by linearizing the airframe normal forces in the pitch and yaw planes to obtain 
the slopes n za and n v p. Since both the pitch and yaw use the same aerodynamic data, 
let n y fj — —n za or n ZC[ — nx, where X is either a or (the sideslip angle). Then compute 
the slope at each frame. The linearization of the feedback network (see Figure 3.17), 
reduces to that shown in Figure 3.19. 
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Fig. 3.19. Linearized load factor control system. 



In Figure 3.19, Kj is an integrator gain and G is a gain whose function is to force 
the response of the reduced-order system to behave like the original system. Note that 
this system is identical to that of a digital aircraft yaw axis system. 

Using the methods of classical feedback control theory, the inner-loop parameters 
are as follows: 



1 2tco~i -K[ 1 / ~ Z ZZ 

Y= 2 +2 V 4 ^ w - 4 ^i + *U’ (3.60a) 

1 2,ta)-\-K\ 1 / Z Z 9 

jr= 2 + (3 ' 60b) 



orT { T 2 

ni 



(3.60c) 



K i = 2C,cd. 



(3.60d) 



Next, we need the expressions for £ and co. In order to solve for these two parameters, 
we first write the second-order differential equation for the angle-of-attack response 
at launch with a forcing function of order 0. Thus, 



d 2 a 
dt 2 



+ 2£cu 




+ tt> 2 a: = 0. 



(3.61) 



It is well known that the response in a for this system will be a damped oscillatory 
motion that decays to zero. Equation (3.61) has the solution a — re 1 '. Substituting 
results in 



A . 2 T 2 £ cvX T co~ — 0. 

Solving for the eigenvalue X yields two solutions as follows: 



(3.62) 



2c i,2 = — j co- 



Cd. 63) 
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where j is a complex number. Now let 

r\ — a — jb for Ai, 
ri — a + jb for A 2 . 

Substitute both solutions into (3.61) and add the equations to obtain the result for 
a(t). Thus, 



a(t) = e i ’ mt cos cotJl — ^ + bsincotJl—t; 2 ^ . (3.64) 



Differentiating a(t), we obtain 



da(t ) 
clt 



= e 



bcoyj 1 — £ 2 — a^co^j cos cotj 1 — tj 2 



— — C 2 + btjco \ sin cot^jl — £ 2 



(3.65) 



Given the initial conditions ao and dao/dt, a and b can be solved for by setting t — 0. 
Thus, 

a = ao, 

do + C<w«o 

b= 

<»y/ 1-C 2 

Therefore, the time responses for a(i) and da(t)/dt become 



y(t)—aoe ( coscot J 1 — £ 2 H ^ sin wt J \ — tj 2 

\ VW 2 v 



(3.66) 



^ ^ = aoe (cos cuiV 1 — t, 2 ^ sin cot J 1 — 

v Vi-c 2 ) 

W -Scot , 



t/f 



— ap— — e f sincof^/l — £ 2 . 

v i - r 2 



(3.67) 



Next, if a tip-off behavior exists at launch, we can solve for the peaks by setting 
da(t)/dt — 0: 



0 = — crp e sin < vtj 1 — £ 2 . 

yi-c 2 v 

We note that peaks are periodic according to the sine wave function. Therefore, time 
to peaks are given by 




«=0, l,2,...,oo 



(3.68) 
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Fig. 3.20. Roll rate response. 



The time of the first overshoot (i.e., n = 1) is the largest overshoot. Next, we solve for 
the peak overshoot at f pea k for n— 1 : 



t^peak — Otoe 




(3.69) 



Finally, we can solve for £ from (3.69), giving 



f = 









^ 2 + [/w (_^)] 2 



(3.70) 



and for a> from (3.68), yielding 



a> = 



n 

t ] peak \/ 1 



(3.71) 



A few words about the roll axis model are in order. The missile roll dynamics for stabil- 
ity axis roll rate, P, can be modeled after the roll approximation of lateral/directional 
motion. This approximation ignores the coupling effects in the rotary cross terms and 
in sideslip angle fi. The roll approximation for roll rate is given by the following 
simple transfer function: 



TstabOO 1 

T’stabcnnj ( s ) TrS + 1 



(3.72) 



where Tr is the augmented roll mode time constant with units of seconds. It should 
be noted that this represents the lower-order equivalent system of the full closed-loop 
airframe. The filter time constant may be found from linear analysis of the entire 
closed-loop system or by computing it with the augmented stability derivatives. From 
classical control theory, we note that the roll rate response for a step input command 
is as illustrated in Figure 3.20. 
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In Figure 3.20 we note that the time constant Tr corresponds to the time at which 
the roll rate response reaches 63% of its final value. 

Example 3. Assume for simplicity that the missile’s motion is constrained in the 
vertical plane. Furthermore, we will assume that the missile can be modeled as a 
point mass. Therefore, from the missile’s balanced forces shown in the diagram 
below, we can write the equations of motion as follows (see also Example 1 in 
Chapter 2): 




Equations of Motion 



dV 

dt 

dy 

dt 



( \/m)[T cos a — D] — g sin y, 
(\/mV)[L + T sin a] — (g/V) cos y, 



dx 

— = V cos y, 
dt 

dh 

— — V sin y. 
dt 



Aerodynamic Derivative Coefficients 



L = \pV 2 SC L , 

D = \pV 2 SC D , 
Cl = CLalfit — a 0 ), 
Cd = Coo + kCi, 
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where 

g — acceleration of gravity, 
h = altitude, 

k — induced drag coefficient, 
m — mass, 

D — drag, 

L = lift, 

M — Mach number, 

S = reference area, 

T — thrust, 

V — velocity, 

Cd, Cl — drag and lift coefficients, 
respectively, 

C do — zero-lift drag coefficient, 

Cl„ = SCi/da. 

The aerodynamic derivative coefficients Ch a , Cd o, and k are functions of the Mach 
number as follows: 



C Lot = Clc(M), 

Cdo — Cdo(M), 

M = M(V, h), 
k = k(M), 

P = p(h). 

Moreover, the missile mass and the thrust are functions of time; that is, 

/ n—m(t) and T — T(t). 

Note that the AOA can be treated as a control variable (if used in connection with an 
optimization case) that satisfies the inequality constraint 

&min — & — ^ max ■ 



3.4 The Missile Guidance System Model 

This section briefly describes the basic subsystems that form a missile’s guidance 
system. Guidance is the means by which a missile steers, or is steered, to a target. 
A guided missile is guided according to a certain guidance law. In this chapter we 
consider homing guidance systems. A meaningful comparison of homing guidance 
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Missile 




Fig. 3.22. Major subsystems of a missile guidance system. 



systems for missiles requires realistic models for the missile and its target engagement 
geometry model in order for terminal miss distance to be accurately evaluated. This 
model should include the important dynamics and system nonlinearities that influ- 
ence performance, and yet be representative of missiles in general. In the simplest 
form, the principal elements that make up a missile guidance system are illustrated in 
Figure 3.22. 

For a missile, the inputs are target location and missile-to-target separation. The 
desired output is that the missile have the same location as the target. The missile does 
this by using a certain guidance system and flying according to a certain guidance 
law. As stated in Chapter 1 , the type of guidance system chosen is dependent on the 
overall weapon system in which the missile can be used, on the type of threat the 
missile will be used against, and the characteristics of the threat among other factors. 
The various subsystems indicated in Figure 3.22 will be discussed in the subsequent 
sections. It should be noted at the outset that the model developed herein assumes that 
the target and missile motions are constrained to a plane. Consequently, development 
of the missile and guidance models is limited to a single channel. 

A more detailed block diagram for a controlled missile guidance model than the 
one illustrated in Figure 3.22, which includes the equations of motion and aerody- 
namics, is given in Figure 3.23. Note that this model is for a roll-stabilized missile. 

Listed below are the three main problems that the guidance system designer must 
face in the design of a guidance system. 

General Problems of Guidance System Design 

1. Help to maximize the single-shot kill probability ( SSKP ) by minimizing the miss 
distance. 
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Fig. 3.23. A typical roll-stabilized missile guidance/kinematic loop. 



Sources of Miss Distance 

( 1 ) Initial heading error. 

(2) Acceleration bias. 

(3) Gyro drifts (if gyros are used in seeker stabilization). 

(4) Glint (scintillation noise). 

(5) Receiver noise. 

(6) Fading noise (except for monopulse). 

(7) Angle noise (due to varying refraction with frequency diversity). 

Noise components of miss depend on guidance-system response and thus a/y 

and refraction slope R. 

2. Preserve stability of the parasitic attitude loop (to be discussed in Section 3.5). 

(In maneuvers, missile pitching affects seeker boresight-error measurement.) 

3. Filtering. 

(1) Limit power consumption and saturation of the actuators. 

(2) Prevent noise from excessive hitting of dynamic-range limits, such as auto- 
pilot g-limits. 

A more detailed discussion of the above problems is in order. First, in order to maxi- 
mize the single-shot kill probability, the main problem of the guidance designer is to 
minimize the miss distance enumerated above under problem 1 . The seven sources 
of miss distance listed above are statistical in nature. Both the alpha over gamma 
dot of the airframe and the statistically varying radome -refraction slope R affect the 
speed of response of the guidance system and thereby affect the components of miss 
due to noise. Evaluation and partial optimization of total rss (root-sum-squares) miss 
distance can be performed rapidly and efficiently with a digital computer program. 
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The second major problem of the guidance designer is to preserve the stability 
of the parasitic attitude loop (for a discussion see Section 3.5). The third problem 
states that some filtering must be used to limit noise perturbations of the seeker 
and actuators, so that power consumption, saturation, g-limiting, etc., will not be 
excessive. A successful guidance design requires a compromise that meets all three 
major problems listed above. 

3.4.1 The Missile Seeker Subsystem 

In this subsection we discuss missile seekers and the role they play in guidance 
systems. The discussion herein is by no means exhaustive, and is intended to stimu- 
late further research. Missile seeker and radome error analysis has been carried out 
extensively. Basically, the main function of the missile seeker (also known as homing 
eye ) subsystem is to: 

(1) Provide the measurements of target motion required to mechanize the guidance 
law. 

(2) Track the target with the antenna or other energy-receiving device (e.g., radar, 
infrared (IR), laser, or optical). We note here that the antenna refers to any type 
of energy-collecting device. 

(3) Track the target continuously after acquisition. 

(4) Measure the LOS (line-of-sight) angular rate dk/dt. 

(5) Stabilize the seeker against a missile pitching rate clO m /dt (also, yawing rate) 
that may be much larger than the LOS rate dk/dt to be measured. 

(6) Measure the closing velocity V c (note that this is possible with some radars but 
difficult with infrared seekers). 

The typical classical seeker hardware consists of two or three gimbals on which 
are mounted gyroscopes (either conventional spinning mass or laser gyros) and an 
antenna. Most seekers of the gimballing variety have two gimbal axes, namely, yaw 
(or azimuth) and pitch (or elevation), which are orthogonal to the longitudinal axis, 
relying on the missile roll autopilot for roll stabilization. Space stabilization about the 
instrument axes is necessary, although a slow roll rate about the LOS itself is tolerable. 
Numerous gimbal configurations have been used in the past. Occasionally, a special 
application may justify a roll gimbal to roll-stabilize the whole seeker. In an active 
radio frequency ( RF ) seeker or passive IR seeker, two gimbals are commonly used, 
namely, one for azimuth and one for elevation measurement. However, only one gim- 
bal and its associated dynamics are required for a planar-motion model. Associated 
with each gimbal is a servomechanism, whose function is to adjust its angular orien- 
tation in response to the tracking error signal measured by, say, a radar receiver. (Note 
that here we will assume that the seeker consists of a radar receiver, unless otherwise 
specified.) In addition to the gimbaled systems, there are also body-mounted antenna 
systems (i.e., in a strapdown configuration), which do not use moveable gimbals in 
order to position the antenna. Moreover, these systems use either a fixed antenna 
position relative to the missile or electronic beam steering by means of a phased- 
array radar antenna. It is important to emphasize here that the use of electronic beam 
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steering is in many ways analogous to the gimbaled system as far as the guidance 
system operation is concerned. The fundamental measurement obtained from a 
homing sensor receiver is assumed to be the indicated angular position of the target 
relative to the antenna centerline or boresight. Figure 3.24 illustrates a typical missile 
radar seeker. 

The line of sight (LOS), X, is defined as the angle between a line from the center 
of the seeker antenna to the target, and some arbitrary nonrotating (e.g., inertial) 
reference line. Commonly, it is convenient to select this reference equal to the LOS 
position at the beginning of the homing guidance phase. Consequently, X(t) at time t is 
the total change in the angular position of the LOS relative to the initial LOS. Referring 
to Figure 3.24, we note that the angular position of the missile body centerline, 6 m , 
is measured relative to the initial LOS. (Note that () m is identified as the pitch angle.) 
Furthermore, the angular position of the antenna centerline is defined by the gimbal 
angle, On . Therefore, the LOS angle X is given by 

X — 0 m T 0j] T e, (3.73a) 

where e is the true boresight error, that is, the error between antenna center line and 
line of sight to the target. Alternatively, by writing (3.73a) as 



e — X 0f]i Ofi , 



(3.73b) 
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we obtain the expression for the true boresight error. It is important to note that the 
boresight error is a function of both the missile attitude relative to inertial space and the 
angular position of the antenna relative to the missile centerline. Since X or dX/dt is 
the desired measurement for guidance purposes, it is necessary to remove the missile 
motion from the LOS measurement data. One requirement on the seeker subsystem is 
to keep the antenna pointed at the target, so that the error e is always much smaller than 
the beam width of the received energy. Furthermore, in the region of small e, the seeker 
receiver measurement of indicated boresight error is nearly linear. However, if s cannot 
be considered small relative to the antenna beamwidth, the receiver boresight error 
processor operation may become nonlinear. Specifically, if s is allowed to approach 
the half beamwidth of the antenna, the receiver detection circuitry will at some point 
lose lock, and all guidance information will be lost. Therefore, the seeker must track 
the tracker sufficiently closely so that large boresight errors do not occur. Otherwise, 
the nonlinearity of the boresight error position should be considered as an important 
system nonlinearity. The actual form of the boresight error processor nonlinearity is 
strongly dependent on the specific beamwidth, processing scheme (i.e., monopulse 
radar, C W radar, etc.), and detector characteristics of individual systems. In the present 
discussion, it will be assumed that the beamwidth and tracking response of the seeker 
are adequate to keep the boresight error processor in its linear region. 

The radome forms the nosepiece of a missile and covers the R F head assembly 
of the target seeker. More specifically, the radome forms an important part of the 
external contour of the missile and becomes a vital link in the electromagnetic path 
of RF energy reflected from the target to the missile antenna. The same reasoning 
applies to an infrared seeker. In either case, aerodynamic, structural, and electrical 
requirements must be adjusted in order to produce optimum performance. 

Because of the presence of lags in the seeker tracking loop and radome refraction 
effects, the seeker will not point directly at the actual target. Instead, the seeker will 
point to the apparent target. The reader should be cautioned that radar reflectivity 
from a target is affected by the frequency of the radar. An aircraft design that is 
invisible to high-frequency fire-control radar may be plainly obvious to low-frequency 
search radar. The half-wave-length phenomena can be a factor. Parts of an aircraft 
that equal one-half of radar’s wavelength create a resonance that greatly increases 
radar reflectivity. The aberration (or refraction) angle error is the result of nonlinear 
distortions in the received energy as it passes through the protective covering (e.g., 
radome in the case of a radar homing sensor) over the antenna. This distortion produces 
a false boresight error signal, s' which is interpreted as an error in the angular position 
and motion of the target by the guidance system. Referring to Figure 3.25, we note 
that the indicated boresight error s' in the presence of radome aberration or refraction 
angle error 0 r is given by the expression [3], [5], [12] 

s' = X + 0 r — 9 m — Oh , (3.73c) 

where the radome aberration angle error 0 r is in general not constant, but is a function 
of the gimbal angle Of , ; that is, the radome angle error is a nonlinear function of the 
gimbal angle 0/,. Mathematically, this can be expressed by the relation 0 r = f (Oh)- 
It should be noted here that the size of the measurement error, that is, the angle 0 r , 
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depends on the orientation of the antenna with respect to the radome (antenna cover), 
which is fixed to the missile airframe. 

Furthermore, the dependency of 9 r on (),„ couples body motion into the boresight 
error signal, thus forming what is commonly called the parasitic attitude loop (for 
a discussion of the parasitic attitude loop see Section 3.5). This loop can drastically 
alter missile response characteristics and in turn increase the miss distance. This is 
particularly true at high altitudes, where the missile body motion tends to be greatest. 
The aberration (or refraction) angle error is normally a nonlinear function of the 
following factors: (a) the angle between the missile center line and the line of sight 
to the target (also known as the look angle, which is defined as (X — 6 m )); (b) the 
radome thickness distribution; (c) material; (d) radome shape; (e) manufacturing 
tolerances; (f) temperature; and (g) erosion of the surface during flight. In addition, the 
nonlinearity arises from such optical and electrical properties as frequency, standing 
waves inside the radome, and polarization of the received signal. Consequently, an 
accurate model may require a nonlinear, time-varying statistical characterization of 
the radome. From the above discussion, it can be said that the radome error magnitude 
can neither be precisely measured nor predicted. However, since these characteristics 
tend to vary over rather wide limits depending on the particular application and missile 
configuration, a constant refraction error slope model is used to capture the important 
body-coupling effect. 

From the above discussion, we note that radome aberration error is one of the 
errors contributing to the overall miss distance of a radar-guided homing missile. 
Figure 3.26 illustrates the aberration angle as a function of look angle [3], 

The derivation of the radome model can be obtained as follows. Let X be the true 
LOS and X m the measured LOS. Then, 



X m — X T 0 r , 



(3.74a) 
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Fig. 3.26. Aberration angle error as a function of look angle. 



where 0 r = /($/,). Taking the derivative of (3.74a), we have 

d hjji d X. dQy 

dt dt dt 

Now, since 6 r — 

( 3 0 r \ 

where (dOy/dO/,) is the radome error slope; that is, 

3 e r 

R=—. 

3 Oh 

Also, the following relation holds: 

d0 r _ / 3 0 r \ / dd h \ d6y _/ d0 r \ fdd h \ 
~di~~\ Wh) ydi ) ^ ~dt~~\deij \~dt~ ) ' 

The gimbal angle, Oh, can be obtained from Figure 3.23 as follows: 

Oji = ( X 0 W ) s , 

or 

0h=X — ( 0 m + e), 

where (X — (),„ ) is the look angle, and 

dOh dX d0 m ds 

dt dt dt dt 



(3.74b) 



(3.75) 



(3.76) 



(3.77) 
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Therefore, from (3.74a) and (3.74b) we have 



dX m 

dt 



dX d0 r dX ( 3 9 r 
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(3.78) 



Note that in (3.78) it has been assumed that the bias angle 0o is zero. The (1 + R) term 
in (3.78) is therefore the radome model; it shows how the actual LOS is perturbed. 
Moreover, this equation shows that the measured LOS, X m , is corrupted by the radome 
error slope. 

A linear model for the general aberration angle can be obtained using a Taylor 
series approximation of the form 



0,. = 0O + (X-dm)R, (3.79) 

where 0 q is a bias angle and R is the radome error slope. Substituting (3.79) into 
(3.73c) and solving for s' yields [5], [ 12] 

e' = (1 + R)(X - 0,„) + 00-0/,. (3.80) 

The radome error slope, R, which varies from radome to radome, is the main parameter 
in the homing loop. Mathematically, the radome local slope can be expressed as 



R= d0 r /3(A.-0„i). 

The radome characteristics have a predominant effect on missile performance at high 
altitude. This occurs because the radome introduces an angular change between the 
actual LOS and the apparent LOS to the target. The effect of this angular change is an 
erroneous radar-tracking error signal, which will command false missile maneuvering 
and can result in large miss distances (see Figure 3.27). 

Normally, the boresight error is assumed to be negligible as compared to other 
system errors. Also, there is a possible contribution of the refraction error to measure- 
ment noise, when the frequency of the received signal is varied in a pseudorandom 
manner in order to reduce the effect of a potential enemy jammer (for example, when 
the seeker is an active or semiactive radar). This noise can be treated as a contributor 
to range-independent noise. The various noises (to be discussed in Section 3.4.2) 
associated with the seeker must also be considered in the design of a missile. As 
discussed earlier, the LOS angle is the fundamental quantity measured by the seeker. 
These measurements will generally be corrupted by various types of noise, which can 
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Fig. 3.27. Radome error slope ( R ) limitations on maximum launch range. 



be categorized according to the dependency of their rms (root-mean-square) levels 
on the missile-to-target range. The actual noise levels and bandwidths are depen- 
dent on the exact form of the measurement signal processor, target configuration and 
characteristics, environmental conditions, and a number of other system effects. In 
block diagram form, the seeker can be represented as shown in Figure 3.28. It should 
be pointed out that here we consider a single channel; an actual seeker system will 
require the implementation of two or three channels in order to account for motion in 
three dimensions. 

Commonly, the stabilization dynamics comprise the gimbal servo and rate gyro, 
mounted on the antenna. Typically, the stabilization dynamics have a very wide 
bandwidth, in excess of 100 rad/sec. Moreover, the track loop model can be 
represented by simple first-order dynamics, commanding a gimbal rate proportional 
to the measured boresight error. In essence, the loop attempts to drive the boresight 
error to zero, thereby causing the antenna to track the target [5]. 

Using classical feedback control theory, it is easily shown that the linear transfer 
function from dk/dt to s' (assuming unity gain for the radome, the signal processor, 
and the stabilization dynamics) is given by the following first-order transfer function 
[3], [12]: 



s' x 
k (l+jr) 



(3.81) 



Note that at low frequencies (i.e., a> < 1/r), the indicated boresight error is propor- 
tional to the LOS rate. That is, the indicated boresight error s' is scaled by 1/r, 
which forms the desired rate command to the stabilization loop. As we shall see 
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X = Line-of-sight angle 
£' = Measured boresight error angle 
S), = Seeker gimbal angle rate relative to missile 
6 ln = Missile attitude angle rate 
X = Track loop time constant 

Fig. 3.28. Typical block diagram of a seeker subsystem. 



later, this is the desired measurement for classical proportional navigation guidance, 
which commands a missile’s lateral acceleration proportional to the LOS rate. 
Equation (3.81) provides an indication of the important region of boresight error 
linearity. Now, using the fact that s' is proportional to dX/dt in the steady state for 
constant dX/dt, we obtain the following expression for e ' max : 



p — T 
° max L 




(3.82) 



Assuming that r is sufficiently small, s' max can be held within the linear range of the 
received beam width. Figure 3.29 illustrates the resulting seeker block diagram with 
a linear refraction error model. 

From Figure 3.29 it can be seen that the transfer function relating d m to X' is given 
by [3] 



X'/0 m = -R/(l + s r), 



(3.83) 



where X' is the measured LOS angle. Thus, the measured LOS rate is corrupted by a 
term proportional to the body rate. Furthermore, since body rate is a result of com- 
manded acceleration, a loop is formed that can have a destabilizing effect on missile 
attitude, resulting in an increase in miss distance. When R is zero, the contributions 
from the body angular rate input (see Figure 3.29) cancel, producing no effect on s'. It 
is well known that most missiles use some form of proportional navigation as the guid- 
ance law. Although classical proportional navigation guidance uses measurements of 
LOS rate, it is more convenient to use measurements of LOS angle in guidance laws 
that utilize a Kalman filter. In such a case, let us define the measured LOS angle X' as 
follows (see (3.78)): 



X' = {\ + R)X-RQ m . 



(3.84) 
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Fig. 3.29. Block diagram of the seeker model with track loop. 



From Figure 3.29 it follows that 

e' = tsX.'/(1+st). (3.85) 

Now, since the boresight error is an observable quantity, (3.85) can be inverted, 
yielding 

y = [(l + ST )/jr]e'. (3.86) 

That is, A/ can be recovered from the measured seeker boresight error. We can sum- 
marize the above discussion by noting that a typical classical seeker model consists 
of (1) the antenna pointing dynamics, and (2) the parasitic coupling of the missile 
airframe dynamics into the LOS direction as perceived by the seeker. 

Earlier in this subsection we discussed the refraction error and the various errors 
caused by the radome’s nonlinear nature. Furthermore, we noticed that a radome 
(or irdome in the case of infrared seekers) is required in order to protect the seeker 
sensor and to transmit the reflected radar (or infrared energy, as the case may be) 
energy from the target. Regardless of the nature of the seeker sensor used, the 
requirements for the protective dome that the guidance designer must consider are as 
follows: 

1. It must transmit the energy with a minimum loss. 

2. It must have minimum aerodynamic drag. 

3. It must transmit the energy with minimum distortion. Specifically, a change of 
angular distortion with seeker position causes a severe guidance problem with the 
parasitic attitude loop. 

4. It must have satisfactory mechanical properties, such as (a) sufficient strength, 
(b) resistance to thermal shock (e.g., from rapid aerodynamic heating), (c) resis- 
tance to rain erosion at high speeds, and (d) minimum water absorption. 

Figure 3.30 illustrates three conceivable shapes for a radome. 
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1) Ideal aerodynamic with (L/D) — 5. 




2) Ideal electromagnetically with (L/D) = 1/2. 




3) Compromise radome with (L/D) = 3. 

Fig. 3.30. Three conceivable radome shapes. 

For minimum angular distortion, a hemispherical shape (e.g., a hyperhemispher- 
ical shape as in ground-based radar) would be ideal from an electromagnetic design 
point of view. However, the drag of such a shape would be excessive. The aerodynam- 
icist, on the other hand, would prefer the first shape shown in Figure 3.30, because 
this shape minimizes drag. This design shape tends to have rather high peak values 
of radome error slope R. A typical compromise is the so-called tangent-ogive shape 
with a length-to-diameter (L/D) ratio of about 3. Some missiles use much blunter 
domes despite the drag penalty. 

The modeling, evaluation, and compensation for dome error angle effects are 
among the most difficult problems of the guidance designer. For instance, each radome 
from a production run has a different characteristic, which varies with plane of exami- 
nation (defined by the longitudinal axis and the seeker boresight axis), frequency, and 
possibly environmental factors. Preliminary analytical models utilize fixed values of 
radome error slope R , which may be positive or negative. These slopes usually lie 
within the range from —0.1 to +0.1 degree/degree. 

In addition to the conventional and electronically steered (or scanned) array (ESA) 
tracking radar seekers, the following seekers are used extensively in various guided 
missiles (see also Section 3.4.4): 

• Electro-Optical (EO) TV. 

• LADAR (EAser Detecting And /hanging, or Laser Radar). 

• SAR (Synthetic Aperture Radar), as well as Semiactive SAR. 

• Dual Mode LADAR/MMW (Millimeter Wave). 

• HR (Imaging Infrared); also seen as I 2 R. 

• TWS (Track- While Scan) multiple-target tracking radar. 

Conformal RF antennas built into missile nosecones will give better combined 
RF /Electro-optical performance, due in part to removing the mechanical antenna 
gimbal. Next-generation missile seekers will most likely use ESAs. However, elec- 
tronically steered antennas will see application if the cost becomes low enough. More 
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efficient solid-state transmitters, low-probability-of-intercept waveforms, and shorter 
wavelengths such as the 94-130 GHz W -band are other potential RF changes. In addi- 
tion, future seekers using these new technologies will offer high-resolution air-to-air 
and air-to-ground modes, outstanding tracking performance, sophisticated ECCM, 
and high reliability. Visual and infrared seekers will be upgraded to true imaging as 
focal plane arrays and processing chips become cheaper. Improvements over pseudo- 
imaging include better target recognition and clutter/countermeasures rejection. 
Automatic target recognition ( ATR ) may change the rules of engagement for BVR 
air-to-air shots. Instead of requiring positive identification before launch, good ATR 
would make th e friend-or-foe decision itself. 

New seeker technology is evolving with a new generation of advanced missile 
systems. For example, the AIM-9X Sidewinder II design features a high-resolution 
rotate-to-view seeker. In this design, the outer seeker casing, slightly larger in diameter 
than the missile body, rotates 360° to provide a clear viewing path for the seeker, which 
is mounted on a two-axis gimbal fixed to the body. The seeker head can be slaved 
to a helmet-mounted sight, with the seeker window rotating to view what the pilot is 
looking at. It has an off-boresight capability of more than 90°. 

We conclude this section by noting that a new generation of radar is the AESA 
(active electronically scanned arrays). This radar allows tracking of fast, stealthy, 
cruise-missile-size flying targets at hundreds of miles. The AESA will be used to 
upgrade the Joint-STARS airborne ground-surveillance radar. An AESA radar uses 
hundreds or even thousands of small transmitter-receiver (T/R) elements (or modules) 
that allow it to conduct widely diversified tasks simultaneously, including surveillance, 
communications, and jamming. These T/R elements are used to update a radar’s 
computer several times a second, so that target data are much more accurate. For 
example, each array in the AESA radar is made up of about 1 ,000 transmitter-receiver 
elements on the F/A-22 Raptor and several hundred in the JSF. The number of ele- 
ments dictates the power output and range of the radar, which equates to about 125 
miles on the F/A-22 and 90 miles on the JSF. Other aircraft that will use the AESA 
technology are the F/A-18E/F and F-15C interceptors. Some of the new combined- 
technology radars may first be used operationally in the Global Hawk UAV. So far. 
Global Hawk’s, mission has been that of supplying E/O and IR digital photos and 
SAR/GMTI data of vehicles (see also Appendix F). However, HAWK’s new payload 
will not be a fully developed new-technology radar system, but rather an existing 
air-to-air radar scaled down with AESA technology. In some applications the ASARS-2 
(Advanced Synthetic Aperture Radar System) is used to track and observe ground 
targets. ASARS-2 offers a resolution of 1 ft over a 1 sq mile FOV from a range of more 
than 108 nm (200 km) and an altitude of more than 65,000 ft when observing ground 
targets. 

Recently, the Navy tested a combination millimeter-wave radar/RF-homing 
seeker for a follow-on to the current HARM radar- killing missile called the Advanced 
Anti-Radiation Guided Missile. In addition, the USAF is investigating the possibility 
of developing a next-generation HARM- like missile that fits in the F/A-22’ s weapons 
bays. 
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3.4.2 Missile Noise Inputs 

We have heretofore assumed that the LOS rate, dk/dt, can be accurately measured 
by the radar on board the missile. In fact, however, noise within the radar and on the 
target signal limits the accuracy to which dk/dt can be measured, and significantly 
affects the miss distance. In either radar or an infrared system, noise tends to be a 
problem because it increases miss distance. In a radar system, the guidance designer 
would like to have high illuminating power on the target, in order to reduce the 
receiver gain and internally generated noise. The antenna size is made a maximum 
within the constraint of missile body diameter, so as to maximize power reception and 
to minimize the angular beamwidth. Furthermore, a passive infrared seeker is also 
designed to have a maximum aperture to maximize incoming power, and may utilize 
special optical modulation of the incoming radiation in order to amplify without drift 
the weak electrical signal of the infrared detector. Thus, in order to maximize the 
incoming infrared power, a tail attack against a jet aircraft tail pipe is preferred, as are 
moderate ranges and good weather. A small wavelength in a radar reduces the angular 
beamwidth, but the choice of wavelength is limited by problems of power generation 
and environmental absorption, etc. Consequently, a missile radar antenna usually 
has a relatively broad beamwidth, and so it is unable to resolve two closely spaced 
targets by their angular separation until the last moments of intercept. This classic 
problem may lead to a bad miss distance. Therefore, because of its much higher ratio 
of aperture diameter to wavelength, an infrared seeker has a narrower beamwidth and 
much higher angular resolution. Radar illumination may be continuous-wave (C W ) 
or pulsed, depending on which factor of application is governing. In a simple CW 
radar, the closing velocity V c is obtained from Doppler measurements. 

The fundamental effect of noise is to mask or hide the true value of dk/dt. Noise 
can occur due to target effects or receiver (missile) effects. As we saw in the previous 
subsection, the radome contributes a bias error due to the diffraction effects of the 
radome, which is called boresight error. Receiver noise is generated within the missile 
receiver, and the target signal must compete with it. This noise has increased angular 
amplitude at longer ranges, where the signal-to-noise ratio* (SNR) of the target is the 
lowest. Receiver noise consists primarily of thermal noise generated by the antenna 
and receiver electronics on board the missile. The effective amplitude of this noise 
increases with increasing range, because of the corresponding decreasing SNR. There 
are, in general, three types of missile receivers that can be considered. These are: 

Passive: The target supplies the radiated signal. 

Semiactive: The target is illuminated by a source that is not on board the missile. 
Active: The target is illuminated by a source on board the missile. 

Specifically, the receiver will generally include some type of automatic gain con- 
trol, which attempts to keep the receiver signal power nearly constant. As a result, the 

*The SNR is defined as 10 logio (A~/o~), where A is the signal amplitude and is the 

variance of the noise. 
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effective noise level will change with received signal power relative to some reference 
level. Commonly, a normalized angular measurement noise model is defined that uses 
the variance (or power spectral density) of the indicated boresight error, measured at 
a range that yields an SNR of unity as the reference level. 

In addition to the receiver noise, there are four other basic types of noise associated 
with the design of a missile guidance system. These are: 

Glint: Glint, or angle noise, is the phenomenon in which interference by two 
or more sources causes a distortion in the shape of the propagating wavefront. The 
effect of this distortion is to change the apparent angle of arrival of the wavefront. This 
appears to the tracking system as a wander of the apparent target location from its true 
location. Thus, glint is a target-induced error term that introduces an angular error 
in the target tracker. The apparent center of the target moves along the length of the 
target and can occasionally exceed the target dimension. Furthermore, the apparent 
location of the target may lie outside of the target a significant amount of time. (Note 
that the phenomenon of glint is also known as the radar bright spot wander.) Since 
glint is a distance error along the target, the equivalent angular error varies as 1 / R mt , 
where R mt is the total missile-to-target range. Glint is usually described as a Gaussian 
random variable whose main value is at the center of the target and whose standard 
deviation (er) depends on the target span, perpendicular to the LOS angle. A typical 
value for the standard deviation of correlated glint for an aircraft is 

cr = 0.25 S/R mt , (3.87a) 

where S is a characteristic length (or effective target length). The correlation coeffi- 
cient is computed by 

p = exp(—a> g ■ AT), (3.87b) 

where 



Wg = the glint half-power frequency, 

AT — magnitude of time since last call. 

The standard deviation of the correlated glint error is then computed by the following 
relation: 

cr c — a ■ 0 ~ p 2 ) 1/2 - (3.87c) 

The spectral density of glint error is maximum at zero hertz and decreases with 
upward concavity as frequency increases. The glint spectral density is commonly 
approximated with a Lorentzian lineshape as follows: 

Og(&>) = Oo[<Ug/(tt> 2 + (Og)], — oo<&xoo, (3.87d) 

where 



<I> ? (ft>) = spectral density, 

<l>o = zero-frequency value of the spectral density, 
a>g — glint half-power frequency. 
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The glint half-power frequency is given by 

o) g = 4 jt£2S/X, 



where 



co g = glint half-power frequency, 

Q = rotation rate of target, 

S = target characteristic length (or span), 
X — radar wavelength. 



The variance of the glint error is given by the integral of the spectral density (3.87d). 
Thus, 

O ' 2 — j Og(&))t/(W, 

which yields, after defining z = a>/a>g and dz = dco/atg, 

a g 2 = f dz/(l + z 2 ) = 7T®ocog. (3.87e) 



A frequency analysis of the time record of glint from an aircraft target often sug- 
gests that its spectrum approximates white noise passed through a first-order lag of 
the form [3] 

<f> ? = (Kg) /(l + of 7) 2 ) in units of m 2 /rad/sec, (3.87f) 

where 71, is the guidance time constant (typically in the range 0. 1 to 0.25 sec), and 
L 2 is the mean-square value of the glint and is given by A 2 = tt A" -2 / 2 77, (if 71, and 
Lg are known approximately, then A" 2 can be evaluated). 

Range-Independent Noise (R1N): This noise, also known as fading noise, has 
a constant angular amplitude, and is caused by amplitude fluctuations of the target 
occurring at the information frequency in the missile receiver, for example, at the 
conical scan frequency of a conical scanning missile seeker. Range-independent noise 
ay is inherent in the missile receiver. The noise can be modeled as ay ~ N( 0, ry ) , 
that is, zero mean and variance, with equivalent white noise power spectral density 
O given by 

<J> = 2 t/rf 

where if is the correlation time constant and ;y is the variance. Also, RIN may be 
caused by the seeker servo system. (Another type of noise commonly encountered in 
connection with a missile radar seeker noise is range-dependent noise. However, this 
is strictly a ground-tracking radar noise, used in command guidance systems.) 

Scintillation Noise: Scintillation is a phenomenon similar to glint, in that reflec- 
tions from various parts of the aircraft (e.g., target) interfere. In the case of scintillation, 
this affects the amplitude of the received signal. Typically, an aircraft is composed 
of many individual conducting surfaces, or scatterers, each with different scattering 
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properties that vary as the viewing or striking angle changes. In addition, multiple or 
sequential reflections of the radiated signal may occur between the various scatterers. 
These features can strongly affect the resultant value of the target aircraft’s radar 
cross section ( RCS )* signature. The large-amplitude fluctuation of an RCS signature 
with respect to small changes in the viewing angle is referred to as scintillation. 
In general, the amount of scintillation decreases as the wavelength increases. The 
presence of interference in the returned signal causes a modulation of the returned 
signal amplitude. In conical scan systems, any component of amplitude scintillation 
at or near the system scan frequency will be interpreted by the system as a signal 
resulting from an off-axis target. For this reason, frequency components of amplitude 
scintillation near the scan frequency will drive the radar off target. Note, however, 
that this effect is absent from monopulse radar systems, which do not extract infor- 
mation from modulation frequencies. Scintillation errors can be modeled in the same 
manner as glint errors, with the exception of the form of the standard deviation. For 
scintillation, the standard deviation is given by 



where 



a, = 



ByJW {u>g)B n 
E 



(3.87g) 



B = beamwidth, 

B n = equivalent bandwidth for the noise of the tracking loop. 



E — error slope, 

1 Wf, 

W(co g ) = 



2tx an + an 



In the above equation, a> g is the glint half-power frequency, and oj s is the scan 
frequency. As for the glint error, the correlation coefficient is given by 



p — exp (— cog • A T), 

and the standard deviation of the correlated output by the relation 

er c = cr(l — p 2 ) 1/2 . 



The reflection characteristics of an aircraft determine both the RCS level and the 
amount of RCS scintillation and target glint, all of which affect tracking accuracy. 
The angular scintillation noise can also be expressed as 

o? = °ld/* 2 m,’ (3.87h) 



*The RCS , a , is a measure of the size of the body as seen by the radar. The RCS is an 
area and is usually measured in square meters or decibels, with 1 nr reference level, and, 
except for the sphere, is aspect-angle dependent. Specifically, a is 4 jt times the ratio of the 
signal power per unit solid angle (i.e., one unit solid angle is the s teradian, and there are 47r 
steradians in a sphere) scattered in the direction of the receiver to the signal per unit area 
(the signal power density) that strikes the body. 
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where er 2 rf is the variance of the apparent wander distance and R mt is defined as 
before. This equation states that although the magnitude of the wander is essentially 
range independent, the equivalent noise on the LOS angle measurements increases as 
range decreases. 

As previously noted, receiver- and range-independent noise are generally assumed 
to be wide-band relative to the guidance system bandwidth. Angular scintillation noise 
is, in general, a narrow-band source and is often modeled as white noise through 
a low-pass filter with a time constant that depends primarily on the target motion 
spectrum. 

Scintillation can be an important factor in various parts of an engagement. For 
instance, during acquisition, target fades, or periods of low aircraft return, can inhibit 
detection of the aircraft and therefore cause a lengthening of the prelaunch time of an 
engagement. Similarly, a surge in return strength can make an otherwise undetectable 
aircraft visible to a radar. Scintillation is usually less important once radar tracking 
has been established, since a lower signal level is needed to maintain track. However, 
in some cases, particularly in marginal detection circumstances, an aircraft fade can 
cause a radar loss of track. 

In the application of modern optimal control and estimation theory to modeling 
of the seeker and missile/target dynamics, glint and scintillation errors are com- 
monly modeled using filtered Gaussian white noise input in order to produce corre- 
lated noise output. The new noise treatment replaces the fourth-order Runge-Kutta 
approximation with the method of conditional probability density function {pdf). 
This new technique allows random draws of correlated errors to be made directly, 
thus eliminating the need to make white noise input draws and to filter this input 
before output draws are made. The probability density function of the correlated 
output error is found in closed form in terms of the previous value of the correlated 
output and the correlation coefficient. The new method provides accurate statistics and 
satisfies the necessary correlation properties. Its computational simplicity translates 
into substantial savings in computer processing time. The correlated output terms for 
glint and scintillation are computed using the same form of the conditional probability 
density function. The pdf is derived from the spectral density and autocorrelation 
function and is given by the expression 



P(gl\gi) = 



sj2n{\- p 2 )- 



• exp 



2(7,? • (1 -p 2 ) 



(gl-pglf 



(3.88) 



where 



g i = previous value of the error term, 
g 2 — current value of the error term, 
a c — standard deviation of correlated noise, 
p — correlation coefficient. 

This is a Gaussian density with mean p and variance er 2 (1 — p 2 ). Error terms are 
computed using the radar error covariance function RADHVfoy), which returns a 
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normally distributed random number with zero mean and standard deviation a c . The 
expression for computing the glint and scintillation is 

g2 — p-g\+ RADEV (cTc), (3.89) 

where appropriate statistics are substituted into the above expression for glint or 
scintillation. 

Thermal Noise: Thermal noise appears as random signals within the system band- 
pass. Thermal noise is usually viewed as white Gaussian noise, that is, noise with a 
power density equal across all frequencies and with an instantaneous value given by 
a Gaussian probability density function. The signal energy the radar tries to track is 
contained within the pulse repetition interval (PRI) lines within the system bandpass. 
Thermal noise contained within the bandpass is also passed through the signal pro- 
cessing elements and competes with the true signal. Since the thermal noise voltage 
is a random process, it tends to drive the radar off the true target return. The direction 
of this thermal-noise-induced track error changes continuously with time, since the 
noise instantaneous value is continuously changing. If the target return is much larger 
than the thermal noise return (e.g., if the signal lines within the bandpass are much 
larger than the noise level), then the thermal noise will have a relatively small impact, 
and system tracking will not be significantly disturbed. However, if the target return 
is small or if the noise level is large, then the energy contained in the noise can swamp 
the energy return of the target, and the system will drift far off the true target position. 

The total measurement noise variance is the sum of the variances of the individual 
uncorrelated noise components 

a ? = a2 f+ a c+°l+°?h’ (3.90) 

where aj h is the thermal noise variance. 

Other typical error sources of a guidance system are multipath and clutter effects. 
Multipath and clutter occur naturally in the low-angle track situation, and each is 
capable of degrading radar-tracking performance. More specifically, multipath and 
clutter are types of noise signals caused by reflections from terrain surface features. 
Multipath and clutter effects will tend to degrade radar performance is, for example, 
tracking low-altitude targets. Multipath is a result of multiple paths the radar signal 
makes from the radar site to the target and return. Both the specular multipath, which is 
that governed by Snell’s law-type reflections from a flat surface, and diffuse multipath 
(or random scatter from rough surfaces) reflection components are considered in radar 
tracking error analysis. The apparent range caused by multipath effects is 

Range= (|Rrsl + IRrspI + |RspsI)/2, 

where Rjs is the target-to-site vector. R?s p is the target-to-specular-point vector, and 
Rsps is the specular-point-to-site vector. The diffuse multipath, which is a random 
scattering of the radar energy from rough surfaces, can be implemented using Monte 
Carlo techniques. Clutter, which is the radar energy return that has been backscattered 
from the terrain surrounding the target, provides a competition signal to the target 




